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Summary. ‘The non-linear differential equation for the velocity potential of three- 
dimensional steady irrotational supersonic flow past wings of finite span has been in- 
vestigated. It is found that the whole Mach number range from 1 to ~ may be divided 
into two regions (not strictly divided), in each of which similarity laws are obtained, 
with two parameters K, = (M* — 1)'”/r" and K, = A(M’ — 1)”; + is the non-di- 
mensional thickness ratio, A the aspect ratio of the wing, M the Mach number of the 
uniform stream in which the wing is placed. The factor n is given explicitly as a function 
of M and 7; in the lower region of Mach numbers it tends to 1/3 as M — 1, for all r, 
giving the ordinary transonic rule, and in the upper region it tends to —1 as M >, 
for all 7, as in the ordinary hypersonic rule. 

It is shown that both two-dimensional flow and flow over a three-dimensional slender 
body, including axially symmetrical flow, are special cases of the present analysis, 
involving only one parameter K, in the similarity rules. 

I. Introduction. One of the major difficulties in solving problems of flow of com- 
pressible fluid is the non-linearity of the differential equations that govern the flow. 
Over a certain range of Mach number, the differential equations of flow of compressible 
fluid can be linearized for many investigations of practical importance, and the resulting 
equations give valuable information on the flow field of a compressible fluid. However, 
there exist other ranges of Mach numbers, in particular the transonic flow (M = 1) and 
the hypersonic flow (M > 1), where the differential equations cannot be linearized. In 
these cases, we have to study non-linear equations. At the present time, there is no 
general method of solving these non-linear differential equations. For practical purposes 
the study of the effects of non-linearity on the flow must rely mainly on experimental 
results. In order to increase the value of either theoretical or experimental data, it is 
useful to find some similarity laws which will enable a family of possible solutions to be 
deduced from a single result. 

In general similarity laws may be found for bodies moving in the compressible fluid, 
which have at least one dimension perpendicular to the direction of main flow small in 
comparison with that in the direction of main flow. For the linearized theory well-known 
similarity laws have been obtained by Glauert [1] and Prandtl [2] for subsonic flow, 
and by Ackeret [3] for supersonic flow. However, it will be shown later that the similarity 
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laws for linearized theory are arbitrary. This fact has already been observed in the study 
of linearized axially symmetrical flow of compressible fluid [4, 5], where contradictory 
results have been obtained by the use of different similarity laws. Unique similarity laws 
can be determined only when non-linear terms are considered. 

Von Karman [6] was the first to obtain the similarity laws for transonic flows, where 
the fluid velocity is very near to the velocity of sound, for both the two-dimensional 
and the axially symmetrical cases. Spreiter [7] extended the transonic similarity laws 
to wings of finite span. 

The similarity laws for hypersonic flow, where the fluid velocity is much larger than 
the local velocity of sound, were first obtained by Tsien [8], again both for two dimensions 
and axial symmetry. Hayes [9] extended the hypersonic similarity laws to three-di- 
mensional slender bodies. 

Spreiter [7] tried to combine the similarity laws of transonic flow with those of 
linearized theory. He was able to do so because of the arbitrariness of the similarity 
laws of linearized theory. Van Dyke [10] empirically generalized the hypersonic similarity 
laws to the supersonic flow where the fluid velocity is not much larger than the local 
velocity of sound. It is the object of this paper to investigate the possibility of determining 
uniquely similarity laws which link flows over more widely extended ranges of Mach 
number than are covered by the particular cases cited above. 

We begin with the non-linear differential equation for the velocity potential of three 
dimensional steady irrotational supersonic flow for wings of finite span; two-dimensional 
and axially symmetrical flows, and flow past three-dimensional slender bodies, may be 
considered as special cases of this general flow. The relative importance of the various 
non-linear terms is discussed so that the ordinary transonic and hypersonic similarity 
laws are generalized for large Mach number range. It is found that from a known solution 
for a given Mach number and ratio of thickness to chord, a family of solutions may be 
determined; there is however a barrier dividing the Mach number range into two regions, 
so that a solution which falls into the ‘generalized transonic range’ cannot be used to 
give a flow in the ‘generalized hypersonic range’, and vice-versa. 

2. Fundamental equations and boundary conditions. If ® is the velocity potential 
of a three-dimensional steady irrotational flow of compressible fluid, the differential 


equation for @ is 


(a? — &)%,, + (a? — &)4,, + (a? — 0°(6,, — 26,9,6,, 


— 26,6.6,, — 26,6,6,, = 0, (1) 


where x, 7, z are the Cartesian coordinates, subscripts denote partial derivatives, i.e.; 
@, = 00/02 etc., and a is the local sound speed which is determined by the equation: 


i a oe ' , 
_ +t y 5) (db: + P ++ P°) = do. (2) 


Here y is the ratio of the specific heats and a, is the speed of sound in the gas at rest. 
Now if a thin wing of finite span is placed in an otherwise uniform stream of velocity 
V in the z-direction, we may introduce a perturbed velocity potential ¢ such that 


@p = Viz + od) (3) 
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Op O 0 
ee SK. 
Ox OY Oz 
Equation (2) may be written as follows: 
a” y ell 279 2 2 2 
a’ = I sarnd 2 M (2¢, + og; + dy, + ¢:), (4) 
1 
where a, is the speed of sound corresponding to the free stream velocity V and M = V/a,, 
the Mach number in the undisturbed stream. 
Substituting equations (3) and (4) into (1) and retaining terms up to second order, 
one has, with \* = (y — 1) (y + 1), 


\1 — M? — (y + 1)M’9, —- oe M*(¢, + 6h. 


bo|-+ 


4 ‘1 — (y anne 1)M’¢, a * : M’*(¢; + N'6)o 


No) 


she {1 — (y — 1)M’¢, — YT M*('¢2 + obs. 


— 2M*$,¢., — 2M’$.4.. = 0. (5) 


This is the fundamental differential equation for ¢ for the consideration of similarity 
laws, and is a generalization of Tsien’s equation for two-dimensional flow [8]. 

In (5) we have retained terms which, while not important everywhere, may be so in 
certain regions. For example, in hypersonic flow ¢; , ¢: , ¢, are all of the same order 
of magnitude. Furthermore, while we retain in (5) terms which are usually second-order, 
we shall seek and use those terms which give an effectively first-order equation. 


The first-order boundary conditions for ¢ are: 


(1) at infinity 


$. = b, =o. = 0; (6) 
(2) on the surface of the wing, which is represented by 
z = A(z, y), (7) 
the normal velocity component is zero, i.e. 
ae) dh(x, y) 
= 8 
(38 _ Ox (8) 


3. General discussion of similarity rules. We are going to find conditions under 
which it is possible to reduce the differential equation for ¢ and the boundary conditions 


simultaneously into non-dimensional form. We write 
~ dm nm _— a , 
oo fc, a nb, >= gl, ¢ a ¢ m, (9) 


where c is the mean chord of the wing and b its span; thus, A = b/c = aspect ratio. 
l and m are conversion factors which are to be determined. 
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Substituting equation (9) into (6) and (8), we have, 





¢ =o =o =0 at infinity, (10) 
a9) a 
= = . £ 

(22 <0 m Je, n); (11) 


where T f(t, 7) = A(z, y) and T = rc is the maximum thickness of the airfoil section; 
it will be supposed that 7 < 1. 
Substituting (9) into (5), we have 


‘y (+o (Z)() +> (2) ()(% “()') | 
| car n(! + + DM ™)(=) 6 +S (™) (2) (Ss + oF(5) ) oe 
: yr? @ 4 ym { " . ( l \? 
+|-! ty — I) ” + a oo ae ae . l oF Pun\ 4 


: 2 72 2 
+ | -1 +o ~ nM 264 w(™) {y e+ (<) of} ls 


c “ c 


+ 


+ au’(™)() oes, + 2m Las = 0. (12) 


c/\cA Cc, 
In order to get the parameters in the similarity laws we write 
l=cr™", m=cr', (13) 


where n and n’ are factors to be determined. 
Then, for similarity laws to be possible, considering only the linear terms in (12), 
the following necessary conditions are found: 


oa = Ki, (14) 
Ar" =K}, (15) 

or, combining these, 
A(M’ — 1)'” = K, (16) 


where K, , Kj and K, are constants. 
From the boundary conditions, we have 
7'-"-" = constant; (17) 
this can only be true for variable 7 provided that 


n+n =1, (18) 


and the value of the constant is unity. 

Equations (14), (16) and (18) represent the conditions for similarity laws in the 
linearized theory of compressible flow. Only the parameter K, for the aspect ratio (16) 
is unique; the index n in the parameter K, , is arbitrary, since the choice of n’ is arbitrary. 
In order to have unique similarity laws it is necessary to study the non-linear terms. 
From (12) we see that the non-linear terms fall into two groups, one group being the 
more important for transonic-supersonic flow (n > 0), the other for supersonic-hyper- 
sonic flow (n < 0). In the transonic-supersonic flow region the important non-linear 
term is (y + 1)M?(m/c)(I/c)*o{ . For this term to be unchanged with variation of airfoil 
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n'—2n 


thickness, for example, in the immediate transonic region (17 = 1) we must have r 


const., i.e., 
n’ — 2n = 0. (19) 


From (18) and (19), we have 
mn = 4- n’ = §. (20) 


This is the well-known transonic similarity law due to von Karman [6]. 
When M is very large, the important non-linear terms are in the second square 
bracket in (12); both are of the same order of magnitude as (y — 1) M’*(m/c) $%; for 


similarity we require 


M’,r" = const. (21) 
Thus, from (14) and (21), 
n’ = —2n. (22) 
Finally, from equations (18) and (22), 
n=-l, n’'=2 (23) 


This is the well-known hypersonic similarity law due to Tsien. 

It is clear, from an examination of (12), that in order to obtain similarity laws which 
are valid in the intermediate (supersonic) range of Mach numbers we must look for a 
variation of n with Mach number. It is found in the next section that n depends upon 
two parameters, M and 7 (say), but that the value of n tends to 1/3 as M — 1 and to 


— las M —o, for all r. 
4. Generalized similarity law for transonic-supersonic flow. As M increases from 
unity, the non-linear term of lowest order is the second term in the first square bracket 
in (12). 
In order to preserve the form of (12) unchanged up to and including this term, when 
the Mach number is not restricted to be in the immediate neighborhood of unity, the 


parameters required are 


= = = Ki, A(M* — 1)" = K,, 


" 
with n + n’ = 1 from the boundary condition as before, and 
M*,;"-* = K;. (24) 
By elimination of 7, we derive 
n{4 log (K;/M) + 3 log {(M? — 1)/Ki}] = log {(M? — 1)/K7}. (25) 


This shows clearly that n — 1/3 as M — 1. Furthermore, by differentiating and pro- 
ceeding to the limit, it may be shown that 


0 (Ks =D), 


Boa | = | - aad (K; > 1), 
\— @ (K; < 1). 
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It is plausible to reject infinite values of dn/dM at M = 1, in view of the established 
position of the transonic similarity law which requires n = 1/3 to be approximately 
valid in the neighborhood of 17 = 1. Then K, = 1, and the law of variation of n with 
M is 


l M* — 1 M* — 1 
nm = oO _— iy - , " 2 ) 
? 3 log ( K )/ log ( Kew" ). (26) 


Examples of the variation of n with M, for three different values of the parameter K, , 


are shown in Figure 1. 
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We are now able to deduce similarity laws in the sense that a solution $(é, n, ¢) of 
the differential equation of motion may be used to give an infinity of flows in the transonic- 
supersonic region. Given values of the parameters K, , K, , if we choose a value for M 
the corresponding values of 7, A, n are found. 

It will be noticed that the solid curves have been stopped at n = 0. As n — 0 the 
relative importance of the non-linear terms in the second square bracket in (12) increases, 
and at n = 0 these terms are of the same order of magnitude as the “transonic”’ term 
which has been the crux of the above discussion. 

5. Generalized similarity law for supersonic-hypersonic flow. When n is negative, 
the most important non-linear term is the second term in the second square bracket in 
(12). The flows for which this term is retained we call ‘generalized hypersonic flows’. 
In order to have coefficients in (12) which are independent of Mach number and thick- 
ness ratio we must have 


M’ 7" = constant = K;i(n’ ~ 2as M >) (27) 


in addition to (14), (16) and (18). 
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Consideration of (27) and (14) as M —o- shows that n — —1 and K, — K, . But 
K, and K, are constants; hence K, = K, . Elimination of 7 yields finally 
log {(M? — 1)/K?7} 
= — —— sa 28 
, log {M*/K\(M* — 1)} (28) 


The variation of n with M is shown for three values of K, in Figure 2. As before the 


reliability of the approximations made falls off as n — 0. 
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6. Conclusions. Both the differential equation for the velocity potential and the 
boundary conditions for the steady irrotational supersonic flow over a thin wing of 
finite aspect ratio A can be put into a non-dimensional form containing two parameters 
only, when only the most important of the non-linear terms are retained. The parameters 
are K, = (M* — 1)'*/r" and K, = A(M* —1)"’, where M is the Mach number, 7 the 
ratio of maximum thickness to chord and n is a function of any two of M, 7 and K, . 
The similarity laws thus hold for such profiles and Mach numbers as correspond to the 
, K, . There are, however, two important restrictions. The most 


same parameters A 
important non-linear term being different in the two cases n > 0, n < 0, the line n = 0 
divides the n,M-plane into two fundamentally different regions; it is therefore not 
possible to use a result which corresponds to positive n to deduce a flow which corresponds 
to negative n, or vice-versa. There is a region around n = 0 in which neglected and in- 
cluded terms are of approximately the same order of magnitude, so that the reliability 
falls off appreciably near n = 0. Further, the use of the rules should not be attempted 
for large variations of Mach number, e.g. from a supersonic flow to a genuinely hyper- 
sonic flow, as the latter requires the inclusion of a term which is neglected in the present 


discussion. 
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For two-dimensional flow A — ©; the similarity laws then contain only one para- 


meter K, . 


For flow over a three-dimensional slender body, including the axially symmetrical 


flow, A — 0; here also the similarity laws contain only the parameter K, . 


In practical cases, K, is usually larger than one for the transonic-supersonic region 


and smaller than one for the supersonic-hypersonic region, when 7 is small. 
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LOWER AND UPPER BOUNDS TO THE ULTIMATE LOADS OF BUCKLED 
REDUNDANT TRUSSES* 


BY 
E. F. MASUR 
Illinois Institute of Technology 


Synopsis. In the buckled state, statically indeterminate, rigid-jointed trusses 
support loads which are generally in excess of those corresponding to initial instability. As 
buckling proceeds, the loads usually approach limiting values, called “ultimate loads’. 
Two theorems are derived establishing lower and upper bounds to the ultimate loads. 
Elastic behavior is assumed throughout. 

Introduction. If the rigid joints of a truss are subjected to a set of external loads, 
whose values are fixed except for a common multiplier A, , the equations of equilibrium 
of the i-th joint are of the type’ 


: feos ¢:; x,\ _ i ae , 
p» Sii lsin $s , Ss = 0 (¢ = 1, 2, n) (la) 
>~M;,=0 (¢=1,2, ---n) (1b) 


where S,;; and M;, are, respectively, the axial force and bending moment in bar (i, j) 
at joint 7, ¢,; is the inclination of bar (z, 7), X; and Y; are external force components, 
and where the summation extends over all the bars adjoining 7. If now all bars are 
composed of an elastic material and the elementary theory of beams in bending under 
axial forces is assumed to hold, the joint moments M;; can be expressed as linear func- 
tions of the joint rotations 6; [2]. Substitution of these expressions into Eqs. (1b) leads 
to a set of linear, homogeneous equations in the joint rotations of the type 


n 


>> a,;0; = 0, (¢ = 1, 2, --+ n) (2) 


l 
where the a;; form a symmetric matrix whose components are transcendental functions 
of the physical characteristics of the bars, and of the axial forces. The equilibrium of 
the structure is stable only if the matrix is positive definite; a necessary condition of 
neutral equilibrium is given by 
f(do) = | ay; | = 0. (3) 

The lowest positive real root of Eq. (3) is the desired “critical’’ multiplier 4) at which 
buckling occurs. 

By definition, the bar forces S;; are uniquely determined by Eqs. (la) for given 
\» in a statically determinate truss. If, on the other hand, the truss be of m-th degree 
of redundancy relative to its axial force distribution, the most general expression for 
the bar forces is of the form 


S.; = 2d A, Si; (4) 
=0 





*Received Nov. 21, 1952. 
1The effect of the linear displacements is ignored in Eqs. (1) and (2). This is in conformity with the 
experimental results given in Ref. [1] (see the bibliography at the end of the paper). 
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386 E. 
in Eqs. (4), the force systems S;;’ are restricted by the set of equilibrium ‘equations 


fimat . . te 
vir)) COS; , A ‘ : ‘ . = 
pW Ser y+ 8.5.7 = 0, (r = 0,1, --+ mt = 1, 2, --- n) (5) 
‘sin @ ee 
where 6,, is the Kronecker Delta. 
It is convenient to restrict the force systems S;;’ further by the requirement (which 


can always be fulfilled) that 


YS Seaae id JFR 
z FE { Zz = poor ; (ry = 0,1, 2, --+ m) (6a) 


D> A Z j k = 5. : (r. s= :: 2. “ee m) (6b) 


where the summation extends over all the members of the truss.” Eqs. (6a) indicate 
that the force system S;;’ is the actual force system, i.e. the one corresponding to mini- 
mum strain energy, in the unprestressed, unbuckled truss for A) = 1. 

The axial forces in the buckled state are further restricted by the “‘compatibility 
condition” that 


| 


b6L, = (u; — u;) cos ¢;; + (v; — v;) sin g;; , (7) 
in which 651, is the change in the distance between joints 7 and j relative to the unloaded 
state and (u, v) are the cartesian components of the linear displacement vectors of the 
joints. If both sides of Eqs. (7) be multiplied by S,”’ (r = 1, 2, --- m) and summed 
up over all the members of the truss, we are led, after rearranging the summations and 
in view of Eqs. (5), to the system of “virtual work’’ equations 


> Sask, = 0. (rr = 1,2, --- m) (8) 


For the sake of generality it will be assumed that the truss has been prestressed, 


with the initial force system expressed by 


St = > v%S,”. (9) 


g=1 
‘ 


On the other hand, 6, can be written in the form 


y —— Yok 
it, oe (10) 


E,A, 
where 5, represents the shortening of the chord length of the k-th member corresponding 
to its curvature in bending. If we substitute Eqs. (10) in (8) and consider Eqs. (6) and 
(9), we arrive at the following set of equations for any bent state: 


r —r* = D> G6, . (ry = 1,2, --- m) (11a) 
k 


It has been shown elsewhere [3] that, for an actual buckling mode, 


je 8,’ 6, = 4 Zz Z. a;;,,0;,0; = wf,- (r = 0, 1, 2, --+: m) (11b) 


‘ 


2S, .L,. Ex, A; are the axial force, length, modulus of elasticity and cross-sectional area, respectively, 


of the k-th bar. 
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in which a subscript, preceded by a comma, designates the partial derivative with 
respect to the corresponding parameter \, or f,, = 0f/dA, . Thus, during buckling, 

A, — A¥ = uf, , (r = 1, 2, --- m) (12a) 
where » > 0 is a measure of the extent to which buckling has proceeded and is expressed 
by 


1 2% 


~<a (12b) 


24, 


in which A,; is the cofactor of a;; in Eq. (3). It follows immediately from Eqs. (3) and 
(12) that, as buckling commences, the load \, will generally increase, and never decrease. 


Considering all parameters A, (r = 0, 1, 2, --- m) to be functions of u, Eq. (3) can be 


Mm 


written in the form 


F(u) = f[ro(u), Ar(u), --* An(w)] = O (13) 
for all values of uw. Differentiating Eq. (13) with respect to u, we obtain the identity 
7 .&, B. &~ 
du J.0 du + ds Sut du 0. (14) 


On the other hand, by dividing Eq. (12a) by u and letting » approach zero, we are led to 


Ir 
- (u=0)=f,¢=0) (r=1,2,--- m) (15) 
Cy 
It follows, from Eqs. (14) and (15), that as the truss is on the verge of buckling 
Xo ete Se ae 
— = -—(f,.)"' > (f.)’, (16) 
au f=1 


where all the terms in Eq. (16) are to be evaluated at u = 0. 

It can be shown that f,» is negative. In fact, if \, be the smallest positive buckling 
parameter, i.e. the one corresponding to neutral equilibrium, the positive definiteness 
of the matrix [a;;] implies f(A.) > O for X» < Aj . Since, by Eq. (3), f(A) = 0, and if 
multiple roots be excluded, it follows that 

f o(Ao = Xo) < 0. (17) 
An inspection of Eqs. (16) and (17) shows that 


> (u = 0) > 0, (18)° 
where the equality sign applies only to the special case 
f.(A = »*) = 0. (r = 1, 2, --+ m) 
It follows further from Eqs. (12) that for increasing y, and provided the X, all stay 
Actually, this relationship applies to all values of yu, as can be demonstrated by means of energy 


considerations. 
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finite,“ the external and internal forces approach a state governed by the equations 
SOM =O 6 = 1,3,--- | 


in which the superscript designates limiting, or “ultimate” values of the parameters, 
with \> referred to as the “ultimate load’’. 

The determination of A} is of interest to the engineer since it represents the ultimate 
carrying capacity of the truss. Since, however, the solution of Eqs. (19) presents formi- 
dable numerical obstacles, it is easier to estimate the value of \; by means of two theorems 
establishing lower and upper bounds to the ultimate load. 

A stability criterion. Before proceeding to these two theorems, it is convenient to 
establish a stability criterion for the unbuckled truss, which is to be used in the later 
proof. This criterion is based on the assumption that, for a given system of external 
forces identified by the load parameter \, , the potential energy corresponding to the 
unbuckled state is less than that of any neighboring, ‘‘geometrically consistent” bent 
state if the equilibrium of the truss is to be stable, or 

Y’ < fF, (20) 


In what follows, a geometrically consistent bent state is defined as one in which 
the deflection curves of the bars are of sufficient degree of smoothness to make the 
discussion meaningful, and satisfy the geometric boundary conditions of continuity 
at the joints. Geometric consistency also implies the satisfaction of Eqs. (7) and (10); 
it is finally assumed that the force equations (4), (5), and (6), and therefore also Eqs. 


(8) and (lla), are valid. 
With these definitions, the potential energy 


vol=¥ es Seo. ~ xu Sika + Fad, 
2 “7 EA, r 


in which U; is the strain energy associated with the bending of the truss bars, can be 
expressed by Eqs. (4), (5), (6) and (7) in the form 


a Lie J, = 
V= grits LN+ Un— 2d Sil ; 


in view of Eqs. (8), (9), and (10), this reduces to 


: ] 1 \> . 1(0) ¢ 
V=—gPmts LMt Usted Sh. 
Similarly, the potential energy in the unbent state is given by 
fa «tg g ? For 
=s 9 pi 0 2 a Ag - 


Thus, criterion (20) now takes the form 
: , l< m ‘ . (0) 
Po 7 ag D> (Ae — A*(Ae + AN) + Us +A D SiS, > O 
==1 


‘For a discussion of this question see Ref. [3]. 
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which, by Eqs. (lla) and by considering only states in the neighborhood of the un- 
buckled state, is finally expressed as 


V—V*= > S65. + Us > 0, (21) 

k 
where the equality sign applies to the case of neutral equilibrium, with the bent state 
represented by the first buckling mode. The stability criterion (21) can also be demon- 
strated directly. In fact, if the arbitrary set of deflection curves y,(x) be developed in 


a Fourier series 
ydx) = > ec, ni(2). 
T=1 


in which 7.(x) represents the normalized buckling mode associated with the positive 
load parameter Xj , it can be shown, by a series of partial integrations and by virtue of 
the orthogonality of the normal modes, that 


Zz, Side + U, = p c7(X5 a Xo). 
; t=] 


The right side of this equation is positive definite if \y < Xo, ie. when the load parameter 
is smaller than the smallest buckling load. 

Lower bounds to the ultimate load. After these preliminary remarks, we now pro- 
ceed to state the first theorem. In what follows, we shall call a load parameter A, > 0 
“critical” if there exists a set of prestressing parameters A* (r = 1, 2 --- m) such that 
the truss is in neutral equilibrium with its force system defined by (Ao , A%_, «++ AS). 
With this definition, which is analogous to that of a “‘statically admissible multiplier” 
in rigid-plastic limit design [4], we state 


THEOREM I: The ultimate load is the largest of all critical loads. 
From this theorem we deduce the following 


Coro.uary: If an arbitrary internal force system, and an external load system iden- 
tified by Xj , satisfy the force equations of equilibrium, and if the truss, with its forces so 
defined, is in stable equilibrium, then 


Mo > rA- (22) 

It is clear that this corollary establishes an easily calculable lower bound to the ultimate 

load since the stability of a truss for a given force system can be determined by a number 
of stability criteria [5, 6). 

Proor: It can readily be verified that if the prestressing parameters be so selected 


that AX* = A! (r = 1, 2 --+ m), Eqs. (3) and (12a) are identically satisfied by 


1 =r» (r = 0, 1, 2, --+ m) 


for all values of u. In other words, in this “special case’’, the truss buckles similarly 
° . . 5 
to a statically determinate truss under constant external loads and internal forces. 





‘It has been shown previously [7] that this case corresponds to a stationary value of the buckling 


load. See also Eq. (18). 
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Inspection of Eq. (1la) shows that the buckling mode of the special case is such as to 
satisfy 


> Ss =0. & = 1,2, --- m) (23) 


If now the stability criterion (21) be applied to a truss in stable equilibrium under an 
external load 4 , and if the arbitrary deflection curves be so selected as to coincide 
with the buckling mode of the “‘special case’’, it follows from Eqs. (23) that 


V—Vt*=nr > Sh. + Us > O. 


Since, also by Eqs. (21) and (23), 


> Sd, + Uz = 0, 


and since further U’, is positive definite, it follows that 
Ao > Ao. (22) 


Upper bounds to the ultimate load. Before proceeding to the discussion of upper 
bounds to the ultimate load, it is useful to state the following 

Lemma I: For a given set of prestressing parameters, the value of the critical load 
of a truss is not reduced by increasing the moment of inertia of a bar.° 

Proor: As stated before, a truss is in stable equilibrium only if the matrix of the 


coefficients a;; is positive definite, or 
Q=3 > d>a,;6,0; > 0 (24) 


where the equality sign applies only to the trivial case 6; = 0 (¢ = 1, 2 --- n). By re- 
arranging the terms in the double sum, this can be expressed [8] in the form 


Q=}3 Zz. S,L,[(e, coth e, — 1)7'(@, + 04)” + (e.)7'(0, — 6{) coth «|, (25a) 
where e¢, is defined by 
. oY Bs : 5 
€E. = . —T (25D 
41, ' ) 
and where 6, and 6; are, respectively, the rotations of the two joints connected by the 
k-th bar of moment of inertia J; . 

Let us now assume that a truss 7 is in stable equilibrium for a given set of internal 
and external forces, and consider the equilibrium, under the same force system, of a 
truss 7’ which is identical with J as to its geometric and elastic properties, except 
that the moment of inertia of its m-th bar has been increased, i.e. 


la fas 


If this bar be a tension member, it follows from Eq. (25b) that 





6This lemma, which has immediate physical appeal, is analogous to a similar principle of Rayleigh 


dealing with the natural frequencies of vibration of elastic systems. 
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and, by inspection of Eq. (25a), that the value of the m-th term in the series has not 
decreased. The same can readily be demonstrated for compression members, which 
correspond to imaginary values of ¢, and for unstressed members. Consequently, for 
any given set of joint rotations, 
, . 
Q’ > Q; 
in view of Eq. (24), the truss 7” is therefore also in stable equilibrium. 

As a next step, it is easy to establish the following 

Lemma II: The value of the ultimate load of a truss is not lowered by increasing 
the moment of inertia of a bar. 

Proor: This is a direct consequence of the previous lemma and of Theorem I. In 
fact, if the prestressing parameters \* (r = 1, 2 --- m) be so selected as to make the 
critical load \, equal to the ultimate load > for truss 7, it follows from Lemma I that, 
for the same \* , the critical load \j of truss 7” satisfies the relationship 


Ai => A =A- 


On the other hand, it follows from Theorem I that the ultimate load d{" of truss T’ is 
the largest of all critical loads of T’, or 


Ao = Ao 
which, in view of the foregoing, implies 
Mo 2 r- (26) 


An upper bound to the ultimate load can now be established by means of the fact 
that a truss of m-th degree of redundancy can be converted into a rigid-link mechanism 
of one degree of freedom by considering it pin-jointed and by removing / members, 
where 1 < m + 1. In general, it will be possible’ to select these 1 members in such a 
way that the mechanism so created is capable of joint velocities corresponding to a 
shortening of all the eliminated bars. 

If now these / members be assigned compressive forces 


S, = —40E,1,(L,)” (k = 1,2, --- D (27) 


and if the remaining bar forces be so chosen as to satisfy the force equations of equili- 
brium (2a) for \, = i > 0, then, by 


THEOREM II: 
No S Xo’. (28) 


In other words, the load parameter found by the method described above represents 
an upper bound to the ultimate load. 

The proof of this follows from the second lemma. In fact, since Eq. (27) defines 
the buckling load of a column which is fixed at both ends, it is apparent that dj’ can 
be interpreted as the ultimate load of a truss which is identical with the truss under 


7Otherwise, the truss does not exhibit an ultimate load, i.e., the load increases indefinitely during 


buckling. 
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consideration except that all bars other than the selected | bars described above have 
infinite moment of inertia. Thus the second lemma implies the relationship (28). 

It should be pointed out that the practical usefulness of the upper bound defined 
by Theorem II is sharply limited. A more detailed discussion of this question can be 
found elsewhere [3]; suffice it to state here that the theorem is much less broad in scope 
than a similar one establishing upper bounds to the collapse load of a plastic-rigid 
frame. [4] This is due to the fact that it has not been possible to construct an infinity 
of readily calculable ‘kinematically admissible’ collapse modes by any simple process 
analogous to the insertion of a sufficient number of yield hinges. Fortunately, except for 
the purpose of estimating the degree of accuracy attained, the interest of the engineer 
is focused on the lower bound only, which can be made to approach the exact value of 
the collapse load as closely as desired. 

Conclusion. Two principles have been derived which establish lower and upper 
bounds to the ultimate load sustained by a buckled redundant truss. Of these, the lower 
bound permits an approach to the exact value from below within any desired degree of 
accuracy. 
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ON SAINT-VENANT’S PRINCIPLE* 


BY 
E. STERNBERG 


Illinois Institute of Technology 


Introduction. The principle bearing his name was introduced by Saint-Venant [1]' 
in connection with, and with limitation to, the problem of extension, torsion, and flexure 
of prismatic and cylindrical bodies. The first universal statement of the principle is 
apparently due to Boussinesq [2], and reads:’ “An equilibrated system of external 
forces applied to an elastic body, all of the points of application lying within a given 
sphere, produces deformations of negligible magnitude at distances from the sphere 
which are sufficiently large compared to its radius.”’ Love [3] writes:* “According to 
this principle, the strains that are produced in a body by the application, to a small 
part of its surface, of a system of forces statically equivalent to zero force and zero 
couple, are of negligible magnitude at distances which are large compared with the linear 
dimensions of the part.” 

As pointed out by v. Mises [4], the foregoing statements are in need of clarification* 
since the forces applied to a body at rest must be in equilibrium in any event. Only 
when the body extends to infinity, and provided we require the tractions at infinity 
to vanish suitably, is it meaningful to speak of the strains “produced” by a non-equili- 
brated system of forces applied to a bounded part of its surface. Moreover, in this instance, 
the strains produced by a given loading are arbitrarily small at points sufficiently far 
removed from the region of load application, regardless of whether or not the loading 
is self-equilibrated.° On the other hand, the stresses and strains at a fixed point of an 
elastic body, in the absence of body forces, may be made arbitrarily large or small by 
choosing the magnitude of the loads sufficiently large or small. These observations further 
confirm the need for clarification. 

What is meant by the statements cited may roughly be expressed as follows:° if 
the forces acting on an elastic body are confined to several distinct portions of its surface, 
each lying within a sphere of radius e, then the stresses and strains at a fixed interior 
point of the body are of a smaller order of magnitude in « as e — 0 when the forces on 
each of the portions are in equilibrium than when they are not. In this comparison 
we must evidently assume that the forces remain bounded as « — 0. The analogous 
interpretation for distributed surface tractions is immediate. 

It should be noted that such an interpretation is implied in the usual applications of 
Saint-Venant’s principle. Moreover, that this is what Boussinesq had in mind is apparent 
from his efforts to justify the principle. With this objective, Boussinesq [2] considered 

*Received January 5, 1953. 

1Numbers in brackets refer to the bibliography at the end of this paper. 

*See [2], p. 298. (‘Des forces extérieures, qui se font équilibre sur un solide élastique et dont les points 
d’application se trouvent tous a l’intérieur d’une sphére donnée, ne produisent pas de déformations 
sensibles a des distances de cette sphére qui sont d’une certaine grandeur par rapport 4 son rayon.’’) 


‘See [3], p. 132. 
‘See also, for example, Biezeno and Grammel [5], where the traditional statement of the principle is 


discussed in detail. 
‘See the general solution to the problem of a semi-infinite medium bounded by a plane, [3], art. 166. 


‘This interpretation follows v. Mises [4]. 
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a semi-infinite body under concentrated loads acting perpendicular to its plane boundary. 
He showed that if the points of application of the loads lie within a sphere of radius e¢, 
the stresses at a fixed interior point of the body are of the order of magnitude e provided 
the resultant force is zero, and of the order ¢’ in case the resultant moment also vanishes. 
Various energy arguments have since been advanced in support of the principle.” 

In 1945 v. Mises [4], in his illuminating paper on this subject, showed with the aid 
of two specific examples that the usual statements of the principle, when properly 
clarified, cannot be valid without qualifications. The two examples chosen by v. Mises 
are the three-dimensional problem of the half-space and the plane problem of the circular 
disk, each under concentrated surface loads.* On the basis of these examples v. Mises 
proposed an amended principle. 

It is the purpose of this paper to supply a general proof of Saint-Venant’s principle 
as modified by v. Mises. The argument is carried on for the case of piecewise continuous 
tractions and is later extended to concentrated forces; it applies to finite and infinite 
domains of arbitrary connectivity.” 

The dilatation formula of Betti. As a preliminary to the proof, we recall here a 
formula due to Betti,"° which is a consequence of Betti’s reciprocal theorem. Let D 
be a regular’ (not necessarily simply connected) region occupied by an elastic medium, 
and let B be the boundary of D (Figure 1). Furthermore,” let u, e;; , and 7;; be a dis- 
placement field, a field of strain, and a field of stress which within D satisfy the funda- 
mental field equations of the linear theory of elasticity in the absence of body forces. 
If r;; gives rise to piecewise continuous surface tractions T = [X, , X, , X;] on B, then 
the dilatation A° = e§ at a fixed interior point Q(é, , & , &) of D is given by,” 


. 


8 oni 

cA? = [ ~~; oo SS. (1) 
“B 1 ~ 2v 

where » and » are the shear modulus and Poisson’s ratio, respectively. In (1), g is a 

displacement field which is defined as follows: 


g=e’+e”’, g’=-gradRkR"', R=|R|, (2) 


where R is the position vector with respect to Q of a point P(x, , x2 , 23) of D (Figure 1). 
Moreover, g” is that displacement field which satisfies the equilibrium equations within 
D and gives rise to surface tractions on B which are equal and opposite to those associated 
with g’. Thus g is characterized by the requirements that (a) it satisfy the equilibrium 
conditions inside D with the exception of the point Q where it must have the singularity 
appropriate to a center of dilatation, and (b) its associated surface tractions vanish on B. 
We emphasize, for future reference, that g is an analytic function of position on any 
analytic portion of B. Betti’s formula remains valid if D is not bounded, provided that 
R’r;; ~0as Ro, 


7See References [6] to [11] and [18]. 

8More recently, Erim [12], applied v. Mises’ analysis to the half-plane under concentrated loads. 

°The restriction that the region be simply connected turns out to be unessential. 

10See [3], p. 235. 

The term “regular region’’ is used in the sense of Kellogg [13], pp. 113, 217. 

Throughout this paper letters in boldface designate vectors; the subscripts 7, 7 assume the values 
1, 2, 3, and the usual summation convention is employed. 

The symbols - and X designate scalar and vector multiplication of two vectors, respectively. 
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The vector field g plays a role analogous to that of Green’s function in potential 
theory. Equation (1) reduces the determination of the dilatation in the second boundary- 
value problem (surface tractions prescribed) to the determination of g’’. We note that 
the function g is completely characterized by the shape of D and the location of Q, 
and is independent of the surface tractions T. This observation provides the key to the 
subsequent proof of the modified Saint-Venant principle in which we are confronted 
with the task of comparing the effects of certain changes in the loading upon the de- 


- Boundary B 





Region D 











Fic, 1. 


formations and stresses at an interior point of the body. Before turning to our main 
objective, we consider two specific applications of Betti’s formula. 

The sphere and the half-space as examples. Let D be a sphere with radius a and let 
Q be its center. Here, trivially, 








» . Al - 2m 
and by (1), (2), (3), 
Q _ 31 _ 2v) ¥ 
ao = 5 9 R dB. (4) 


If the loading, in particular, consists of two equal and opposite concentrated forces, 
each of magnitude L, applied at the endpoints of a diameter and directed toward Q, 
a trivial limit process applied to (4) at once yields, 

3(1 — 2W)L (5) 


S = _—_ -— . 
oo 4n(1 + v)ya” 








396 E. STERNBERG [Vol. XI, No. 4 


This formula was derived by Synge [14], by Weber [15], as well as by F. Rosenthal 
and the present author [16], in each case by entirely different means. 

Next, let D be the half-space x; > 0 and B the plane x, = 0. In this special instance 
a closed representation of g is available’* corresponding to any interior point Q(é, , £2 , &3). 
Indeed, let 


e= (a, oe t,)” + (X2 cates f)* + (23 + §3)°] - (6) 


so that ¢ is the reciprocal of the distance between P(x, , x. , x;) and the mirror image 
of Q in the plane xz, = 0. Then, as is readily verified, 


ae as i. a. 5 ee. - 
le eta. E Ox,’ 26. ans grad 3. @) 
and, according to (2), (6), (7), on the plane z; = 0, 

g = —4(1 — v)[grad R''],.,<0) . (8) 


Equations (1), (8) imply, 


A? = — (1 — 2) | T - grad R'' dB. (9) 
2m /B 
An elementary limit process applied to (9) yields the formula cited by v. Mises'* for 
the dilatation at Q induced by a concentrated load acting at the origin z; = 0. 

Proof of the modified Saint-Venant principle. Let D, with the boundary B, again 
be a regular region’® of space of arbitrary connectivity. Let S,(k = 1, 2, --- m) bem 
non-intersecting closed subregions of B which lie within neighborhoods of m distinct 
points P;* (position vectors r}“’) of B, each S, , Po” being wholly contained within a 
sphere of radius «, (Figure 1). We note that S, need not be simply connected or even 
connected. 

Let u, e;; , 7;; be a solution in D of the field equations of elasticity theory, the body 
forces being absent, which corresponds to piecewise continuous surface tractions 
T = [X, , X. , X;] on B.”” Moreover, let T vanish on B with the exception of the sub- 
regions S, . 

According to Betti’s formula (1), 


“8S 


a? = YAS], cA(SJ = | g-TaB (10) 
k=1 k 


where A°%[S,] represents the “contribution” from the tractions on S, to the dilatation 
A® at a fixed interior point Q of D. We observe that a A°[.S,] possesses individual phy- 
sical significance only if the tractions on S, are self-equilibrated, unless D extends to 
infinity. Furthermore, if D is bounded, either the tractions on each S, are self-equili- 
brated, or there are at least two S, for which this is not true. 

We now examine A°[S,] and, for the sake of convenience, henceforth write S, Po , fo 


“This result is apparently due to Cerruti. See [3], p. 239. 

See Equation (4) of [4]. 

16See Footnote 11. Again D need not be bounded. 

"If D is not bounded, we require r*7;; - 0 asr ©, wherer = |r|, and ris the position vector of a 
point P of D. 
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in place of S, , Pj”, r)"’. Let S and P, be contained in an open simply connected sub- 
region S* of B (Figure 1) which admits a parametrization of the form,’* 


r=r(a,6), ff X%#O0, (6) in =* (11) 


where r = [z, , £2 , Z| here is the position vector of a point P of S*, =* is an open simply 
connected region of the (a,8)-plane, and r(a, 8) is assumed to be at least twice contin- 
uously differentiable in =*. Thus B is assumed to have finite and continuous curvatures 
in S*. The mapping (11) defines a regular curvilinear coordinate net on S*. It is con- 
venient to require that 

r, = r(0, 0). (12) 


Finally, suppose that the closed subregion = of =* is the antecedent in the (a,8)-plane 
: PI p 


of the subregion S of S*. 
From (10), by virtue of the regularity of g on B, we have, on expanding g(a, 8) in 
a Taylor series (possibly with a remainder term) at (0, 0),"° 


cA°[S] = g’ - [ Tdo +g. - [ Ta do + g3- [ TB do + -:-:: (13) 
/S “8S “8S 
where, 
[ do= | |r. X t | da dB. (14) 


Before drawing any general conclusions, let us apply (13) to the example of the 
half-space. Here D is the region x, > 0 and B the plane x; = 0. We may choose the origin 
ct; = Oat P, , employ B as S*, and adopt the parametrization, 


al | 


a=$7; 8B =Ze. (15) 
By (8), in this instance, 


gla, B) — = 402 g —- Gd, £2 bane B, &3] 
(16) 
R = [(, — a)’ + (& — 6)? + &]'” 


Substitution of (16) into (13) yields, except for differences in notation and a constant 
factor, the expansion derived by v. Mises by other means, provided the integrals in 
(13) are replaced with the corresponding finite sums.” 

We now return to (13) and to our main objective, which is to examine the order of 
magnitude of the dilatation at the fixed interior point Q of D in relation to the size of 
the region S, under various assumptions regarding the tractions on S (e.g., if the trac- 
tions on S are or are not self-equilibrated). To this end we first recall the mathematical 


The subscripts a, 6 denote partial differentiation with respect to the argument indicated. The 
existence of such a regular parametrization is assured in the small, provided S* is sufficiently smooth. 
Note that the embedding regions rs , belonging to the various S; , P\", in general, require different 
parametrizations. 

"The superscript zero attached to any function of (a, 8) refers to its value at (0, 0). 

20See Equation (7) of [4]; this equation gives the mean normal stress rather than the dilatation at Q, 
and applies to concentrated forces. The transition from distributed tractions to concentrated forces will 


be discussed later. See Equation (29) of this paper. 
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meaning of the concept of “order of magnitude”. If” 
| f(z)/x’? | <M (M independent of z) as x—0, then f(x) = O(2’), (17) 


that is, f(x) is said to be of the same order of magnitude as x’. It is clear from (17) 
that the question as to the order of magnitude of A*[.S] in (13), with respect to the 
radius ¢, of the sphere enclosing S and P. 
not a variable. 

The foregoing question becomes meaningful, however, if we ask what happens in 
the limit as the region S is contracted to the fixed point P, of B. To make this idea 
precise, consider a one-parameter family of closed subregions S(e) of S* such that for 
every «in 0 < e < « , S(e) together with P, lies within a sphere of radius e, S(¢,) = S, 
and the maximum diameter d(e) of S(e) is a monotone increasing function of «. Next, 
let u(e), e:;(€), 7:;(€) (0 < € < &) be a one-parameter family of solutions in D of the 
field equations of elasticity theory, the body forces being absent, which satisfies the 
following conditions: u(e,) = u; 7,;(€) gives rise to piecewise continuous surface tractions 
T(e) = [X,(e), X.(€), X3(6)] which vanish on B with the exception” of the subregions 
S(e); T(e) remains bounded as e — 0. 

Writing A°(e) for A®[S(e)], we have from (13), (14), 


» , has no meaning since ¢ is a number and 


cA%X(e) = g° . | T(e) do + g°. . [ T(ea do 


+er- | TObdot+---, (18) 
~ S(e) 
[ do=[ |r. Xt | dads, (19) 
* S(e) ¥ Te) 


where =(e) is the antecedent of S(e) in the (a,8)-plane. Furthermore, in view of the 
regularity of the mapping (11), and from (12), 


d(e) = max (a? + p°)'? =0 as e— 0 and 56 = Oe). (20) 
Let ; 
F) =| Ted, MO=[ rxTOde OK<eSa). (21) 
* S(e) #” Ste) 


Thus, F(e) and M(e) are the resultant force and the resultant moment in O(r = 0) of 
the tractions on S(e). Equation (18), in conjunction with (19), (20), (21), now yields 
the following conclusions: 

(a) F(e) = 0 (0 < € < &), in general, implies A°(e) = 0(e’). Exceptions are possible, 
however. In the example of the sphere, we see from (4) that A° = 0 whenever the surface 
tractions are purely tangential, provided Q is the center of the sphere. 

(b) A%(e) = O(e) or smaller if F(e) = 0 io < €S &). 

(c) A%(e) = Oe) or smaller if 


F)=0, | T@ade=0, | TeBdc=0 OK<e<e). (22) 
~ S(e) 


~ S(e) 





21This limited form of the definition is sufficient for our purposes. 
2Recall that S(¢) stands for S;(e) (k = 1, 2, «++ m). The piecewise continuity of T does not refer to 
its dependence on ¢ for which no continuity requirements are imposed. 
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Within terms of 0(e*), Equations (22) are equivalent” to the 12 scalar conditions, 


0 (0 <<s < €). (23) 


II 


[ XQ do =0, / X(oz, de 
S(e) 


~ S(e) 
Hence, if the tractions on S(e) satisfy (23) then A°(e) 
(d) Equations (23) imply, 


Fe) =0, Me =0 O<eX< @), (24) 


0(e‘) or smaller. 


but the converse is not true, in general. Within terms of 0(e*) the equilibrium conditions 
(24), in view of (21), (11), are equivalent” to 


F(e) = 0, x [ 


Se 


T(a do + 12 X / Ti)pdeo=0 O<e<e). (25) 
) S(e) 


Again, these conditions are met if (22) hold, but (22) do not follow from (25). Thus, 
if the tractions on S(e) are self-equilibrated A°(e) = O(e*) or smaller. This conclusion 
contradicts the interpretation of the traditional statement of Saint-Venant’s principle 
cited in the Introduction, according to which the order of magnitude of A°(e) should 
always be smaller when the tractions on S(e) are self-equilibrated than when they are 


not. 
(e) Suppose, in particular, the tractions on S(e) are parallel, so that 


T(a, 8; ©) = kTfq, 8; © in r(6) (0 <€< &), (26) 


in which k is a fixed vector and T' a scalar function. Here (23) are satisfied if and only if 
(24) hold for every choice of k, that is, if and only if the system of parallel tractions 
remains in equilibrium under an arbitrary change of its direction, the magnitude and 
sense of the tractions being maintained (“astatic equilibrium”). Therefore, in the event 
the tractions on S(e) are parallel and in astatic equilibrium, A°(€) = O(e*) or smaller. 

(f) Suppose (26) holds and, in addition, 


k- (r. xX rs) + 0, (27) 


so that k is not parallel to the tangent plane of S* at P, . Then (25) imply (22). Hence, 
if the tractions on S(e) are parallel to each other, self-equilibrated, and not parallel to a 
tangent plane of S*, then A°(e) = O(e*) or smaller. The analogous conditions for concen- 
trated forces were satisfied in the special example investigated by Boussinesq [2] and 
described in the Introduction. 

We have so far considered only the order of magnitude of the dilatation. According 
to Lauricella,™ the strains e,2 at an interior point Q of D, in the absence of body forces, 


admit the representation, 
= | T+ gab, (28) 
B 


which is analogous to that given by Betti’s formula (1) for the dilatation A°. Here the 
g,; are displacement fields which satisfy the equilibrium conditions inside D, with the 
*%In case S* is plane, and for the parametrization (15), the equivalence is exact. 
*4See [3], p. 216. 
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exception of the point Q where they have certain prescribed singularities,”” and which 
give rise to vanishing surface tractions on B. Clearly, the g;; are again regular on any 
regular portion of B, and the previous argument applies without modification.. Therefore, 
the conclusions listed under (a) to (f) remain valid if A°(e) is replaced with e,9(e), and 
hence with 7,5(e). If a rigid displacement of the whole body is excluded, the conclu- 
sions also apply to the displacement u°(e), as follows from Somigliana’s representation 
of the displacement field.*° Specific illustrations of these general conclusions were given 
by v. Mises in [4].?’ 

Extension to the case of concentrated forces. Remarks. ‘The theorem proved in the 
preceding section is readily extended to concentrated forces. With reference to (13), 


let the concentrated forces T, = [X;", X;", X;"] (n = 1, 2, --- N) be applied at the 
points A, of S*, the points of application together with P, lying within a sphere of 
radius ¢, . Let r™’ = [z,"’, x,"’, x;"] = r(a, , B,) be the position vector of A, . Now 


° . ‘ ‘ » - . e ° ¥( - 
consider an S consisting of N non-intersecting closed subregions S™ of S* such that 
each S*”’ is simply connected and contains A,, in its interior. Proceeding to the limit in 
(13) as S™ is contracted to A, while fs... T ds — T, , we obtain,” 


\ N 
cA°[S] = g°- DT, +82 - Do Tron 


+g3->T.B,+-::. (29) 


The conclusions previously reached, therefore, remain valid for concentrated forces 
provided the integrals in (18) to (25) are replaced with the corresponding finite sums, 
and provided O(e’), O(e*), O(c’) are replaced with O(1), O(e), O(c’), respectively. 

At this place we discuss an example which may serve to clarify the implications of 
the theorem established earlier. Consider a bar of the general shape indicated in Figure 2. 
Let the bar be acted on by the two equal, opposite, and collinear concentrated loads, 
each of magnitude L, the region S consisting of the two points of application A, and A, . 
According to conclusion® (d), in this instance 7,4(€) = O(e) or smaller, where Q is, 
say, the fixed interior point shown in Figure 2. Since we may choose the two ends of 
the bar as close together as we wish, a careless interpretation of (d) may lead to the 
absurd prediction that the stresses at Q, for arbitrarily large fixed magnitudes of the 
loads, can be kept as small as we wish. The paradox is resolved by observing that the 
shape of the bar is given once and for all, and that a definite gap, however small, exists 
between the two ends of the bar. The statement 7;“(e) = O(c), in view of the definition 
(17), merely implies | 74(e) | < Meas ¢ 0, where M is a positive number independent 
of e«. Thus, by contracting the load region sufficiently, say toward A, , while maintaining 
the loading within the «sphere in equilibrium, | 7,4(e) | can be made arbitrarily small. In 
this process of contraction, however, the end of the bar which carries the point A, , 
will eventually cease to lie within the contracting «-sphere so that the entire character 

%g,.,, for example, at Q has the singularity appropriate to a force-doublet the axis of which is parallel 
to the z)-axis. 


See [3], p. 245. 
These examples refer to concentrated forces (see the next section of the present paper). 


*%Note that S here consists of the points A,(n = 1, 2, --- N). 
Modified for the case of concentrated forces. 
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of the loading changes, that end being now free from loading. This example provides 
additional evidence for the vagueness of the traditional statements of Saint-Venant’s 
principle, quoted in the Introduction. 

Hoff [17] pointed out important limitations inherent in certain conventional engi- 
neering approximations which are usually based on an appeal to Saint-Venant’s principle. 
These observations are consistent with the results established here. For Saint-Venant’s 

















Fig. 2. 


principle is a statement about relative orders of magnitude and does not tell us anything 
about the extent of the region within which a self-equilibrated system of tractions, 
applied to a portion of the surface of an elastic body, “materially” influences the stress 
distribution in the body. 

In conclusion, it may be well to repeat v. Mises’ remark [4] to the effect that the 
present theorem does not preclude the validity of a stronger Saint-Venant principle 
for special classes of bodies, such as “thin’’ plates or shells and “long”’ cylinders. 


BIBLIOGRAPHY 


1. B. de Saint-Venant, ‘““Mémoire sur la torsion des prismes,”” Mém. Savants étrangers, Paris, 1855. 

2. M. J. Boussinesq, ‘‘Application des potentiels,’”’ Gauthier-Villars, Paris, 1885. 

3. A. E. H. Love, “A treatise on the mathematical theory of elasticity,’ 4th ed., Dover Publications, 
New York, 1944 

1. R. v. Mises, “On Saint-Venant’s principle,’ Bull. Amer. Math. Soc., 51, 555 (1945). 

5. C. B. Biezeno and R. Grammel, “Technische Dynamik,” J. W. Edwards, Ann Arbor, Michigan, 
1944. 

6. R. Southwell, “On Castigliano’s theorem of least work and the principle of Saint-Venant,’’ Phil. 
Mag. (6) 45, 193 (1923). 

7. J. N. Goodier, “‘A general proof of Saint-Venant’s principle,” Phil. Mag. (7) 23, 607 (1937). 

8. J. N. Goodier, “Supplementary note on ‘A general proof of Saint-Venant’s principle’,”’ Phil. Mag. 

7) 24, 325 (1937 

9. O. Zanaboni, “Dimostrazione generale del principio del De Saint-Venant” Atti Acc. Naz. Lincei 
25, 117 (1937). 

10. O. Zanaboni, “Valutazione dell’ errore massimo cui da luogo l’applicazione del principio del De 
Saint-Venant,”’ Atti Acc. Naz. Lincei 25, 595 (1937). 

11. J. N. Goodier, ‘“‘An extension of Saint-Venant’s principle, with applications,” J. Appl. Phys. 13, 167, 
1942. 








402 


12. 


13. 
14, 


. E. Sternberg and F. 


E. STERNBERG [Vol. XI, No. 4 


K. Erim, “Sur le principe de Saint-Venant,” Proc., Seventh Int. Cong. for App]. Mech., London, 

1948. 
O. D. Kellogg, ‘“‘Foundations of potential theory,”” The Murray Printing Company, New York, 1929. 
J. L. Synge, “Upper and lower bounds for the solution of problems in elasticity,’ Proc., Royal Irish 


Ac. (A) 53, 41 (1950). 


. C. Weber, ‘Kugel mit normalgerichteten Einzelkriften,” Z. Angew. Math. Mech. 32, 186 (1952). 


osenthal, ‘‘The elastic sphere under concentrated loads,” J. Appl. Mech. 19, 


413 (1952). 
N. J. Hoff, ‘“‘The applicability of Saint-Venant’s principle to airplane structures,” J. Aero. Sci. 12, 


455 (1945). 


. P. Locatelli, ““Extensione del principio di St. Venant a corpi non perfettamente elastici,’”’ Atti 


R. Ace. delle Sci. Torino 75, 502 (1940). 














403 


ON THE EXPANSION OF FUNCTIONS IN TERMS OF THEIR MOMENTS* 


BY 
H. 8. GREEN AND H. MESSEL 


University of Sydney 


Summary. A general method is devised for the reconstruction of functions of a 
continuous variable from their moments. An analogue is given for functions of a discrete 
variable. An application is given to the solution of partial differential equations with a 
given initial condition. 

1. Introduction. Recently, in some work connected with the lateral spread of showers 
of cosmic-ray particles in their passage through the atmosphere, the authors were faced 
with the problem of finding the radial distribution function of the particles with respect 
to the shower axis. This function was known to be the solution of a partial differential 
equation, from which, however, only the moments of the function could be conveniently 
determined. It was then necessary to devise a method for reconstructing the function 
from its moments. From the physical nature of the problem considered by the authors, 
it was known that such a function existed; however, situations could arise in which 
even this knowledge was not available. Criteria for the existence and uniqueness of 
distribution functions corresponding to a given set of moments having been the subject 
of extensive study [1], and analytical procedures for the determination of the distribution 
functions have been described. However, it has also been pointed out [2] that these 
methods have little value in practice. The only method known to the authors which is 
well adapted to application has been given by Spencer and Fano [3], but this was not of 
sufficient generality for our purpose. We therefore developed a method of considerably 
wider applicability. 

The expansion of a function in terms of its moments is a problem which arises not 
only in methematical physics, but also in many branches of statistics. In the hope that 
the method may prove useful to workers in fields other than our own, we now give a 
brief account of the theory in its most general form. 

2. Expansions in terms of the 6-function. The use of the 6-function has been well 
established in quantum mechanics [4] and pulse theory; however, the rigour of the 
mathematical procedures in which it is used has sometimes been questioned [5], and we 
therefore state at the outset the unambiguous meaning of an equation of the type 


ro) 


f*(x) = w*(z) DY a(x) 8 (2), (1) 
k=0 
where the superfix represents the number of differentiations of the 6-function with 
respect to the arbitrary variable x. The first k derivatives of the function a,(x) must 
exist for x = 0, but w(x) may have any kind of singularity there. Then (1) will be held 
equivalent to the assertion 


r’ dx g(x) f*(x) _ y (—1*4?) 


oe w*(x) 


q(x) = g(x)a,(2) (2) 


*Received January 28, 1953. 
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where q(x) is an arbitrary function regular at = 0, and a and b are any positive numbers 
less than the radius of convergence p of the power series 


g(x) = >> g(O)a*/k!. (3) 
Let 
f(x) = f*(x)/w*(az); (4) 
then since 
| q(x) f(z) dx = >> qg“(0) | a‘ f(x) dx/k!, (5) 
one has, by comparison with (1) and (2), 
f(z) = > (-1)'fu 6 (x)/kI, (6) 
fu = | a‘ f(x) dx. (7) 


Thus, any function, integrable in the range —a < x < b can be expanded in this range 
as a series of derivatives of the 6-function, with coefficients which are proportional to 
the moments (7) of the function f(z). This result is easily extended to functions of any 
number of variables. 

3. Expansions in orthogonal polynomials. The result (6) requires some elaboration 
in order to provide a practical method for the determination of f(x) when the coefficients 
fy of the series are known. The method which we shall adopt is to make a formal 
expansion of 6“’(x) in series of orthogonal polynomials;* but first we obtain the ex- 
pansion of the function f(z) which anticipates the result. Let 


) 


f(z) = w(2) 2 frSn(X); (8) 
where S,(z) is a set of polynomials which will be defined presently; and if f(z) has 
singularities at z = 2, , 22, +--+ , 2% , Where | z, | < a, | 2, | < b, let w(x) be a ‘‘weight- 
function” with a singularities of the same type, so that f(x) /w(x) is regular for | x | < a, 
|x| < b. The polynomials S,(x) are defined by 


W co W ove Win | 
} 
Wes w Win) | 
S,(z) = tee se ves es for n> i, (9) 
u Ww u ' 
l 7 see - 4 


*For the general theory of such polynomials, c.f. Szegé [6]. The cases of Legendre and Laguerre 


polynomials are known [1]. 
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where S,(x) = 1, and 
ob 


Wa = | w(x)a* dx. (10) 


Then these polynlmials satisfy the orthogonality relations 


ab 





| w(x) S,,(z)S,(x) dx = Nidan , (11) 
where the normality constants are given by 
N, = A,A,-: ; 
W io Wo1) ee Win) 
Wi, W 2) soe Wn+1) 
A, = for n> 0, (12) 
Wir Wine) ee Won) 
A., =1. 
There is obviously an expansion of zx‘ of the form 
z= dD &.1Si(2); (13) 
i=0 
hence, for values of x less than the radius of convergence p of the power series 
f(x) /w(2) = > v2" (14) 
k=0 


which exceeds both a and b, there will exist a convergent expansion of the type (8), 
with coefficients given by 


f, = >> vike.n (15) 


k=n 
The coefficients f, are most readily obtained by multiplying (8) by S,,(x), and inte- 


grating from —a to b; thus: 


+b 


Nf, = | f(x) S,(a) dx 





Wray W 4) tee Wn) 
Wi W 2) [+e W (n+1) 

Win-i Wn) es Wi2n-1) 

i 0 F 1 rae Fin 
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The function f(z) has then been expanded in a series of orthogonal polynomials, the 
coefficients of which are linear combinations of its moments f,,) . It is necessary that 
the quantities A, defined by (12) should be positive, if a monotonic function f(x) exists. 
The conditions for this have been discussed in detail by Shohat and Tomarkin [1], 
Chap. 1. 
The result (8) may also be deduced from (1) by expanding 6“’(x) in terms of the 
orthogonal polynomials S,(x), thus: 
co 
6°°’(a2) = w(x)(—1)"k! z. Aicey Si(x)/N, (17) 
l=k 
where A,,;) is the minor of w,,,,) in the last row or column of the determinant A, . Sub- 
stituting (17) into (1), and interchanging the order of the summations, one then obtains 
(8), which is therefore one way of interpreting the formula (1). 
This procedure is easily adapted to distribution functions F'(n) of a discrete variable n. 


Here we consider expansions of the type 


F(n) = Win) > F,S,An) (18) 
where 
W vo) W ay eee W ir) 
S,(n) = for > 1, (19) 
W r=—1) W or) lig W iee—1) 
1 n eee n” 
S,(n) = I] 
and 
Wir) = >, Wn)’. (20) 


n=O 


It is necessary for the existence of the W,,) that W(n) should decrease at least expo- 


nentially for large values of n. 
The determinants S,(n) satisfy the orthogonality relations 


>> Wn)S,(n) Sn) = Nedq,r 5 (21) 
where the N, are defined by equations precisely analogous to (12); 
The coefficients F, in (18) are determined by the relation 


> S,(n)F(n) = NF, (22) 
n=0 
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assuming the convergence of the series. Substituting for S,(n) from (19), this gives 





W wo) 7 W ve) 

NF. = | for q>1, (23) 
| W wq-1) er W 2e-1) 
| Fo Sin Fw) 

Nolo = Feo) ; 

where 
Fi = >> F@)n’* (24) 
n=0 


are the moments of the distribution, assumed to be finite. 

In order to ensure that the series (18) should converge, it is necessary to impose 
rather stringent conditions on W(n), depending on the nature of the function F(n). 
In practice, if the weight-function W(n) is chosen to approximate fairly closely to 
F(n), the series is either convergent, or else is asymptotic and equally suitable for com- 
putational purposes. 

If one has no previous knowledge of the nature of the function f(x) or F(n), but is 
given only the moments f,,) or F:,) respectively, defined for the range 0 = a < x (orn) < 
b < o, the weight-function should be chosen in the following manner. The asymptotic 
behaviour of f,,) for large k is compared with that of the moments of a weight-function 
of the form 

w(x) = exp (—Bz*) (B > 0, A > 0), (25) 

namely 
Wey = (AB“*?/4)“"P6(k + 1)/A} ~ (k/BeA)”4 for 
Way ~ ve */(k +1), dfinite, a=0 


A numerical comparison of these expressions with the known values enables one to 
choose a proper value of A, and thus to determine the nature of the singularity of f(x), 
if any, at the origin. Any other singularity will be revealed by a failure of the series (8) 
to converge for large values of x; and in difficult cases it may be necessary to change 
the origin. A similar procedure may be followed to identify the “singularities” for a 
discrete variable. 

The value of B in (25) should be determined as follows. Assuming that numerical 
values are available for the first m + 1 moments, the series (8) is terminated at the 
mth term, and B chosen so as to give correctly the value of the (m + 1)th moment, 


b= oa a=0 (26) 











thus: 
W (0) Wn) 
m 
S ime) _ pe (27) 
=0 
" Wn-1) W (2n-1) 
W (m+1) nit oa W(m+n+1) 
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This procedure for the determination of B requires only the solution of an algebraic 
equation at the (m + 1)th degree. If the weight-function is well chosen, it is our ex- 
perience that a good approximation to f(x) will result from the use of the first few 
moments only. 

4. Solution of differential equations. We shall demonstrate the utility of the method 
described above in the solution of partial differential equations with a special type of 
boundary condition. 

Consider the general equation 

O.f(t, x) = O.f(t, 2), (28) 


where O, and O, are operators of the form 


l 9 
== = q(t, x) 
p(t, x) at 
. ah 
0. = >. a,(z) —, (29) 
k=0 Ox 


with the initial condition 
f(0, 2) = 6(2) (30) 


By an iteration procedure a solution can be found of the form 


f(t, 2) = > Y¥,At, 2) 31) 


where 
i) ‘4] 
— {q(t, x)V,+:(¢,2 )} = p(t, x) ym ate) —— VAt 2). 32) 
ot -* — ar 
One may choose ¥,(¢, x) = 0; then g(t, x) ¥, (t, z) must be independent of ¢ and, according 
; 4 q | 
to (30), 
V,(é, x) = g(0, x) 5(x)/q(t, x) (33) 
Furthermore, 
= ig - 0° 2 2 
W(t, x) = iq, Xs | p(t, x) a a,(x) {q(O, x) 6(x)/q(t, x)} dt, (34) 
/0 Ox 


ete. 
The solution of (28) obtained in this manner is of the form (1), and can be reduced 


to the form (6), by use of the general formula 


nee — k'—1) ! (kel) a 
a(x) 6 (2) = a Mk = pi Go ‘tz. 35) 


Thus, the moments of the function f(t, x) with respect to the variable x are obtained 
immediately as functions of ¢; and, using the procedure described in the previous section 
for the reconstruction of a function from its moments, an explicit solution is obtained 
of the differential equation which satisfied the boundary condition (30). 
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The solution f’(é, x),f,28) for the initial condition 
f'(0, x) = A(z) (36) 


is then easily obtained in the form 
it 2) = | f(t, 2 — 2)N@) a’. (37) 


An almost identical method of solution may be devised for equations of the type 
0, f(t, Xr, gan La) = er x a Se Lm) 


which, in the application to cosmic ray shower theory, we have actually [7] succeeded 


in solving. 


REFERENCES 
1. J. A. Shohat and J. D. Tamarkin, The problem of moments, American Math. Soc., Mathematical 
Surveys 1 (1943). 
M. G. Kendall, The advanced theory of statistics, Vol. 1, Griffin, London (1943). 
L. V. Spencer and U. Fano, Phys. Rev. 81, 464 (1951). 
P. A. M. Dirac, The principles of quantum mechanics, Oxford, Clarendon Press (1947). 
J. V. Neumann, Math. Grundl. d. Quantenmechanik, Berlin, Springer, 1932. 
G. Szegé, Orthogonal polynomials, American Math. Soc. Colloquium Publications 23. 
H. Messel and H. 8. Green, Phys. Rev. 87, 738 (1952). 


a> orm W OD 


~I 




















411 


LINE LOAD APPLIED ALONG GENERATORS OF THIN-WALLED CIRCULAR 
CYLINDRICAL SHELLS OF FINITE LENGTH* 


BY 
N. J. HOFF, JOSEPH KEMPNER, anp FREDERICK V. POHLE 
Polytechnic Institute of Brooklyn, Brooklyn, N. Y. 


Summary. Donnell’s differential equations of the thin circular cylindrical shell are 
integrated in the case when the loads are radial forces or circumferential moments 
distributed sinusoidally along a generator. Closed form expressions are obtained for 
the displacements, internal moments, and the membrane stresses. In addition, loads 
distributed uniformly along a segment of a generator and concentrated loads are dis- 
cussed and radial forces are combined into a longitudinal moment. 

1. Introduction. One of the most common elements of structures and machinery is 
the thin-walled circular cylindrical shell. When other elements are attached to it, forces 
and moments are likely to be transmitted to the shell across the areas of contact because 
of gravity and inertia effects, and in consequence of thermal expansion. These forces 
and moments can often be represented, with an accuracy sufficient for engineering 
purposes, as loads distributed along a short segment of a generator. For this reason it is 
of practical interest to investigate the deformations and the stresses of thin-walled 
circular cylindrical shells loaded along generators. 

The problem of loads along generators was first solved by Finsterwalder [1] in 1932 
in his investigation of the disturbance stresses arising at the free edge of a circular 
cylindrical shell-type roof structure. In his approximate theory the longitudinal moment 
M, and the torque M,, in the shell were disregarded and a single partial differential 
equation in the circumferential moment M, was derived. This approach was further 
simplified by Schorer [2] in 1935 whose differential equation in M, consisted of only 
two terms. In the same year Dischinger [3] derived a rigorous solution. In order to 
reduce the large amount of work necessary to obtain numerical results from the rigorous 
theory, Aas Jakobsen [4] developed an iteration procedure in 1939 through which the 
equilibrium conditions could be expressed with sufficient accuracy in terms of the radial 
displacement w alone. In 1941 he showed [5] how this approach could be used in the 
calculation of the effects of concentrated loads. 

The concentrated load problem was attacked independently by S. W. Yuan [6] in 
1946 who solved Donnell’s single eighth order differential equation in w for an indefinitely 
long shell. The surface loading was represented by a Fourier series in the circumferential 
direction and by a Fourier integral in the longitudinal direction. The method of images 
was then used to obtain a solution for the cylinder of finite length. In the same year 
Odgqvist [7] gave a closed form solution for the deflections of finite cylinders subjected 
to sinusoidally distributed line loads and infinite series solutions for the deflections of 


*Received Jan. 29, 1953. This work was performed under a consulting contract with the Knolls 
Atomic Power Laboratory, Schenectady, N. Y., operated by the General Electric Company for the 
United States Atomic Energy Commission. The authors are indebted to the company, to the AEC, 
and to Mr. Daniel R. Miller, the project supervisor, for their permission to publish this article. 
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finite cylinders subjected to concentrated loads. These results were. obtained from 
Schorer’s equation. 

Although the work of German and Scandinavian civil engineers seems to indicate 
that the moment M, and the torque M,, are unimportant in the balance of forces and 
moments in a reasonably long shell-type roof subjected to distributed loads, the same 
conclusion need not necessarily hold when line loads or concentrated loads are applied 
to the structure. Indeed, in the experiments performed by Schoessow and Kooistra 
[8], in which forces and moments were transmitted to a shell by means of comparatively 
large-diameter pipes, the ratio of the maximum longitudinal bending stress to the 
maximum circumferential bending stress was found to be close to one-half. Therefore 
some doubt arises regarding the range of applicability of Schorer’s equation. But Odqvist’s 
paper, based on this equation, is the only one to present results from which the stresses 
in the shell can be obtained without extensive calculations. This situation prompted 
the development of the new solution given in this paper. 

In the derivations that follow the stipulation of the conditions of equilibrium in the 
forms presented by Love [9], Fliigge [10], Biezeno and Grammel [11], and Timoshenko 
[12] is replaced by the simplified statement proposed by Donnell [13] and recently 
recommended by Batdorf [14]. In addition to the eighth-order differential equation in 
the single dependent variable w, used and re-derived by Yuan, two fourth order equa- 
tions are satisfied rigorously. The dependent variables are u and w in one equation, 
and v and w in the other; w is the radial displacement, u the longitudinal displacement, 
and v the circumferential displacement. Closed form solutions are given for the dis- 
placement, moment, and membrane stress quantities arising from sinusoidally distributed 
line loads. Series of these solutions are capable of representing the effects of loads constant 
over segments of generators as well as those of concentrated loads. 

2. Basic equations. In the absence of surface and body forces Donnell’s [13] equili- 
brium conditions of an element of a thin-walled circular cylindrical shell can be stated in 


the following form: 


V'w + 4K‘w. 22 = 0, (1) 
Vu = Worse — Wee 5 (2) 
Vv = (2+ vwsee + Wives (3) 
where the subscripts following a comma indicate differentiation, 
4K* = 12(1 — »*)(a/h)?,  D = Eh®/(12(1 — v)], (4) 
and VY’ is Laplace’s operator 
V* = (0°/dz*) + (8°/d¢’). (5) 


The non-dimensional distances and displacements are defined by the equations 


x = (xz*/a), u = (u*/a), v = (v*/a), w = (w*/a), (6) 


and z* is the distance measured in the axial direction along a generator, ag that measured 
around the circumference, and u*, v*, and w* are the displacements in the axial, circum- 
ferential, and radial directions, respectively, as shown in Fig. la. The remaining geo- 
metrical and physical quantities are a, the mean radius of the shell; h, the thickness of 
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the shell; L, the length of the shell; 2, Young’s modulus of elasticity; and », Poisson’s 
ratio. 


Fic. la. 


In Donnell’s approximation and with the present sign convention the strains ¢ and 
the curvatures « of the median surface of the shell are defined as 


é =U., kz = (1/a)w,22 , 
Rip Ug O'R, Ky = (1/a)W, o¢ ’ (7) 
Yzre _ Ui» + U2 ’ Keo - (1/a)w,., P 


The membrane stresses can be given as 
[E/(1 — v*) ee + vey) 
[E/(1 — Poe. + €,) 
Tee = [E/21 + )}re¢ = [E/21 + »))(u,, + 2.2). 


The moment-resultants per unit length are indicated as right-hand vectors in Fig. 
lb. They are 


II 


[E/(1 ia! v*)\(u,s i vw), 
[E/(1 — v*)\ou,. +2,. — w), (8) 


Oz 


I 
I 


Te 


o Qy Myx 
se 
+9 
1+" ox 


whe? 


Fie. 1b. 
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ll 


M, = —D(k. + vx.) —(D/a)(w.22 + WW. oo); 


—(D a)(W, o¢ + VW 2x); (9) 


Il 


M, — Divx, + xy) 


Il 


M,, = —M,, = (1 — »)Dx., = [(1 — ») D/alw.., . 


The shear forces per unit length which accompany these moments and act in the 
radial direction are 


Q. (1/a)(M,.. + M,y..¢) = —(D/a’)(w. sre + W.2¢¢) 


] 


(10) 


Il 


Q, = (1/a)(M,., — Mug,2) = —(D/a?)(w,o¢¢ + W.2x9) 


If the equivalent shear (1/a)M,,,, corresponding to the distributed torque is added 
to Q, in agreement with Kirchhoff’s suggestion (see, for instance, p. 90 of Reference 
12), the total effective shear force per unit length is 


Q..et. = —(D/a’)[w. cee + (2 — Yw.ezy] (11) 


3. Sinusoidally distributed radial force symmetric with respect to tr = 0. It is 
assumed that the edges z* = (L/2) are simply supported. This condition implies that 
the radial displacement w vanishes and so do the axial membrane stress and the axial 
moment resultant, but not the membrane shear stress. If the expressions just given 
for these quantities are considered, the simple-support condition can be reduced to the 
following requirements: 


U2=v=w=vw,,, = 0, when x = +(1/2)(L/a). (12) 


The load is applied along the generator g = 0 in the radial inward direction and is 
2V cos nz with n = mr(a/L) and m an odd integer. With —z < ¢ < 7, the solution is 


assumed as 


Wor = We) + w(2x — g) + w(2e + ¢) + wl4r — ge) + Wl4r +o) +-::, (13a) 


Uror = Uy) + ul(2er — o) + U2r + 9) + Ul4r —G) + Ul4r+et+:::, (13b) 


Vor = Vy) — v(2r — g) + v(2r + ¢) — vV(4r—g) + V4r+y)—:::, (13¢) 
where the dependence on z is not shown but is understood, and 

wy) = Ae”’® cos nz, (14a) 

u(y) = Be’* sin nz, (14b) 

v(y) = Ce’* cos nz, (14¢) 


where A and 7 are assumed real and B, C, and p may be complex. These displacements 
can be imagined to be contributed by two fictitious sheets of width L, each beginning 
at ¢ = 0 and winding around the actual shell indefinitely many times, one clockwise 
and the other counterclockwise. Only those solutions will be considered in this analysis 
in which the real part of p is negative. 

It can be seen that the assumed deflections satisfy the boundary conditions at 
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x = +(1/2)(L/a). In the circumferential direction they are continuous everywhere 


except at ¢ = 0. By a suitable choice of the conditions along the generator g = 0, con- 
tinuity can be enforced there and the prescribed line load can be applied to the shell. 
Worked-out numerical examples have shown that in most cases the first term of the 
infinite series suffices for a satisfactory approximation to the deflections, and two terms 


are needed only infrequently. 
The boundary conditions at g = 0 are: 


(1) For reasons of symmetry 
w, = 0 when g = 0. (15) 
(2) To eliminate gaps caused by the circumferential displacements 
v=0 when g = 0. (16) 


(3) The shear stress r,, acts in the same direction along the edges of the two sheets 
at ¢ = 0. In the absence of a longitudinal line load, 7,, must vanish. Hence 


u,tv, =O0 when eg = 0. 
But v,, vanishes at ¢ = 0 in consequence of (16). Therefore the requirement reduces to 
u,.=0 when g = 0. (17) 


(4) Finally the total effective shear force Qy,,, , given in (11), must be equal to 
one-half the applied radial load. However the second term in the brackets vanishes in 


consequence of (15) with the result that 
(D/a*)w..¢¢ = V cos nz. (18) 
4. Enforcement of the differential equations. Substitution into (1) of the expression 
assumed for w yields 
[(p’ — n°)* + 4K'n‘|w = 0 (19) 
which must be satisfied identically. Hence 
p =n? + V2 (-1)“nK. 
If the fourth root of —1 is taken as (V2, ‘2)(1 —7) then 
pi =n +(1 — ink, (20) 
and if p, is written as 


a ey + 7B, ’ 


where a, and £, are positive real numbers, one has 


ai — Bi = n(n+ K), (20a) 
a8, = (1/2)nK, (20b) 
B, = (1/2)(nK/a,). (20c) 


Finally 
ai = (n/2){(n + K) + [(n + K)? + K’]'”}. (20d) 








416 N. J. HOFF, JOSEPH KEMPNER, AND FREDERICK V. POHLE [Vol. XI, No. 4 


The square root is taken with the positive sign because a, must be real by assumption. 
When the choice is 


(—1)'* = (V/2/2)(-1 — 0), 
the value of p is designated as p, . The above equations are replaced by 


D2 = —a2 + 1B. ’ 


and 
as — 8; = n(n — K), (21a) 
QB. = (1/2)nK, (21b) 
B. = (1/2)(nK /a.), (21c) 
az = (n/2){(n — K) + [m — K)’ + K’*]"’}. (21d) 


The conjugate solutions, corresponding to the remaining two fourth roots of —1, 
namely 
mam? 1B; , a = ae 1B. 
need not be discussed here in detail because they are governed by the same relation- 


ships as p; and 7p, , respectively. 

Of course the solutions for a, and a, of (20d) and (21d) can be negative as well as 
positive. The negative values are rejected because they lead to positive real parts of the 
exponent p and thus yield displacements which increase without bounds as ¢ increases. 

If the deflections assumed in (14) for u and w are substituted in equilibrium con- 


dition (2) one obtains 


Bip’ — nn’)? = An(n’ + p’). (22) 
Consequently 
B = Avn*[1/(p? — n’)?] + Anlp’/(p’ — n’)’). (23) 
When p = 7, one has 
(p, — n*)? = —1(2K’n’), 1/(p; — n®)? = +[t/2K?n’], 


pi/(p; — n°)? = (1/2Kn) + i[(n + K)/2K°n]. 


Equation (23) becomes 


B = (A/2K’){K + i[(1 + »)n + K]} when Pp=DpP. (24) 
When p = p, the equations following (23) become 
(p, — n’)? = 1(2K’n’), 1/(p; — n’)? = —[i/2K’n’], 


pi/(pa — n'y? = —(1/2Kn) — il(n — K)/2K*n] 


and thus 
B = (A/2K’){—K — i[(1 + »)n — K}} when Pp=?P.. (25) 
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Substitutions from (14a) and (14c) into (3) result in 
C(p* — n’*)’ = Aplp* — (2 + »)n’). (26) 
Manipulations similar to the ones shown lead to 
C = (A/2K’n){ —a,K + B,[(1 + »)n — K] + i{8,K + a,[(1 + »)n — K]}}} (27) 
when p = p, and to 
C = (A/2K'n)laxK — Bal(1 + »)n + K] + i{-6.K — ool +o)n + K]}} (28) 


when p = pz. 
Instead of making use of all the four exponents p, to p, , it is more convenient to 


give the radial deflections w in the following form: 
w= [(A,R + ADe”"* + (A3R + A,De”’*] cos nz. (29) 


Here the expressions in parentheses are operators with R signifying the real part and J 
the imaginary part of the complex function upon which the operation is performed. 
Equation (29) can be written explicitly as 


w = [Ae “** cos Big + Ae *** sin Big + Ase “** cos Bap + A,e™*** sin B.¢] cos (nz). 
(30) 


From (14b), (24), and (25) one obtains 
u = (A,R + A,J)(1/2K*){K + i[(1 + v)n + K]}e~***(cos Big + isin B,g) sin (nz) 
+ (A,R + A,J)(1/2K°){—-—K — i[(1 + »)n — K]}e™***(cos Bg + isin Bg) sin (nz). (31) 


When the operations indicated are performed, (31) becomes 


u = (A,/2K*)e***{K cos Bye — [(1 + »)n + K]sin B,¢} sin nz 


(A,/2K’)e*'* {{(1 + »)n + K] cos Big + K sin B,g} sin nz 
(32) 


+ (A;/2K*e***{—K cos Be + [(1 + »)n — K] sin Bg} sin nz 
+ (A,/2K’)e™**°{—[(1 + »)n — K] cos Bg — K sin 8,9} sin nz. 
The circumferential displacement is found in a similar manner: 
vy = (A,/2K’n)e***} {-—a,K + B,[(1 + »)n — K]} cos Bie 
— {8,K + a[(1 + »)n — K]} sin B,g} (cos nz) 
+ (A,/2K*n)e~*"*} {8,K + a[(1 + »)n — K]} cos Bip 
+ {-aK + 6,[(L +n — K]} sin 8,¢}(cos nz) 
+ (A3/2K*n)e"***| {a2K — B[(1 + »)n + K]} cos Bp 
+ {82K + az[(1 + »)n + K]} sin 6,¢}(cos nz) 
+ (A,/2K*n)e~***{ { —B.K — a[(1 + »)n + K]} cos Bry 
+ {aK — B,[(1 + »)n + K]} sin Bg} (cos nz). (33) 
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5. Enforcement of the boundary conditions at ¢ = 0. From (29) 
= [(A,R + A.D)pe”* + (A;R + A,l)p.e”**] cos (nz). 


At ¢ = 0 the function w,, must vanish. Hence the requirement contained in (15) leads 


to 
(A\R + Alp, + (AR + Ap, = 0. 
This can also be written as 
—a,A, + B, A, — a,A3 + BA, = 
From (16) and (33) one obtains 
{—a,K + 6,[(1 + »)n — K]}Ai + {8,K + a[(1 + »)n — K]}A, 


+ fa.K — 8,[(1 + »)n + K]}A; + {-8.K — a.[(1 + yn + K]}A, = 0. 
To satisfy the conditions expressed by (17), u,, must be calculated: 


= (1/2K*)(A,R + A,J){K + i[(1 + »)n + K]}p.e” sin (nz) 


(34) 


(35) 


(36) 


+ (1/2K*)(A,R + A,I){—K — i[(1 + dn — K]}p.e”” sin (nz). 


Evaluation of this expression for g = 0 results in 


{—a,K — 6,[((1 +n + K]}A, + {8,K —a,[(1 +n + K]}A, 


+ {a.K + 6.[(1 + »)n — K]}A; + {-—8.K + a,[(1 + on — K]}A, = 0. 


From these conditions A, , A, , and A; can be expressed in terms of A, : 


M- Sea=t 2. 


442 at ey % pee 4) 
si (2. + bs), 
obi ee bs) 14 ‘ 


In the calculation of the shear, w,,,, is needed. From (29) one obtains 


W,ece = (AiR + A,D)pie”*® + (Ask + Aal)p2e”**] cos (nz). 
At ¢ = 0 this equation becomes 
Wooo = (AiR + A2l)pi + (AR + A4lD)p2] cos (nz). 
Evaluation of the operations indicated in (41) and substitution in (18) result in 


V = —(D/a’){(a? — 3a,83)A1 + (81 — 38.0%) A2 


+ (a3 — 3a282)A; + (63 — 38.02) A,}. 


(37) 


(38) 


(39) 


(40) 


(41) 


(42) 
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If the expressions given in (38), (39), and (40) are substituted into (42), manipulations 
yield 


OF ins 4Kn(D/a?)( 2+ #) 4, ; (43) 


a, — Bz 
This can be solved for A, : 
(44) 








ee (3 —&) 
A+ = TD 4Kn \oi + Bi!" 


The other three constants can also be expressed directly in terms of the applied load: 








5 ng a ( ne bs) 

a * -E a a? + B/’ 7 
mp a (+4) 

A, = D 4Kn \o? + B3/’ (46) 
_. ae (+ bs) 

i, = D 4Kn \az + B3/" (47) 


The expressions for the displacements u, v, and w are then as follows: 


Va’ ” Bieeaia i 
u= 8K'nD ‘3 +B {[28,K + n(l + »)(a, + 8,)] cos Big 


+ [2a,K + n(1 + »)(a, — B,)] sin Bg} 


—a2:3¢ 


e 


a: + B 





+ {[—28.K — n(1 + »)(az — B2)] cos Bop 


+ [—2a,K + n(1 + »)(@2 + B;)] sin ace} sin (nz), (48) 


Va" - ass rt 
v= Pras fe “**{n(1 + v) cos Big + [n(1 + vy) — 2K] sin Big} 


+ e-*"{—n(1 + ») cos By + [n(1 + ») + 2K] sin B.g}} cos(nz), (49) 





w= ear ‘$ : B {—(a, — B,) cos Big + (a, + B:) sin Big} 
+ JT H {(az + 82) cos Bap + (a2 — Bz) sin ae} cos (nz). (50) 


6. Constant radial force over segment of generator. The results obtained can be 
used in the calculation of the displacements and stresses of cylinders subjected to 
symmetric radial forces of all kinds. For instance, when the inward radial load is stipu- 


lated as 
(P*/25) when |z*|< 6 andas 0 when (L/2)>|2*|>6 (51) 
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it can be represented by the Fourier series 


> P*™ cos (mrax/L) 
=. oo 
with the coefficient P*‘”’ defined as 


p*'™ = (2P*/mm5) sin (mri/L) m= 1,3, -°> (52) 


The displacements u, v, and w can then be calculated from (48), (49), and (50), re- 
spectively, if V is replaced by (1/2)P*‘” and the resulting expressions are summed 
over m. Naturally the summation indicated in (13) must also be carried out if the first 
term alone does not give sufficient accuracy. However, this is seldom the case. The 
summation causes no difficulties because the right-hand members of (13) are geometric 
series. 

If the coefficient b,, is defined as 


b,, = (P*aL/4x’ DK 8)(1/m’) sin (mrb/L) (53) 


m 


substitutions and simplifications yield the following expressions for the mth components 
of the moments and the membrane stresses: 


su 


MS” /{(— D/a)nb,, cos nz] = -3 rary {{[—2vKa, + n(1 — v)(a, — B,)] cos Bie 
a TT Pi 


+ [2vKB, — n(l — v)(a, + B,)] sin Big} 


+ 3 +e {[—2vKa, — n(1 — v)(az + B.)] cos Be 

+ [2xKB, — n(l — v)(az — B.)] sin Bw}, (54) 
MS” /[(— D/a)nb,, os nx] = ; CB {{[—2Ka, — n(1 — v)(a, — B,)] cos Bie 

aT Ki 


+ [2K6, + n(1 — v»)(a, + B,)] sin Big} 


é as¢ j 
+ -3 s {[—2Ka, + n(1 — v)(@, + B,)] cos Bay 
a> + By 
+ [2K6, + n(1 — v)(a, — B.)] sin B29}, (55) 
M!™/[((D/a)(1 — v)nb,, sin nz] e ***[—cos Big + sin Big] 
+e “**[cos Bao + sin B,¢], (56) 
oS” /[(E/2K*)nb,, cos na] 
” is . ‘ , an ‘ 
=-5 = {[na, + (n 2K)8,] cos Bie + [(n + 2A)a, — nB,] sin B,¢} 
a T 3; 
ibs cr ' Ok 9k i 57) 
4. — = {[—na, + (n — 2K)8,] cos Bp + [(n — 2K)az + nf] sin 6,9}, (57) 


Qa + Bb; 
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-ai¢ 
é 


a, /\(E 2K’)n’ b,, cos nx] = -3——s3 {—(a + B,) cos Bip + (—a, + 8,) sin Biv} 
A, + Bo 
+ -5=—s3 {(a: — 62) cos Bap — (a2 + 62) sin Bp}, (58) 
az + B2 
r\™ /\(E/2K*)nb,, sin nx] = —e *** {cos Big + sin Big} 


+ e“** {cos By — sin B.9}. (59) 
When the load is concentrated at x = 0, the distance 6 approaches zero. In that case 


Pp = (2P*/L), (60) 


and 


b,, = (P*aL/4xDK)(1/m). (61) 


7. Sinusoidally distributed radial force antisymmetric with respect tox = 0. The 
load is applied along the generator g = 0 in the radial inward direction and is 2Vsin nz 


with n = mx(a/L) and m an even integer. The solution is assumed in the form of the 
series (13) with 
w= Ae’® sin nz, (62a) 
u = Be’* cos nz, (62b) 
v = Ce’* sin nz, (62c) 


where again A and n are assumed real while B, C, and p may be complex. The boundary 
conditions at x = +(1/2)(L/a) and at ¢ = 0 remain unchanged except that 


(D/a*)w..,, = V sin nz. (63) 


Solution of this problem leads to the same expressions for the constants A, to A, as 
obtained earlier and presented in (44) to (47). Similarly the expressions given in (49) 
and (50) for the displacements v and w are unchanged except for the replacement of the 
factors (cos nx) by (sin nx); also in expression (48) for u the factor (sin nx) must be 
replaced by (—cos nz). 

8. Loading by longitudinal moment. When the radial inward load along the gen- 


erator ¢ = 0 is stipulated as 


M*/5° when 0 < z* < 3, 
—M*/s? when -—-5< 2* < 0, (64) 
0 when (L/2) > | z* | > 4, 


where M* is the applied longitudinal moment to which the line loading is statically 
equivalent, the loading can be represented by the Fourier series 


>> Ms” sin (mrax/L). 


m=2,4,°** 
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The coefficient M**”’ is defined as 
M*'” = (4M*/mr8)[1 — cos (mri/L)] m= 2,4,---. (65) 


The deflections u, v, and w can be calculated from (48), (49), and (50), respectively, if 
V is replaced by (1/2)M*'”’, the changes in sign and in the trigonometric expressions 
stated in the preceding article are carried out, and the resulting expressions are summed 
over m. Naturally the summation indicated in (13) must also be carried out if the first 
term alone does not give sufficient accuracy. However, this is seldom the case. 

If the coefficient b,, is now defined as 


b,, = (M*taL /2x°DK&)(1/m*)[1 — cos (mri/L)] m = 2,4, --- (66) 


the mth components of the moments and of the membrane stresses are again given by 
equations (54) to (59) provided all factors (cos nz) are replaced by (sin nx) and all 
factors (sin nx) by (—cos nz). 


When the moment is concentrated at x = 0, the distance 6 approaches zero. In that 
case 
M\™ = (24M*/L’)m (67) 
and 
b, = (M*/4DK)(a/L). (68) 


9. Sinusoidally distributed circumferential moment symmetric with respect to x = 0. 
This problem differs from the symmetric sinusoidally distributed force problem only 
in the boundary conditions at ¢ = 0. The solution is assumed in the form 


Wier = We) — w(2r — ¢) + W(2r + ¢) — W(i4r —G) + w4rt+ot+::-, (69a) 
Uroe = UG) — U(2r — g) + Ul2er + ¢) — ul4er — go) + U4r +e) —:::, (69b) 


Vior = V(g) + v(2r — ¢) + 0(2r + ¢) + v(4r —g) + V(4r+9)+::-, (69c) 


wy) = Ae’* cos nz, (70a) 
u(y) = Be’* sin nz, (70b) 
v(v) = Ce’* cos nz, (70c) 


where the dependence on z is not shown but is understood, and A and n are assumed 
real while B, C, and p may be complex. 

Along the generator ¢ = 0 an external moment 2M cos nz is applied with n = mz(a/L) 
and m an odd integer. The boundary conditions are 


w=u=0,=0 when g=0 (71) 
Mo = M cos (mrax/L) (72) 


and M,, is the internal circumferential moment per unit length of the shell along the 


cut at ¢ = 0. 
It follows from (71), (31), (82), and (33) that 


A, + A; = 0, (73) 
2A,K + A,[(1 + »)n + K] — A,[(1 + »)n — K] = 0, (74) 
2AwK + A,|(1 + »)n + K] — A.l(1 + )n — vK] = 0. 


| 
“~~ 
~I 
or 
ww 
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Solution of (73) to (75) yields 


A, = A; = A, = —A,. (76) 
The fourth boundary condition (72) along g = 0 gives 
A, = —(1/4x)(ML/KD)(1/m). (77) 
Hence the displacements can be written as 
w = —(1/4r)(ML/KD)(1/m)[e-*** (cos B,g — sin B,¢) 


— e ***(cos By + sin B.¢)] cos nz, (78) 
u = (1/8x)(ML/K*D)(1/m){e~*** {(1 + »)n cos Bye + [(1 + »)n + 2K] sin 6,9} 
—e “**{(1 + y)n cos By — [(1 + »)n — 2K] sin B.~}} sin nz, (79) 
vy = (1/84)(ML/K*D)(1/mn){e~*** {[(1 + »)n(a. — B:) + 28,K] cos Bip 
+ [(1 + »)n(@, + Bi) — 2a,K] sin £,¢} 
— e“**{[(1 + »)n(@2 + B2) + 282K] cos Bay 


— [(1 + »)n(ez — B2) + 2a.K] sin B.y}} cosnz, (80) 


10. Constant circumferential moment over segment of generator. When the external 
circumferential moment is distributed according to the rule 


M*/25 when is?is8 
(81) 
0 when (L/2) > |a*|>6 
it can be represented by the Fourier series 
>> M+” cos (mraz/L) 
m=1,3,°** 
where the coefficient mMs™ is 
M*‘” = (2M*/x6)(1/m) sin (mri/LZ). (82) 


The displacements w, u, and v can be calculated from (78), (79), and (80) if M is re- 
placed by (1/2)M*‘” and the resulting expressions are summed over m. Naturally the 
summation indicated in (69) must also be carried out if the first term alone does not give 
sufficient accuracy. However, this is seldom the case. 

If the coefficient b,, is defined as 


b,, = —(M*/4n’)(L/DK8)(1/m’) sin (mxé/LD), (83) 
the moments and the membrane stresses in the shell can be given in the form 


M\™ /[(—D/a)b,n cos nz] = e *** {[—(1 — »)n + 27K] cos Bp + (1 — »)nsin Bie} 


+ e***{[(1 — »)n + 20K] cos By + (1 — »)nsin By}, (84) 
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MS” /[(—D/a)b,n cos nz] = e~***{[(1 — vn + 2K] cos Bye — (1 — v)nsin B,¢} 
+e “**{[—(1 — v)n + 2K] cos Bw — (1 — »)nsin £9}, (85) 
MS" /{(—D/a)(1 — v)b,,n sin nz] = e7**° { —(a, + B,) cos Big + (a, — B;) sin Big} 
+ e~*** f(a, — Bz) cos Bap + (a2 + Bz) sin By}, (86) 
oS” /|(E/2K’)b,,n cos nz] = —e~*** {n cos Big + (n + 2K) sin B,¢} 
+e” ***{n cos Bag — (n — 2K) sin B.¢}, (87) 
oS” /[E/2K’)b,,.n* cos nx] = e *'* {cos By + sin Big} 
+ e***{— cos Be + sin Be}, (88) 


so AE 2K’)b,,n sin nz] = e “** f(a, — B,) cos Bie + (a, + 8,) sin Big} 


Tze 


+e “**{—(a, + B.) cos Bog + (a, — 8) sin B29}. (89) 
When the moment is concentrated at z = 0, the distance 6 approaches zero. In 
that case 
M‘\™ = (2M*/L) (90) 
and 
b, = —(M*#/4rDK)(1/m). (91) 


11. Convergence. When 6 is finite in the series representing the loads, all the dis- 
placements, moments, and stresses converge. This is true even when 6 = 0, x ¥ 0, and 
the load is a concentrated radial force P*. In this case, when m is sufficiently large, the 


quantities of practical interest behave at ¢ = 0 in the following manner: 


w ~ (1/m’*) cos (mrax/L), w., ~ (1/m’*) sin (mrax/L), 
M, ~ (1/m) cos (mraz/L), M, ~ (1/m) cos (mzaz/L), (92) 
a. ~ (1/m’*) cos (mrax/L), ao, ~ (1/m’*) cos (mraz/L), 


tz, = M,, = 0. 


When ¢ = 0 and zx = 0, the series for M, and M, diverge, but the remaining series 
converge. 

When the load is a concentrated longitudinal moment at x = 0, g = 0, the behavior 
at ¢ = 0 for large values of m is: 


w ~ (1/m’) sin (mrazx/L), w,, ~ (1/m) cos (mrazx/L), 
M, ~ sin (mrax/L), M, ~ sin (mrax/L), (93) 
a, ~ (1/m’) sin (mrax/L), a, ~ (1/m’) sin (mzraz/L), 


tT, = M,, = 0 
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The displacements and the membrane stresses again converge, but the moments 
M, and M, are given by divergent series representing zero when x ~ O and an in- 
definitely large quantity when x = 0. The series converge for all values of ¢ ¥ 0. 

Finally, when the load is a concentrated circumferential moment at ¢ = 0, x = 0, 
one obtains at ¢ = 0 for sufficiently large m: 


w= WW. =| 7, ~*~ 66 = 0, 
WV. ~ cos (mrax/L), M, ~ cos (mrazx/L), (94) 
M... ~ sin (mrax/L), Teg ~ (1/m’) sin (mrax/L). 


Here again the series for the moments diverge in consequence of the representation 
of the load by a divergent series. The series for the stresses converge, and all the quanti- 
ties are represented by convergent series when ¢ ¥ 0. 
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PLASTIC FLOW IN A DEEPLY NOTCHED BAR WITH SEMI-CIRCULAR ROOT* 


BY 
ALEXANDER J. WANG 


Brown University 


Summary. The unsteady motion problem of a circular-notched bar pulled in 
tension in plane strain is considered. The theory of perfectly plastic solids is used. Large 
strains are analyzed so that the material can also be considered as plastic-rigid. The 
basic equations governing stress and velocity are integrated independently in the char- 
acteristic plane. The results are used to construct the boundary change in a step-by- 
step manner. The problem is greatly simplified because at each step the new free boundary 
of the plastic region can be approximated by a circle. The final shape of the boundary ° 
of an initially semi-circular notch is presented when plastic flow has reduced the initial 
connection at the root to a line contact between the shanks. 

1. Introduction. We consider the plastic flow in a deeply notched bar under tension 
with semi-circular root. For a deep notch, the plastic flow is localized in the vicinity 
of the root of the notch and the parts remaining elastic prevent appreciable lateral 
contraction, thus allowing us to consider the plastic flow to be in plane strain. The 
present paper follows a series of other papers on the deformation of notched bars under 
tension with V-shaped roots [1]** and rectangular roots [2]. The former is a quasi-steady 
case in which the configuration maintains geometrical similarity. The latter is an un- 
steady case in which the analysis can be readily made by building up some known 
solutions of slip line fields. The present treatment is another attempt at an unsteady 
case. Furthermore, since actual specimens have circular fillets at the corners, the result 
can also be extended to interpret such tension experiments. Thus this solution may 
help to examine more completely the singularities created by sharp corners in the pre~ 
vious cases. 

The present paper limits itself to the case when the plastic region extends only to 
the circular part of the notch, Fig. 1 where CjBC, is the circular part and C,D, CjD’ 
are the linear parts of the free boundary. The case when the plastic region extends to 


U 


D f ” 
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the linear part of the notch and the case of a circular fillet will be considered in later 
work. 

We shall use plastic rigid theory. It is valid here as well as in most metal forming 
problems, since plastic flow needs to be analyzed in the regions of large flow only. We 
shall also consider the material to be perfectly plastic, i.e., plastic flow occurs at a con- 
stant stress limit. It is not equivalent to neglecting work-hardening, but rather to averag- 
ing its effect over the field of flow. The basic theory of such analysis has been fully 
discussed in the recent literature [3], [4], [5]. Accordingly, only a brief résumé of the 
final equations is given below. 

2. Equations governing plastic flow in plane strain. In the plastic region the stress 
has to satisfy the yield condition and the equilibrium equations. In plane strain these 
three equations lead to a pair of first order hyperbolic equations having two orthogonal 
sets of real characteristics commonly called slip lines. The equations expressed in terms of 
these lines as coordinates (called the canonical form in the theory of partial differential 


=< 








Fic. 2. 


equations) take on a particularly simple form. Referring to the notation of Fig. 2, the 
equilibrium equations become 


p + 2ke = constant along an a-line, 
(1) 
p — 2ke = constant along a #-line, 
and so 
dy 
=, = 0 (2) 
0a OB 
Equation (2) governing the slip line field can also be expressed in terms of the radii 
of curvature R, S of the a and 8 lines respectively, giving the alternative relations 
dS + Rdg = 0 along an a-line, 
(3) 
dR — Sdg = 0 along a £-line. 

The condition of incompressibility and the relation between stress and strain-rate 
lead to a similar set of equations for velocity components. They have the same char- 
acteristics. If u and v denote the velocity components along the a and £-lines, respectively, 
we have 

du — uvudg = 0 along an a-line 
(4) 
dv + udg = 0 along a #-line. 
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We shall use the convention that R, S are positive in the directions of increasing 6 
and a, respectively. If the curvatures of the a, 6 net are as shown in Fig. 2, R will be 


positive and S negative, then 








¢ =a+ 8+ constant, (5) 
p 
_= <2 + 8+ constant, (6) 
as oR = 
a. = —R, dB = §. (7) 
We then have for (3) and (4) 

oR a’s 
“” ia” ©) jaa ~* (8) 

eu av 
da 0B oe da OB is (9) 


In the following, (8) will be referred to as the stress equations and (9) as the velocity 
equations. They are all of the form of the telegraph equation [6]. 
The transformation from the characteristic plane to the physical plane is achieved 


by integrating the following relations: 


Ox Ox 


= Ss = 8 si 
= R cos ¢, 9B sin ¢, 
(10) 
Oy , oy 
=— = Ss —_ = — S 3 2. 
Ao R sin g, ag S cos ¢ 


3. Analytical solutions. Finding the plastic flow inside the plastic region requires 
the solutions of equations (8) and (9) with their respective boundary conditions. The 
solutions of (8) will give the shape of the slip lines at any instant. Since the slip lines 
are the trajectories directed along the plane of maximum shear stress, we can find the 
stress distribution throughout the plastic region. The solutions of (9) will in turn give 
the strain-rate distribution at any instant. The two sets of equations can be solved 
separately in the characteristic plane; thereafter they are transformed back to the 
physical plane. 

When the unsteady motion of the initially semi-circular notched root is investigated, 
it is found that the free-surface of the plastic region maintains very closely the form of 
a circular are. To carry out the complete analysis at all stages in the deformation, it 
is therefore only necessary to evaluate the stress and velocity fields for a free boundary 
having circular form. This basic problem is considered in detail below. 

First of all, we shall determine the constants in equations (5) and (6). By symmetry 
we may consider one half of the bar. Let the free surface of the plastic region subtend 
an angle 26, of a circular arc, and take the origin at the center of the bar, Fig. 1. Equa- 


tion (5) becomes 


e=atB+i. (11) 
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On the free boundary p = —k. This determines the constant in (6), giving 
1 
ap = tt B- 4&5 (12) 
The free boundary BC therefore has the equation 
—atB= 0.. (13) 


On the boundary we also have the relation, 
at+B= #6. (14) 
Next, we look for the solution to the stress Eqs. (8). Because of symmetry it is 
necessary to solve for only one of the terms F# or S. Hence, for S 
as ss 
Saag + 5 = 9: (15) 
The boundary conditions on BC are: 
R= -—S = +2'a. 


Using (7), we have 


as —. mye 

Oa 
Making use of the fact that BC is a circular arc, we find 

os 1 9 

—— «= 3°'6. 

08 
The boundary conditions are therefore, 

, as " as 
S = —2""a, da = —2"a, ap = 2%. (16) 


This is a Cauchy problem since two boundary conditions are given on a curve which 
is nowhere tangent to a characteristic. The solution is obtained by Riemann’s method 
[6]. The Riemann’s function in this case is Jo[2(¢ — a)’"(n — B)'”’]. 

The solution is obtained by considering the general formulation of Green’s theorems 
which in this case reduces to the following 


os ~y Od « — 
§ {J oF a8 + 8 ia ~ 


With the path PTQ as indicated in Fig. 3, we have, 


, as ,. ado 
i « I {J og Ot 85! dah 

















7%. 
= 

~ 
ll 


E+6-¢ 
—2'"a — [ Jo[2i(—E — B + 8.)'"(8B — n)*]2'%a dB 
g fone 1/2 
“[teesea)" 
1-%. g ae 0 
= —2'a{1 +P; a P 1} 


= ~2"al1 +2 > (—)* eet [T,(2B) + Tn.s(2B))}, (17) 





m=0 


; _ 
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where 2B = 6, — » + &. The last step is obtained through the help of Sonine’s integral 
representation of Bessel’s function [7]. 

We know, however, from consideration of symmetry that the slip lines form logarith- 
mic spirals with radii of curvature 


S = —-2'?7 = —2) ae" = —2'"a¢"*. (18) 


The numerical identification of (17) and (18) is readily established, but the reduction 
of (17) to (18) through Bessel function recurrence relations appears to be quite com- 
plicated. 








B 

Q fc 
-0+B=0¢ ; 
} 

Ae -JPlien) 
a 

Cc’ A 

Fia. 3. 
By symmetry, 
RE, 1) = —S(—1n, —8). (18a) 


The pressure distribution in the plastic region can be found easily, being especially 
simple across the line of symmetry AB. On AB, —a@ = 8 so that 


o,= —pt+k (tension) = 2k(1 + 6, + 2a). 
With 
xz = ale” — e””), 
6, =a+tB=2a+t 6, 


6, = log ( _ 2). 
a 


| ee oe. ee 
Qk log (. 2) 2’ 


ai 1 + log (. _ 2) |. (19) 


These relations are plotted in Fig. 4. 


We obtain (see Fig. 4) 


Oy 
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Next, we look at the solution to the velocity equations (9). By symmetry we may 


consider u only. Thus, we have 


92 
ae pe © (20) 




















With the boundary conditions 


u=sing along AC, 
ou 

v= —cosg = - along AC’, (21a) 
0a 


where U has been taken as the unit of velocities. 

The last relation is obtained from Eq. (4). At the point A, we see that there is a 
jump in u across AC’ of magnitude 2'’’. Integrating the last relation we find that, along 
AC", 

u = —(sing — 2'”’). (21b) 
Note that the boundary conditions on the plastic rigid boundaries AC and AC’ depend 
only on ¢ and not on the exact shape of the curve. Therefore, the result can be applied 
for any shape of plastic rigid boundary, the particular discontinuity of 2'”” in u requiring 
the body to be symmetric at least about one of the x and y axes. 

The problem of solving (20) and (21) is a problem of the Riemann type since we 
two characteristics. The Riemann’s function 


are given one condition on each of the 
‘ , tifa ‘72 “io 
is again J,[2(¢ — a) ‘(mn — 8) ’’]. Using 


F “ a. 0 
> E an dB + u oJ da | = () 
d 0g Oa 
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with the path PQAT as indicated in Fig. 5, we have 


F Ou j OJ 5 
upl(é, n) = ulé, O) - | Jo dp — a de 


JQA 0p Jar Oa 


2"? — sin (: + *) + [ Jol2E'*(n — B)'”"] cos (s r ) 6 


“0 


4 
0 1/2 1/2 
+ [ | 2" — sin (« + *) | AS ol2n E —9)"] da. 
t 4 da 


B 








Integrating the second integral by parts and changing the variables, we obtain 
unl, 0) = 212-7] + [ Lol2(—¥)'"] c08 (, +5- z) dz 
+ [* 1l2ne)"*] 00s (¢ + % + 2) ae 


2-/97,(2(—ém)'7)1 + 9 — +277 Do L2[2(—En)"”] 


p= 
n \ [ n £ | 3 | £ n | 
/ ms 1 > Ss = (£) s 1 < M : ‘99 
AE) Lap +1 2p — | \7 +17 Zp — | (22) 
The last step is obtained through the help of Sonine’s integral representation of Bessel’s 


function [7]. The convergence of the series in the range considered is very good. By 


symmetry we have 





In the ste p-b 


useful. On the | 


Hence, the radia 


U » ale’ 


-2! 


v(é, 9) = —u(—y, —8). 


(22a) 


ame 


step construction of the boundary, the following expressions are 


uundary 


t= 1(0 — 8@.), n = 3(0 + 8.) 


\ elocity 


(1+ 6)1,((@ — &)'*] 


vg — ey’ +l (2 + a) (8 — apt) “e (? — 0, 
6 — 0.) \4p° — 1 6+ 86, 


yog2)) es 
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and the tangential velocity 


r -1/2 - 2 2\1/2 0 6.\” 6 — 0.\" 
U, = 2° (ut) = > I2y[(6. — 6) | (2+) ~ (2= 4, | (24) 





Uy 
nd , o_O 
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These quantities are plotted in Fig. 6 for 6, ranging from 0 to 7/2. At the point B, 
(U.)p = 0, 


(U,)p = 2?ug = 10(0.) + 2 D> (—)’Z2p+1(0.)- 


p=0 


Since sinh z = 2 > 05.5 Iny+:(z) (being transformed from the expression for sin z [7]), 


we have 
2"7u, = 1,(6.) + sinh 06. — 4 >> Iu,-1(0.)- 
p=l 


The convergence of the last expression is particularly rapid; for 4 significant figures 
one term in the summation is enough, i.e., 


2 °us = I,(0.) + sinh 6, — 4/;(8.) (25) 
The transformation functions are obtained in the following manner. With 
~y /2 6-—-8+a /2 Oe- 
S = —2ae’ s fafa 


we can integrate Eq. (10) obtaining 


a= 2a / eft *** cog (a +B+ =) da + f(8) 
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or 
= —2'a / er ie (. + 8+ *) dB + g(a) 
or 
a = ae’ *** cos(a + p)+C,. 
Ata = 8 = 0,2 = 0soC, = —ae™. Therefore, 


ax = ae’‘[e*~* cos (a + 8) — 1]. 
(26) 


Similarly, 


y = ale’**** sin (a + 8)]. 


They appear to be quite simple, however, with a, 8 and @, all varying with time, the 
transformation functions for a generic ¢ are very complicated. 

4. Approximate solutions. Since the analytical solutions above are not in closed 
form, the dependence of @, , a, 8 on t becomes very complicated and it is much simpler to 
analyze the unsteady motion by approximate means. If we represent a boundary curve 
in intrinsic coordinates R and S and if we know the normal and tangential velocity com- 
ponents along this curve, then the subsequent shape of the boundary will depend only on 








the normal velocity component. Referring to Fig. 7, we have, after the time At, 


1 a 


~ 55, OU: + vu? | + O(At), (27) 


R, =R+ af v, 
where we denote differentiation with respect to @ by a prime. 

For an originally circular boundary, dR/d@ = 0; U, + U?’ was numerically com- 
puted for @, = 2/4. The result shows that if At_is taken as 1, the deviation of the new 
boundary from a circle with the same original center is within +2.0% and if At is taken 
as 0.1, the deviation is within +0.4%. Furthermore, we can change the position of the 
center of the circle with the resulting error of 0.33% for At = 1 and 0.14% for At = 0.1. 
In the following step-by-step construction, a circle is drawn through the new positions 
of B and C with the center remaining on the z-axis. The validity of approximating the 
new boundary by such a circular are is checked all along the step-by-step process for 
6. ranging from 7/2 and zero. The maximum error is about 0.2%. It is considered to 
lie within our error of graphical construction. 

5. Step-by-step process. After a small time increment Af, an originally circular 
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arc boundary BC will take the form B,C, , Fig. 8. We shall pass a circle through B, 
and C, with center O, on the z-axis. Hence, 


(r + Ar)’ = | cos 6. + (#2) At — an, | + sin 6. + (4x) f, (28) 
P c . 


r+ Ar+ Am =rt+ 2)7uz, Al. (29) 





B2 


a ne a Lugar 








For the values at C, we can make use of the boundary conditions. Then, 


=) , # (: ) 
= Ss _ } _ ) “OS . =— = 
(= u, Sin q é ; v. CO Pr 0 


= cos 0. + sin 8, 


(W) 14 sin¢ 4 
di). ~ sin 6, cos 6, . 


Solving (28) and (29) for Az, and Ar, we obtain 


Q'*uz — 1 — sin 8, 


At = ae ~At+ O(AL), (30) 
1 — cos 6 
—2'*u, cos 86. + 1+ sin 8, ” : 
Ar = At + O(A?), (31) 
| — cos @ 
mm ‘=e . 
2 toe AX» 
r-+ Ar 
Axo ’ 2 "up ' . 
= 9, + (1 — oe) — = at + O(AP). (32) 
r r 


; ais 
Since these algebraic equations are expressed in terms of @ and 2’’“u, only, where 


2'/*1,, is represented by Eq. (25), it is easier to work with these equations than to inter- 
polate the results from Fig. 6. 

An example for an initial 6, = 2/2 is carried out. The free boundary is found after 
a small time increment Aft by using eqns. (30), (31) and (32). The first step is to replace 
6, by 1/2, put 2'up(x/2) = 3.644 in these eqns. and find the position of the new center 
and the values of the new radius and of @., . For this new value of 6,, we can calculate 


the corresponding 2 "uz . The new values of 6. and 2'’’w, are then substituted in the 


three equations to find out the next @, and so forth. The magnitude of the time incre- 
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ment is chosen to provide as rapid a procedure as possible in conjunction with satis- 
factory accuracy. In the present example At = 0.050 with r = 1 at that instant. The 
accuracy of such a step is checked by taking At = 0.025 for two steps and comparing 
the result with that of one step of At = 0.050. The resulting points are indistinguishable. 
However, similar comparison between two steps of At = 0.05 and one step of At = 0.10 
shows that Af = 0.10 yields considerable error. Therefore, the step of At = 0.05 is 
taken and the resulting error is no greater than the unavoidable graphical errors. The 





15 to? 2.149909 L.of jsym. 


t = 1.1490, 






t= 0.5400, 





wm 





{ L.of sym 
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final shape along with two intermediate steps are shown in Fig. 9. From Fig. 8 the 
initial discontinuity in slope is calculated as follows 





eae ties (r + Ar) sin (0, + AO.) — rsin 6, — At sia 
—— ae + (r+ Ar) cos (6. + A@.) — r cos 8." 


For 6. = 2/2 and 6, + Aé, = w/2, sin (0, + Ad.) = 1 and cos (6, + Aé,) = —Ad, 


Ar— Al 
tanv¥ = 
Ary — rAd, 
26 — 1.6 
_ _2.04At — 1.64400 _ a” = .2509, 
(1.6446? 4+ 3.644)at 3.986 
7 = l “hg a. 


Note that the boundary forms a recess at the point C which is also the case for a V-notched 


har. The final angle of 63° 26’ can be found in a similar manner. 
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6. Remarks. We note that the problem is conveniently separated into two mutually 
independent ones, viz., stress and velocity solutions. These are then combined to give 
the solution to the over-all problem. However, the combination is complicated by the 
fact that both a and @ are functions of time. Complicated integral expressions arise and 
so step-by-step method is adopted. Fortunately, the approximation of the new boundary 
by a circular are is good, making the computation reasonably simple. 
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ON THE STABILITY OF SOME FLOWS OF AN IDEAL FLUID WITH 
FREE SURFACES* 


BY 
J. L. FOX AND G. W. MORGAN 
Brown University 


1. Introduction. Steady state plane flows of an incompressible inviscid fluid with 
free surfaces were originally studied by Helmholtz [2] and Kirchoff [3], and have since 
been thoroughly reported in the literature. Their work was an attempt to improve the 
classical solutions of flow around sharp corners which are physically unacceptable because 
they involve infinite velocities at the corners. Helmholtz and Kirchoff reasoned that as 
the velocity becomes large, the pressure in the fluid decreases to the value at which the 
fluid goes over into the vapor state. This gives rise to a so-called cavitated region bounded 
by a “free surface’ over which the pressure is assumed to be maintained constant and 
uniform. 

One also deals with steady free surface flows in the case of jets flowing in an ambient 
constant pressure atmosphere. Several examples of such problems are treated by Milne- 
Thompson [4]. 

Very little work has been done in the past concerning time dependent flows with 
free surfaces. Lord Kelvin [5] discussed the vibrations of a hollow columnar vortex flow. 
Some unsteady free surface flows under the influence of external forces, such as gravity 
waves in water, are discussed in Lamb [6] and the investigation there is extended to 
include the effects of surface tension and viscosity. 

tecently Ablow and Hayes [1] developed a theory of the small perturbations of the 
two-dimensional flow of a perfect fluid in the presence of a free surface without external 
forces. They then used their theory to study two specific problems, namely the flow 
around a hollow vortex and the flow through a Borda mouthpiece. 

The present investigation will concern itself with an extension of the work of Ablow 
and Hayes to some free surface flows of jets as well as to a number of generalizations of 
problems treated in [1]. Our primary concern will be to obtain information concerning 
the stability of these flows. 

2. Resume of basic theory. The basic theory underlying the methods used in this 
report has been discussed in detail in the work of Ablow and Hayes [1]. For the sake of 
convenience however, a brief outline of the important results will be given below. 

A. Assumptions. We shall be dealing with perturbations of steady state flows 
which do not fill the entire plane. They can be conveniently divided into three categories: 
(1) flows which are cavitated due to the fact that there is a minimum pressure the fluid 
can sustain; (2) jet type flows; (3) flows which are a combination of (1) and (2). In all 
cases there exist in the steady flow free surfaces along which the pressure remains constant 
and uniform. 

The fluid is assumed homogeneous, incompressible and inviscid. Both the steady and 
perturbed states are assumed to be irrotational and two-dimensional. 
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All quantities are written in non-dimensional form through the use of a characteristic 
length, pressure and velocity in such a manner as to make the steady state velocity along 
the free surface of unit magnitude. 

Under the assumptions made the flows must satisfy Bernoulli’s equation in the form 

e+ pq + pe = Clb, (2.1) 


where 7p is the pressure, p the density, q the velocity, ¢ the velocity potential and C(t) 
is a function of time alone. The dot indicates partial differentiation with respect to time. 
Furthermore, we can introduce a complex potential f and a complex velocity w 


such that 


= af (2.2) 
dz 
where w = u — iv and wu and v are the Cartesian velocity components in the z = x + iy 


plane. 
B. Basic flow. The basic steady flow satisfies the steady form of (2.1) 


~ 


Po + $PQo = const., 


where the use of a zero subscript denotes the basic steady flow. Since 


df. 
Wo =a = aand = qo = | Wo I; 
we can write (2.3) as 
Po + Fpwows = const., 
where the star indicates the operation of taking the complex conjugate*. 
C. Perturbation relations. We shall now give the steady state basic flow a small 


perturbation in terms of a small real parameter ¢ in the form 


Z2=2 + &(%, bt), 
rz. £) = FolZo) 7 Ello. t 
wz, t) = Woz + €W,(Z ft), 
plz, t) = polo) + epilZo , 2), 


where f and w are analytic functions of z) . All subsequent relations will be linearized 
by neglecting terms of order e and higher. Hence all perturbations are small perturba- 
ure correct only to first order in «. It is convenient to perturb the 


tions in that they : 
although the perturbations in f, w and p are given in terms of 


independent variable Zp , 
the fixed point 2 . 
The perturbations given above are not independent since we can derive the following 


relations from (2.1) and (2.2) 


Woz1 + W, = Wo! , (2.4) 
», + Relw,w* + fi — wei] = 0, (2.5) 
I : 


—— 


*The star is used instead of the usual bar, for typographical reasons. 
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where the prime indicates partial differentiation with respect to f, , and Re denotes the 
real part. Thus we see that only two of the four perturbations are independent. 

We note that, when properly chosen, two different sets of perturbations, e.g., (23 , fs) 
and (z, , f4) may represent the same physical perturbation. Their difference, namely 
2; = 23 — 24, f1 = fs — fa, will leave the flow unchanged and the perturbation (2, , f;) 
will be called an invariant perturbation. 

We define a stationary perturbation (z. , f2) to be one in which any given physical 
perturbation is evaluated at a fixed point 2, of the basic flow, i.e., one for which the 
space variable is not perturbed (z. = 0). We can now find, corresponding to a given 
perturbation (z, , f,;), a unique stationary form by superposing on (z, , f,;) the invariant 
perturbation z,; = —2,. 

Using (2.4) and (2.5) we can derive the following relations between the stationary 
perturbations f, , w. and p, : 


W. = Wof (2.6) 
P2 = —pRe [wowt ft + fo). (2.7) 


In this formulation there is only one independent perturbation quantity, say f. , 
restricted only by the condition that it be admissible under the boundary conditions of 


ad 
29 


the problem. 
In subsequent work, for the sake of compactness, we shall not change the name of a 


function after a change of independent variable, e.g., we shall write 
f(Z) = flwoleo)] = f(wo). 

D. Free surface condition. There are two conditions that must hold on the free 
surface. First, the free surface pressure remains constant, and second, a particle originally 
on the free surface remains on the free surface in the perturbed state. From (2.5) we see 
that we can satisfy the first condition by demanding that 

Re [f + fz + wiz,] = 0 (2.8) 
where we have used the relation between f, and its stationary form f/f, 
fe = fi — Wh - 
The second condition can be shown to imply that 
Im [f3 + (woz,)’ + (woz,)'] = 0 (2.9) 
where Im denotes the imaginary part. 

We can satisfy (2.8) identically if we set the expression in the bracket equal to 

. o . . 7 a . a) . 
ix (Zo , t), where the function x is real on the free surface and otherwise arbitrary. Solving 
for z, to obtain the perturbation in the space co-ordinate which takes the basic flow 
free surface into the perturbed free surface, we have 


“1 


z= + fix — Difel} (2.10) 


where the operator D[ ] = @[ J]/dfo + A[ J]/dt. Now substituting for z, from (2.10) 
in (2.9) we find, after some reduction, that the free surface boundary condition is 
Im {D[f. — wD[f.]] — fr} = 0 (2.11) 


; 


where w = W)/W5 . 
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Adopting the notation 
H = Lif] = Dif. — oD{f.ll — f; (2.12) 
(2.11) becomes with wy as the independent variable 
H(w.) = (H(w»))* on wows = 1. (2.13) 


E. Symmetry considerations. In cases where the basic flow has a velocity distri- 
bution which is a symmetric function of z) (symmetric basic flow) certain simplifications 
can be made in the problem. In functional notation a symmetric function F(v) satisfies 


F'() = (F'e*))* 
while an antisymmetric function satisfies 
F°(v) = —(F*(v*))*. 
We can combine both relations in a convenient notation 
Fy) = +(F"(*))* 
where the first superscript corresponds to the upper sign, and the second superscript to 


the lower sign. 

Certain operations performed on a symmetric or antisymmetric function preserve 
these properties. It can be shown, for example, that the operations of differentiation and 
integration are symmetry preserving, i.e., the symmetry or antisymmetry of the function 
remains unchanged. Also, if 


Fv) = +(F**(v*))* 
and we transform to a new variable u such that 


u(v) = (u(v*))* 


then 
F'(u) = +(F**(u*))* 
i.e., the function retains its symmetry or antisymmetry properties in the wu plane. 

The important consequence of these considerations is embodied in a theorem which 
will be stated here without proof (for proof see [1, pp. 26 et seq.]). 

Symmetry Theorem: If the basic flow is symmetric, any perturbation can be repre- 
sented as the sum of a symmetric and an antisymmetric perturbation each of which 
satisfies all boundary conditions and so is an admissible perturbation in its own right. 

F. Separation of time dependence. At this point in the development of the theory, 
the only restriction placed on the perturbation in the potential is that it shall satisfy all 
applicable boundary conditions. We shall attack the problem by assuming solutions of 


the form 
fo = Gilwe® + G,(wi)e". (2.14) 


We anticipate that this choice of time dependence will lead to an eigenvalue problem for 
the determination of the functions G, and G, . We expect to find that the boundary con- 
ditions can be satisfied only for certain specific values of \ and that a general solution 
will be a sum of all such elementary forms of f, . Our primary concern will be to determine 
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the magnitude of the real part of all admissible \ since this will indicate the stability of 
the flow. We shall have 


unstable perturbations for Re{A} > 0, 
neutrally stable perturbations for Re{rA} = 0, 
stable perturbations for Re{r} < 0. 


It might seem sufficient to assume f, in the form 
fe = Gwe“ (2.14a) 
since if \ and A* were both eigenvalues, both would be found among the admissible 
values of \. The form (2.14) has been chosen because it is found that the elementary 
form (2.14a) is not capable of satisfying all the boundary conditions, whereas the form 
(2.14) can represent an admissible perturbation. 
Substitution of the form for f, from (2.14) in our previous expression (2.12) for the 


operator H gives 


H(w.) = L,[G,Je! + L*[G.]Je"', (2.15) 
where 
we G ; dG 
IG] = re as = T(w’ »). , 
L(G] > a + 2du in + AG@’ + Aw) (2.16) 


Our free surface boundary condition (2.13) becomes 
L,[G,(wo)] = (Ly+[G2(wo)])*. (2.17) 


In the event that the basic flow is symmetric we can decompose f, into symmetric 
and antisymmetric components and write 


fe=fth, (2.18) 


where 


fo? = Gie' + Ge", 
and we can show that 
Gr"(wo) = +(Gi*(w))*. (2.19) 


With (2.19) we can eliminate G, from our free surface boundary condition (2.17) which 


then becomes 


L,(Gi"°(w,)] = «| (2). (2.20) 


This condition is to be applied on the free surface wow} = 1 or wi = 1/wo . We can use 
analytic continuation, however, and demand that it hold over the entire wy plane. 

The remainder of this paper will be devoted to solving the perturbation equation 
(either (2.17) or (2.20) for several types of problems. The analysis gives rise to relations 
which are exceptionally long and cumbersome. Due to limitations of space, only the 
essential features of the solutions will be presented here. For further details the reader 
may refer to the report mentioned in the footnote to the title. The subscript 0 in wp , 
which denotes the basic flow, will be dropped in the work that follows. 
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3. Jet impinging on a finite plate. One of the flows which the authors intended to 
study was that which results when a plate of finite width is placed at right angles to an 
infinite stream which is uniform at infinity and which has an infinite cavity bounded by 
constant pressure surfaces downstream of the plate, (the Helmholtz plate problem). 
Since the velocity is assumed to be the same no matter how one approaches infinity, the 
points corresponding to upstream and downstream infinity in the physical plane map 
into one and the same point in the hodograph (w,) plane. In the analysis of the perturbed 
flow one is then confronted with the problem of applying boundary conditions pertaining 
to different regions in the physical plane at a single point in the hodograph plane. Math- 
ematically the result of this is that the perturbation quantities exhibit a very irregular 
behavior at this point; it is found, in fact, that one is led to a differential equation which 
has an irregular singular point there. 

To overcome this difficulty an attempt was made to differentiate between upstream 
and downstream infinity by artificially separating slightly the source and sink repre- 
senting these points in the hodograph plane. After applying the boundary conditions 
the points would be made to coalesce and it was hoped that the results so obtained would 
be the same as if the original problem had been solved directly. 

In the course of this study it was realized that the separation of upstream and down- 
stream infinity could be accomplished in a more direct and rigorous manner by con- 
sidering the Helmholtz plate as the limiting case of a related and much more general 
physical problem, namely that of a jet of finite width which impinges on a finite plate 
and is thereby divided into two jets which diverge from the edges of the plate and tend 
to become straight jets inclined at +8 to the horizontal as they approach downstream 


=—8 
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infinity (see Fig. 1). The limiting case of infinite jet width then constitutes the proper 


mathematical formulation of the Helmholtz plate problem. 


A. Basic flow equations. We begin by writing the potential of the basic flow in the 


hodograph plane (see Fig. 2). This flow may be thought of as arising from the presence 
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of unit sources at w = +1 and'sinks of strength one-half on the unit circle at w = +a 


and -+-a* where argument of a = 6, the jet inclination at downstream infinity. The 
potential is then found to be 
; (w’ — 1)* 
f.(w) = log ———s— 3.1 
Sow) -(w — a)(w — a*’) (3.1) 
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Fia. 2. 


Since we are only interested in that part of the hodograph plane which corresponds 
to the physical flow, it is convenient to introduce a transformation ¢ = w* such that the 
right half of the unit circle in the w plane goes into the entire unit circle in the ¢ plane. 
This transformation has the advantage of reducing the number of singular points on the 
unit circle. 

The potential in the ¢ plane becomes 

fol3) = log 7 sae (3.2) 
where 
b = a’. 

[he asymptotic inclination of the downstream jets, as characterized by a parameter 
8 = sin 6, can be related to the dimensionless ratio (d/D) of plate width to original jet 
width. This can be accomplished by considering the integral 


»4/2D l ‘d 
, d= (4). 


“0 


In terms of hodograph variables we have 
l (dfo 
dz = —|——] dw 
w \dw 


2" ae = [1 (te) 
7 = : >= 2(d ‘D) 
| dz | ad du 1/2(d/D). 


“0 -0 


Thus 


Now using (3.1) for f,(w) which involves the unknown sink position a, we finally arrive at 


; {+28 
(d/D) = 7/21 —-— V1 = 8°) + 8/2 log ( + 8) (3.3) 
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The limit cases 8 = 1 and 8 = 0 correspond to the impinging of a jet of finite width 
on an infinite plate, and the impinging of an infinite stream upon a finite plate, re- 
spectively. These, in turn, correspond to making the sink position a approach 7 or +1, 
respectively. The steady state form of the potential, (3.1), in these limiting cases, goes 
over into the known form of the potential with (d/D) = 0 or © respectively. 

B. Derivation of the perturbation equation. As discussed in the basic theory we try a 
solution of the form 


fo = Ge“ + Gi" (3.4) 


Our first objective is to narrow the possible choice of functions G, and G, by investigating 
the functional properties which f, must possess in order to satisfy some of the conditions 
of the problem. 

(1) Wall streamline condition. The boundary condition on the imaginary axis (map 
of the plate) in the w plane is that the perturbed flow have no component normal to the 
axis. It is convenient to transform this condition to the real axis. To do this we rotate 
the hodograph plane by putting 7» = tw. The boundary condition will now be satisfied 
if w./w is real on 7 real. Denoting differentiation with respect to 7 by the subscript 7 


and evaluating w./w we have that 


an + 1)(n* + b(n’ + D*) ee 

Qn7{2(n? + b(n? + b*) — (n? + D[(y? + BD) + (n’ + 0*))] fan 

must be real on 7 real. Since the coefficient of f,, in the expression above is itself real for 

real 7, we infer that f,, must be real on 7 real and hence a symmetric function of 7 in 

ln| <1. 

(2) Analyticity and symmetry considerations. Since the basic flow is everywhere 

regular in | 7 | < 1, f2, must also be regular in | 7 | < 1. Upon integrating /,, we then 
find f.(n) is both regular and symmetric in | 7 | < 1, ie., 


fn) = (f.(n*))*. (3.5) 
Since the basic flow is symmetric we may decompose the perturbations into symmetric 
and antisymmetric components, i.e., 


f,=fit+f? where f2°(0 = 4(f2°(¢*))*. 


These relations hold in the ¢ and w planes since the transformation from w to ¢ preserves 
symmetry. 
a 
In the 7 plane (y = i¢?) the above becomes 


l 


fo(n = +(f2°(—*))*, 


where the superscripts still refer to the symmetric and antisymmetric components of 
f2 referred to the ¢ (or w) plane. Using the relation (3.5) we have 


fi"(n) = +f37(—7). 


Thus, a ¢-plane symmetric or antisymmetric perturbation is represented by an even 
or odd function of 7, respectively, and we may write 


f= Lian, f= Dba". 
0 
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Since ¢ = —7’ these relations transform in the ¢ plane to 
f=F(O, fr = ¢?R(D). 

where F* and F” are regular functions of ¢ in | ¢ | < 1. We can also write 

fx" = Gite + Gye" 
where we know that, by virtue of the symmetric basic flow, 

Gy(w) = +(Gi°(w*))*. 
Since ¢ = w” is a symmetric function of w, the above symmetry relation also holds in 
the ¢ plane (see 2E), i.e., 


G2°(9) = £(G"(*))*. 


The four functions G{* and G3° can be replaced by two functions g** by the following 
relations which are found to be consistent with the above considerations. 


Gi = g‘(d), G, = (g°(¢*))* 
—(g"(¢*))* 


where g’ and ¢~*g* are regular functions of ¢ in | ¢ | < 1. 


(3.6) 


G= 90), G 


(3) The free surface boundary condition. For a symmetric basic flow the free surface 
condition in the ¢ plane is 


1G") = Lo] (2) | (3.7) 


where the form of the differential operator L, in the ¢ plane is found from (2.16) and 
(3.1). It is not reproduced here because of space limitations. 
If we substitute the forms (3.6) for Gi* in (3.7), we arrive at the following relations 


win oe 
ws) = W(2) (3.8) 
h(2) = -m(3) (3.9) 
where 
h(s) = Llg"(9)] (3.10) 


The form of the differential operator Z, and our knowledge of the functional behavior 
of g* and g* in| ¢ | < 1 together with (3.10) show that, in | ¢ | < 1, h’(S) is a regular 
function of ¢ and that h°(¢) behaves like ¢* multiplied by a function with a simple pole 
at ¢ = 0. For convenience we shall represent h* and h’ in the following manner 


wi ot F (¢ — 1) | 
mo = at Dal Se as, 


ee (¢ — 1) | 
Wi) = 2 lg — DG — I" 








(3.11) 
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with a, and }, unknown constants. These forms for h* and h* possess the required behavior 
in | ¢ | < 1, and in addition satisfy the functional relations (3.7) termwise. 

Our aim is to find functions g‘*(¢) satisfying the relation (3.10) with h**(¢) having 
the forms given in (3.11). To do this we regard (3.10) as an inhomogeneous differential 
“*(¢). We note that the inhomogeneous term /**(¢) is known in form only 


equation for g 
7°°(¢) will retain some arbitrariness which, for a given problem, 


and hence the solutions ¢ 
ought to be determined by the initial conditions. 

C. Solutions of the perturbation equation Equation (3.10) is a second order linear 
differential equation of the Fuchsian type [7] with regular singular points at ¢ = 0, 1, b, 
b* and ©. We can first find the complementary solution about the origin, which we can 


denote as AK‘ + BR by the standard methods of solution in series. Knowing the 
complementary solutions, the particular integral is readily found by the method of 


variation of parameters, viz 


ai ‘ . ‘ r(l ‘ ‘ 
HOVK di ‘o(C)K '(o) dé 
PI= KY | A K® | ES (3.12) 
W(K’'’, K Wik ok) 
where 
l | l , 
pg) = | —T —_ |n*(x), 
L ¢t¢ 1) 2 (% 0 2¢(¢ b*) 
and W(K'‘’.. K”’) is the Wronskian of the two solutions which can be found from the 
differential equation and is At *(é Bae b)°"(¢ b*)** where A is a known constant. 


is an arbitrary ordinary point of the differential equation which we shall 
This solution will be valid up to the nearest singular point, i.e., within 
will have to be applied not only within, but 


The point £5 
choose as ¢ l 
the unit circle. Since boundary conditions 
also on the unit circle, in particular at the singular points, we must find solutions valid 
at these points. To do this we study the indicial equation and hence the form of the com- 
plementary solution appropriate to the point in question and then match this form of the 
solution with the solution around the origin by the process of analytic continuation. A 
new form of the particular integral may then be found. This method must be continued 
until we obtain solutions near all points at which boundary conditions must be applied. 


To conserve space, we shall write the form of the solution around the origin only. 





It is 
eK + BSK” + Dd Er’Kys: , (3.13) 
where A‘;, BY and / sre unknown constants. The term FL; Ay,, is a particular 
integral when A‘(¢) in (3.12) is replaced by, (see (3.11 
od (—1 
} 
; | ( i) ( b J 
Similarly Z°K%,.. is a particular integral when h’(¢) in (3.12) is replaced by 
. 1; 
4 l | 4 1) | 
4 L(¢ bic — b 


The index R. which is an integer used to designate the solutions, is chosen such that 


E** is the first non-zero E**. 
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D. Application of the boundary conditions. At this stage the symmetric and anti- 
symmetric perturbations found as solutions of the perturbation equation satisfy the 
wall streamline and free surface boundary conditions. We shall now investigate the 
restrictions to be placed on these solutions in order that they may satisfy all remaining 
boundary conditions. We shall first discuss the antisymmetric solutions. 

(1) Anti-symmetric solutions. (a) The edges of the plate (¢ = —1). 

We demand that the inner surfaces of the jets diverging from the edges of the plate 
continue to originate there in the perturbed state. Since z, is the perturbation in the 
space coordinate which takes the basic flow free surface into the perturbed free surface 


(see Eq. (2.10) in 2D), this condition will be satisfied if 


I 


Ww 


7 tix — fe — f2;} = 0 at {= —]. 
Using (3.4) and (3.6) we find this implies g;(— 1) = 0. The solution gz valid near ¢ = —1 
will satisfy this if an equation of the form 


C,Ai + C.Bi = 0 (3.14) 


is satisfied with C, and C, non-zero constants. 

(b) At the point ¢ 0. We have all ready specified the behavior of g*(¢) in| ¢ | < 1. 
We can insure that the solution gg possess this behavior if ¢ ‘7° is a regular function of 
tast—>0 

This boundary condition can be satisfied by equating to zero the terms in the solu- 
tion gz valid near ¢ = 0 which do not have the required } order behavior. This results 
in two further linear homogeneous equations to be satisfied by the unknown constants 


which are of the form 


Br + da,E? = 0 (3.15) 
and 
>, 6.7 = 0, (3.16) 


where a, and 6, are non-zero constants. 

(c) Downstream infinity (¢ = b and ¢ = b*). In the basic flow the jets diverging from 
the edges of the plate asymptotically become straight uniform jets. One can investigate 
the behavior of a straight jet when subjected to small disturbances using the approach 
described in Lamb’s Hydrodynamics (Chapter IX). Such a flow is found to be neutrally 
stable and any disturbance is propagated downstream unchanged with the velocity of 
the jet. We shall demand that our perturbed flow behave like a straight jet as we approach 
downstream infinity. Hence, an observer moving with the asymptotic jet velocity will 
see no change in the velocity perturbation w, . If we express this fact by means of the ma- 
terial derivative of w, we find that the boundary condition demands that lim;_.,(¢ — b)~"g 
and lim ;.4°(¢ — 6*) 9" exist. 

An examination of the solutions valid about b and b* shows them to be derivable 
from each other by a simple interchange of the roles of b and b* and hence only one of 


the conditions has to be applied. 
It is found that the terms of the complementary solution satisfy the boundary con- 
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dition for any value of \. By investigating the behavior of the integrals making up the 
particular integral, it can be shown that the term in (¢ — b)~*gz with the least real part 
of the exponent of (¢ — 6) has a non-zero coefficient and hence the boundary condition 
demands that the real part of the exponent be greater than zero. This implies that 


R+Ref{rA}<l-ni. (3.17) 


where ,,,. is the greatest r for which we must have E* ¥ 0. From (3.16) we see that, 
since E> = 0, r,,,. must be at least as great as one. Hence, the inequality becomes 


R+Re {ar} <0 (3.18) 


(d) Upstream infinity (¢ = 1). The basic flow originates from a source point at 
upstream infinity. Since we do not wish the perturbations to alter the fundamental 
nature of the entire flow we shall demand that the perturbation velocity w, and the 
perturbation pressure p, vanish at upstream infinity. These conditions will be satisfied 
if (¢ — 1)g? and g’, respectively, vanish as ¢ — 1. 

Since the point ¢ = 1 is a regular singular point of the differential equation, however, 
our knowledge of the behavior of the solution at such a point shows that these conditions 
are redundant and only one, say the latter, will be applied. 

We must show that in the neighborhood of ¢ = 1 the expression for gz consists of 
terms which go to zero as ¢ — 1 or that any terms which do not go to zero have a co- 
efficient which is zero. If one examines the solution gz near ¢ = 1, one finds that the 
dominant term from the complementary solution behaves like (¢ — 1)~” while the 
dominant term contributed by the particular integral is either (¢ — 1)~™ or (¢ — 1)7**', 
depending upon the relative magnitude of R and Re {}. As yet, R and Re {\} are re- 
stricted only by the inequality (3.18). We can now proceed to investigate the various 
manners in which the boundary condition at ¢ = 1 may be satisfied depending upon 
which of the following inequalities is applicable: 


(z) 2R+ 1 < Re {-—2\}, 
(72) 2R + 1 = Re {-—2\}, (3.19) 

(722) 2R+1> Re {-2r}. 
Consider, for example, Case (i). The term (¢ — 1)*“*" has a non-zero coefficient; 
hence, the boundary condition demands 2R + 1 > 0 which gives R = 0, 1, 2, --- and 
then the inequality shows Re {A} < —}. In a similar manner one finds in every case 


that Re {A} < 0. 

(2) Symmetric solutions. The application of the boundary conditions to the symmetric 
perturbations proceeds in the same manner as in the antisymmetric ones. We again find 
no admissible perturbations with Re {A} > 0. 

E. Admissible complementary solutions. One may consider the possibility of satisfying 
the boundary conditions with solutions containing only terms of the complementary 
solution. In both the antisymmetric and symmetric case one can readily show, however, 
that there are no such solutions capable of satisfying all the boundary conditions. 

F. Conclusions as to the stability of the basic flow. Although we have not found a one- 
dimensional continuum of eigenvalues \, our present analysis has been able to restrict 
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the possible values of \ to the left half of the complex \ plane. This shows that any 
admissible perturbation has Re {A} < 0 and we conclude that a jet impinging normally 
upon a finite plate gives a neutrally stable or stable flow configuration. 

The results of this analysis are valid for any finite but non-zero ratio of plate width 
to jet width (d/D). By making this ratio as small as we please, we can consider the flow 
of a jet of arbitrarily great width past a finite plate. But this is a proper physical in- 
terpretation of the Helmholtz plate problem since one never has a truly infinite stream in 
reality. From this point of view, the stability conclusions reached in this section apply 
to the Helmholtz plate problem. 

On the other hand, we find a different situation if we make the ratio (d/D) very 
large but still finite. In this case the point £ = —1 continues to be an ordinary point of 
the perturbation differential equation and the boundary condition z, = 0 must always 
be satisfied there. However, if one actually sets d/D = o, the basic flow of this section 
becomes that of a jet impinging on an infinite wall. Now the point £ = —1 is a regular 
singular point corresponding to downstream infinity and the boundary condition z, = 0 
no longer applies. A reexamination of the work of the present section will show that it 
was the presence of this additional boundary condition which precluded the existence 
of any non-trivial unstable perturbations. Thus, as far as stability is concerned, the 
limiting case (d/D) = © corresponding to the impinging of a finite jet on a plate of 
extremely large width does not correspond to the case of a jet impinging on a truly infinite 
wall. The latter problem is equivalent to that of perturbations which are symmetric 
about the vertical axis of two equal and opposite jets impinging upon each other. It is 
considered in Section 4. 


& 





4AI\ 
Ma 


Ze - plane 
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4. Equal and opposite jets. In this section we shall consider the flow of two equal 
and opposite two-dimensional jets which impinge upon each other (Fig. 3). 
The potential of the basic flow in the hodograph plane is readily found as that due 
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to two sources at w = + 1 and two sinks at w = + 7, all of the same strength (Fig. 4). 


The flow possesses a vertical axis of symmetry as well as a horizontal one and the 
perturbation problem is conveniently divided into two problems, that of perturbations 
symmetric and antisymmetric about the vertical axis of symmetry. These correspond 
in term to demanding that w,/w be either real or purely imaginary, respectively, on the 


imaginary axis in the hodograph plane. 








Fia. 4. 


Now if one proceeds as in Section 3 to consider the boundary conditions on w,/w 
mentioned above along with analyticity and symmetry requirements and the free surface 
condition, one is again led to perturbations which are solutions of an inhomogeneous 
differential equation. In terms of a new variable ¢ = w* the equation can be reduced to a 
standard form known as Heun’s equation [10] and [11] possessing regular singular points 
at ¢ = 0,1, —1, and o~. 

Again proceeding as in Section 3 we can apply the remaining physical boundary 
conditions to the solutions of the perturbation equation. It is found that for perturba- 
tions which are symmetric about the vertical axis of symmetry, there exist non-trivial 
solutions with 0 < Re {A} < 1. All other admissible perturbations have Re {A} < 0. 

5. Generalized orifices. Each of the flows to be considered in this section represents 
the draining of an infinite reservoir through an orifice (Fig. 5). The sides of the orifice 
are made up of two semi-infinite planes inclined to each other at an angle of 2 7/n radians, 
where n = 2”, p = 0,1, 2, --- . When p = O the configuration becomes the Borda mouth- 
piece which has been treated by Ablow and Hayes, and hence this section is essentially 
a generalization of that problem. 

The similarity of these flows for various n is seen by examining the potential of the 
basic flow in the hodograph plane. The flow in the physical plane is mapped into a sector 
of the unit circle in the hodograph plane bounded by radii inclined at an angle of + r/n 
radians to the positive real axis (Fig. 6). 

To find the potential for which the circular arc and radii are streamlines, we image 
the sector to cover the entire unit circle. The basic flow potentials for all n are similar 
in that they may be thought of as arising from the presence of a source at the origin and 
sinks on the unit circle at the n roots of unity. Thus, the potential may be written as 

oe (5.1) 


ow) = - 5. 
f (w) log (w” — 1) 
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Ze- plane Ze- plane 
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But if we transform to a new variable { = w” the potential becomes 


1 f 

folS) = 5 log ep? (5.2) 
which provides a single representation for the potential of all the flows considered here. 
We are now in a position to investigage the stability of the flows for all n by a single 
analysis. 

In this problem the wall streamline and free surface conditions lead to an inhomo- 
geneous second order differential equation governing the perturbations which can be 
put into the standard form of the hypergeometric equation having regular singular points 
at ¢ = 0, land ~. 

After writing the solutions to the perturbation equation, the application of the 
boundary conditions, which are similar to those applied in Section 3, shows that all 
admissible perturbations have Re {A} < 0. One difference noted here is that it appears 
possible to satisfy all boundary conditions with an antisymmetric perturbation made 
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We- plane We- plane 
(a) n=] (b) n=2 
We- plane 





up only of terms from the complementary solutions. This is possible if ) takes on the 
values 


) 1/2 N = 0,1, 2, --- 
A= -(2N+1)+ E (2N + »| , 
n = 1, 2, 4, 8, --: 
We note that for NV = 0 and n = 1, \ = —1 + 2! which has a positive real part. 


Thus, for all n but n = 1 we have only stable perturbations, while for n = 1, corre- 
sponding to the Borda mouthpiece flow, there exists an isolated unstable perturbation. 

6. Hollow vertex bounded by cylindrical walls. This problem is a generalization 
of a problem treated by Ablow and Hayes who considered the unbounded flow about a 
hollow vortex. 

The basic flow is a cyclic irrotational motion with circular streamlines bounded on 
the outside by a solid circular wall and on the inside by a concentric circular hollow 
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vortex forming a constant pressure surface. The basic flow potential in the hodograph 
plane is simply 
fw) = 7 log w. (6.1) 
This flow has an antisymmetric velocity distribution and hence we may no longer 
decompose an arbitrary perturbation into symmetric and antisymmetric components 
which satisfy the boundary conditions separately. 
Due to the simplicity of the basic flow potential, however, one can assume the form 


of the perturbation potential as 
fr = G,(w)e™ + G.(w)e', (6.2) 
where the forms for G, and G, may be deduced from a knowledge of the basic flow be- 


havior as 


G,(w) = B, log (w) + > a,w’, 


G,(w) = B, log (w) + >> bw’, 


with B, , B, , a, and b, arbitrary constants. 
Applying the boundary conditions that the circular wall remains a streamline and 
that the perturbation satisfies the free surface condition (2.17) we find an infinite set of 


eigenvalues 


\, = —i(n + [N(n)]'”) 


” 1 — a” 
N(n) = nf + ‘| 


and a is the ratio of the radius of the hollow vortex to that of the wall; A, is purely 
imaginary and we conclude that the basic flow is a neutrally stable configuration. 

Corresponding to each value of n we can now find an elementary perturbation f2, 
and can generate any admissible perturbation by summing the elementary solutions 
for all n. With f, known we can evaluate the perturbation of the free surface from (2.10). 
Aside from the arbitrariness introduced by the function x the perturbation is made up 
of a wave pattern whose components travel at angular velocities equal to 1 + [N(n)]'"/n, 
i.e., they either lead or lag the basic flow by an angular velocity of [N(n)]*/n. 

An interesting analogy with the propagation of gravity waves in water is found when 
the ratio of the depth of fluid between the vortex surface and the walls, to the radius of 
the vortex surface is made very small. At the same time we consider a fixed but arbitrary 
ratio of depth of fluid to wave length y of the disturbance. The propagation velocity of 


the disturbance is found to be approximately 


for all n, where 


2rh 


? = — tanh — 
Y 


c 
2ra 


where h is the depth of fluid between the free surface and the wall. This is identical with 
the velocity of the propagation of a gravity wave in water if (1/a) is replaced by the 
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acceleration of gravity g. The term (1/a) represents approximately the centrifugal 
acceleration throughout the fluid for the case of small depth. 

This problem has also been solved by Lord Kelvin in a different manner and after a 
change of notation the results presented here are seen to agree in detail with those found 


by Kelvin 
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CLOSURE WAVES IN HELICAL COMPRESSION SPRINGS WITH 
INELASTIC COIL IMPACT" 


BY 
J. A. MORRISON? 


Brown University 


Summary. This paper deals with the problem of spring surges taking into account 
coil closure. This may occur in many cases of compression springs subject to impact. 
Inelastic coil on coil impact conditions are assumed. The simple theory of spring surges 
is adopted wherein only the motion of the spring wire parallel to the axis of the spring is 
considered and the assumption made that each element of the spring satisfies the force- 
longitudinal strain relation of the whole spring, before coil closure occurs. The basic 
theory of coil closure with inelastic impact conditions has been given by Lee [1]*. The 
spring is initially at rest and unstrained with one end fixed and the other (the impact 
end) is given an impulsive velocity and then either maintained at this velocity or de- 
celerated at a constant rate. The case of a mass attached to the impact end is also con- 
sidered. Conditions are obtained for which partial or complete closure of the spring 
occurs. 

1. Introduction. We are concerned with the motion of a helical compression spring 
as depicted diagramatically in Fig. 1, the end x = / being fixed. Following [1] we use as 


000 00 , 
MADAMAKAS 





u 


food .. 


Fic. 1. 











an independent space variable the position x along the unstrained spring. The dis- 
placement of an element of the spring in position z is taken as u(z, ¢), so that the position 
of the element at time ¢ is given by (x + u). If f is the compressive force transmitted 
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across the section x, then the equation of motion of an element dz of the spring is 


1 
U,, = ——f., 1 

‘ m** (1) 
where m is the mass of the spring per unit original length. This theory is analogous with 
that of longitudinal waves in a linear elastic material and here f replaces the nominal 
stress and m the density. The nominal compressive strain of an element of the spring is 
¢ = —u, . Before closure occurs there is a linear relation between the stress and the strain 


given by 


f = Ke, (2) 


where £ is a constant depending on the dimensions and material of the spring. 
Substituting in (1) gives the linear wave equation for uw, 
Ure — Coz, = O, (3) 
where ¢, = V E/m is the constant velocity of wave propagation. 
The velocity v of an element of the spring is given by v = u, . The characteristics in 
the (z,/)-plane are the straight lines (x + cot) = constant. The corresponding character- 
istic relations are coe = v = constant, or mcv - f = constant. 


Now let us introduce the dimensionless variables 


where / is the length of the unstrained spring, v, is the initial impulsive velocity of the 
impact end, fo = mceovo and €) = U/t . 
Then the characteristic lines in the (X,7')-plane are 


X + T = constant, (4) 


and the corresponding characteristics relations are 


V + > = constant, or V + © = constant. (5) 


We now consider closure of the coils. This takes place when the approach of sections 
of adjacent coils with the same angular position is equal to the initial separation of the 
coils. This condition is expressed analytically in the form 


u(x) — ux +p) =p-—d, (6) 


where p is the pitch of the coils and d is their thickness. 

Since we are considering in this paper only inelastic coil impact conditions so that 
the coils remain in contact after closure (until separated by the elasticity of the spring) 
and we are concerned only with the average effects around a single coil it is permissible 
to replace the finite difference condition in (6) by a derivative (see [1] for a fuller dis- 
cussion) so that coil closure is to be associated with a fixed upper bound for the com- 
pressive strain ¢e, giving the condition 


2= —u, < (p — d)/p = Enz - (7) 
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or 
&<K (8) 
where K = «,,,,/€) . Hence we have 


@ = X, for =< K 
(9) 
@ = K, for =>K 


since an increase in the compressive force may take place when coil closure has occurred 
but no increase can take place in @ according to (8). 
The relation between = and @ is shown in Fig. 2. Considerations of continuity and 


z 











RR — — — - 


Fia. 2. 


momentum across wave fronts of closure yield the conditions 
| Av | = c-Ae 
(10) 
Af = mc-| Av | 


where A represents the discontinuity of the variable due to the passage of the wave, 
and c is the magnitude of the wave velocity, relative to the unstrained spring. 
In terms of dimensionless variables Eq. (10) becomes 


| AV | = Q-A® 
(11) 


Ar = 2-| AV | 


where 2 = c/c, . Equations (4), (5), (9) and (11) are used in determining the solution 
of the problem for prescribed boundary conditions. 

2. Impulsive velocity at impact end maintained. In this section we consider the 
case in which the end z = 0 is given an impulsive velocity v) at time ¢ = 0 and this 
velocity is maintained, until complete closure of the spring has taken place. The spring 
is taken to be initially unstrained and at rest. Until closure occurs the configuration in 
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the (X,7)-plane is as shown in Fig. 3, X = 0 corresponding to the impact end and 
X = 1 to the fixed end of the spring. Throughout the motion W = V |,.,. = 1. 

A wave of velocity, stress and strain discontinuity is propagated along the spring 
with constant speed c. , commencing at the impact end at time ¢ = 0 and being reflected 
alternatively at the ends. Using Eqs. (4) and (5) we find that 


V=3{1-(-)'}, 2=r (r>0), 


t 
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T#3 
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in region r of the (X,7)-plane. It is seen that after the passage of the discontinuity 
wave front the stress at a certain position in the spring is increased by an amount f, and 
remains unchanged until the wave again passes. The corresponding velocity is alterna- 
tively vp and 0. 

Now closure does not occur so long as 2 < K. Since = = r in region r it follows that 
closure must eventually occur, as is otherwise evident. Further, when closure commences 
it does so from an end of the spring. 

Letn < K < n+ 1(n > O, integer). Then closure commences at 7’ = n and from 
the impact end if n is even, from the fixed end if n is odd. We first consider the case when 
n is even. Since we are assuming inelastic coil impact conditions the coils do not separate 
after closure so that at T = n + 7 a portion of the spring 0 < X < x is moving asa 
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rigid body with uniform speed equal to that of the impact end. It follows that the stress 


throughout this portion is constant. 
Now the coils in front of the closure wave are at rest and the stress is given by = = n. 


Hence in (11) we have | AV | = 1, and Aé@ = (K — n), using (8) and (9). AlsoQ = dx/dr. 
So from (11) we obtain 
2 = 1/(K — n), A> = 1/(K — n). 


Hence the closure wave velocity is constant as shown in Fig. 4(a) (for the case n = 2). 
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The shaded portion depicts coil closure. Complete closure of the spring has taken place 
at r = (K — n), and then T = T” = K. 

The stress at X = x is given by n + AD = {n + 1/(K — n)}. If K is close ton 
this stress is large, but so also is 2. (Note 2 > 1). We have in fact, 


»T=#K 


| Afdt 


II 


Ina 
fo | A> dT = ag mlvy 


0 « T=n 0 


II 


(weight of spring) X (impact velocity). 


When n is odd closure commences from the fixed end. Ahead of the closure wave 
the stress is given by > = n and the velocity is v) , so V = 1. In this case the closed 
portion is at rest, so again we have | AV | = 1 and A® = (K — n). Hence the closure 


wave velocity is given by Q = 1/(K — n), as before, and the stress in the closed portion 


is again {n + 1/(K — n)}. 

The (X,7)-plane is shown in Fig. 4(b) for the case n = 1. 

3. Impact with constant retardation. In this section we consider the case in which 
the end x = 0 is retarded at a uniform rate a, after being given an impulsive velocity 


vo at time t = 0. The spring is again taken to be initially unstrained and at rest. Until 
closure occurs the configuration in the (X,7)-plane is as shown in Fig. 3. During the 
motion the velocity at the impact end is given by W = (1 — AT), where A = al/covo > 0. 
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As in the previous case a wave of velocity, stress and strain discontinuity is propagated 
along the spring with constant speed c, , commencing at time ¢ = 0 from the impact 
end and being reflected alternately at the ends. By making use of Eqs. (4) and (5) we 
find expressions for V and = as given in the Appendix. We find that, if closure does occur, 
then it takes place from an end of the spring and moreover from that end at which 
reflection of the wave of discontinuity has just occurred. 

We next investigate the time at which closure occurs, if at all. Let =* denote the 
value of = at JT = (r — 1) + 0, (r > 1), at the end X = 0 if r is odd and at the end 
X = 1 if r is even. 

Then 

eee = {(2k + 1) - 2k"A} | 
4 k > 0. 
Lhe = +1 )(1 — ka) 


In particular, =} = 1, =} = 2, ie. independent of A. Closure occurs when @ = > = 
K(0 < K < &), Tr if 0 < K < 1, closure commences immediately from the 
impact end and if 1 < K < 2 from the fixed end after the passage of the first wave. 
These results are independent of A, so we now need consider only 2 < K <~o. The 


= 0 was discussed in the previous section and it was shown that closure always 


case A 
took place. 
Case X > 1/2 
Here =* < 2, (r > 1), since 
(242 ~T Livi) = _ (22; +1 — 22 *) - (1 ga 2kX), (k = 1). (12) 

Hence if K > 2 no closure occurs. 
Cas = = <A < 2 integer 

ase 2+) = ‘ 2 ger. 
Here >* increases strictly up to 2%,,, and 2* < Tk,., all r. [See (12)]. 


But 2%,.. = 2(p + 1)(1 — pad). Hence, if K > — + 1)(1 — pa”), no closure occurs. 
If 2 < K < 2(p + 1)(1 — pa), then closure does occur. In this case there is an integer q 
such that 2 < g < (2p + 1) and 2% < K < 2%,, . Hence closure commences a 
T = q + 0, and from the impact end if g is even, from the fixed end if g is odd. 

We now determine under what conditions complete closure of the spring occurs. The 
distance moved by the impact end is given by y = vot — 4at’, whence t = 1/a{u) — 
Vu — 2ay}. Complete closure occurs if, and only if, y ikea the value le... , i.e. 
if, and only if, v5 — 2ale,,,. > 0. This condi*ion may be written in the form 2\K < 1. 


Closure from the fixed end. We consider here the case where closure commences from 
the end X = 1 of the spring at 7 = (2k + 1) + 0. Hence we have 


{(2k + 1) — 2k’A} < K < {(2k + 2) — 2k(k + IA}, and (1 — 2k) > 0. 


The closure wave has a velocity not less than c, and we suppose that the wave front has 
reached X = (1 — x) at T = (2k + 1) + +r. The configuration in the (X,7)-plane is 
shown in Fig. 5 (for the case k = 0), as before the shaded portion depicting coil closure. 

The closed portion of the spring (1 — x) < X < 1 is at rest, so the stress throughout 
is constant and equal to that at X = (1 — x). Now ahead of the closure wave front 
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> = and Eq. (11) gives AS = 2’A®. Hence closure ceases when 2 = 1, the coils 
then barely coming into contact. (Q = dx/dr). 
Just ahead of the closure wave front we have 


D> = 6 = {(2k + 1(1 — Av) — (2h + 2K + 1DA4+ AI — x} 
and 
V = {(1 — Az) — (2k + Idx}. 
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Just after the wave front has passed, 6 = K, V = 0. Substituting in the equation 
| AV | = QA, putting Q = dx/dr and integrating we obtain 


Afx” + 2(2k + l)x7r + 7} — 27+ 2{K — (2k + 1) + 2k(k + lA}x = 0 (13) 


(r=0 when x = 0). 
Equation (13) is that of a conic in the (X,7’)-plane. If k = 0 the conic is a parabola and if 
k > 1 a hyperbola. 


Case k = 0. Closure commences at T = 1+ 0,1 < K < 2. The equation of the parabola 


1S 











Mx + 7)? — 27 + 2K — 1)x = 0. (14) 
This may be written in the form 
fy -2= 0K fy, =m 
l 2V2d V2 4KX V2 
where 
ass a 1 1 


=(x+7), V= (r — x). 
V2 V2 
The parabola is shown in Fig. 6, with K < 2. Its axis is parallel to the V-axis, and its 
vertex A is at U = (2 — K)/2V/2\, V = —(2 — K)?/4 V2 KX. The tangent to the 
parabola at A is parallel to the U-axis and so at A we have dx/dr = 1. If x > latA 
then complete closure occurs, while if x < 1 only partial closure occurs. But at A 


= —_(U — V) = (4— K*/8R. 


— (V2 
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Hence if K* < 4(1 — 2 AK) complete closure of the spring occurs on this closure wave, 
while only partial closure occurs if K? > 4(1 — 2\K). Now 1 < K < 2. Thus we obtain 
the necessary and sufficient conditions for complete closure in the form 
(0<r<3/8 <x <3 
) Pi vente: or (15) 
1 < K < 2{V/4\?7 + 1 — 2a} \ < (4 — K’)/8K 

; 7, 








The case \ = 0 is included here and it is seen that if K = 2 this is the only ease in which 
complete closure can take place. From Eq. (14), {1 — A(x +7)} =+ V1 — 2r\Kx, 


since r = 0 when x = 0. The stress behind the closure wave (l — x < X < 1) is 


~ . KiU=-rAet7} _  KV1— 2Ky _ 
a [K = {1 —_ A(x T)} | {K = /1 — 2nKx} 


As x increases this decreases. We also find that 
dx V1 — 2rKx 


o= = a 
dr {K — V1 — 2AKx} 


The case of partial closure after the first wave has been considered. Here 1 < K < 2, 
K? + 8\K — 4 > 0. Closure ceases when 
(3K — 2)(2 — K) 


a 


(4 on Ee) 
f K ‘ . 


x BK ”* 
It may be shown that the spring then starts opening up at the constant rate c, (relative 
to the unstrained spring). If, moreover, closure does not commence from the impact 
end when the discontinuity wave is reflected there, then no more closures occur and the 
configuration in the (X,7')-plane is as shown in Fig. 7. The condition for no further 
closure is K? — 8\K + 4 < 0, so together with the above conditions we have 


I} <i <2 


and , = : : 
8S\K 


<r < 5/8 JA > 5/8 
or or ) +2 (16) 
; , ; : ‘ K* 
Qin — V4n7—- 1} < K <2 li<K<2, a> 24. 


If K? > 4| 1 — 2\K:‘|, then partial closure occurs from the fixed end after the first wave 
and closure commences from the impact end after the second wave. Complete closure 
may, or may not, eventually occur (according as 2AK < 1, or 2\K > 1, respectively). 
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Fia. 7. 


Closure from the impact end. We consider here the case in which closure commences 
from the end X = 0 of the spring at time 7 = 2k + 0, (k > 0). Hence we have 


2k{1 — (k — 1A} < K < {(2k + 1) — 2k*a}, and (1 — 2k) > 0. 


The closure wave has a velocity not less than co and we suppose that the wave front 
has reached X = x at J = 2k + +. The closed portion of the coil is moving as a rigid 
body with the velocity of the impact end and hence it is being decelerated at a constant 
rate. The configuration in the (X,7)-plane is shown in Figs. 8 and 9 for the cases k = 0, 


k = 1, respectively. 


la—— XK — » 
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Fia. 9. 


By use of Eq. (11) we find x, r are connected by the relation 
Ar’ + 4kAxr + 2{K — 2k(1 — kA)}x — 27 = 0. 
If k = 0 this is the equation of a parabola and if k > 1 that of a hyperbola. We consider 


the case k = 0. 


(17) 
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Since the portion 0 < X < x is moving as a rigid body with uniform retardation a 
we obtain 
> = Zlx-, — Ax — X), 0O< X <x. (18) 


Hence the stress in the closed portion decreases as X decreases. From Eq. (11) we find 
ee — saci 
2Ix-x = & (1 — Ar) = K (1 — 2\Kx), 


using Eq. (17). Hence at a pt. in the closed portion 2 decreases as x increases. Complete 
closure takes place if, and only if, (K*° + 3AK — 1) < 0. We also have the condition 
0 < K < 1. Hence we obtain the necessary and sufficient conditions for complete closure 


in the form 
[A >0 lO<K<1 


‘ d r 9 

: l Jono? LA < sa 1—-K) (1 ) 
O< zk < {V9rX + 4 — 3Q}, O<jA <> 

2 3K 


4. Mass attached to impact end of spring. In this section we consider the case in 
which a mass M is attached to the end x = 0 of the spring. The spring is initially at rest 
and unstrained, the axis of the spring being horizontal, when the attached mass is given 
an impulsive velocity vp at time ¢t = 0. 

The equation of motion for the mass M is 


dw 
= —usS 20) 
_ BS, (2 
d7 
where W = V |,.,, S = = |,-. and u = ml/M. Until closure occurs the configuration in 
the (X,7)-plane is as shown in Fig. 3. A wave of velocity, stress and strain discontinuity 
is propagated along the spring with constant speed c) , commencing at the end x = 0 
at time ¢ = 0 and being reflected alternately at the ends. 
Using Eqs. (4), (5) and (20), we find 
In region 0:> = V = 0 
In region 1:5 = V = e “*‘"* (21) 
e , s 7 u(T - 
In region 2: 2 = 2 ” cosh {u(X — 1)}, 


V = 2e°*"7-” sinh {u(X — 1)}. 


Since, in region 1, 7 > X, we have = < 1. (u > 0.) Hence, if 0 < K < 1 closure com- 
mences from the impact end immediately and if K > 1 no closures occur in0 < T < 1. 


For 1 < T < 2, if 1 < K < 2 closure commences from the fixed end at 7 = 1 and if 
K > 2 no closures occur. If 2 < K < (2 + e ™) closure commences from the gun end 
at T = 2. We consider these three closures in detail. 


—_ 


Closure from fixed end after first wave. 1 < K < 2. Closure commences at T 
We suppose that at 7 = 1 + 7 the closure wave has moved along the spring to - 


PS 
I 
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(1 — x) the portion (1 — x) < X < 1 being closed and at rest. The configuration in the 
(X,7T)-plane is as shown in Fig. 5. Now just ahead of the closure wave front we have 
@=V=2=e%*"*”, Setting 2 = dx/dr we obtain from Eq. (11) 


ee d : oe 
~atren) ico}. 
Integrating, 
Kux + e***” = 1 (22) 


since initially x = 0 = r. Hence 


r= hog {te} - x 

rm 1 — Kux 

Also, 2 = (1 — Kux)/{(K — 1) + Kux} which decreases as x increases. We find 
> = K(1 — Kux)/{(K — 1) + Kux}, for (1 — x) < X < 1. Complete closure occurs 
on this wave if, and only if, (1 — Ku) > {(K — 1) + Ku}, ie. K(1 + 2u) < 2. But 
1 < K < 2. Therefore, 


1<K<2 ‘each 


2 


(2 — K) 
x — eee 
(1 + 2y) 


te E <K< 
If complete closure occurs then the time taken is given by (1 + 7) = 1/u log (1/(1 — 
Ku)) = T’. But during the period 0 < T < T’, W = e-*’. Therefore f{f W dT = K 
so that the gun end has moved through a distance le,,,, , which is correct. 

If complete closure does not occur then closure ceases when 2 = 1, i.e., when x = 
(2 — K)/2Ku = x* and then r = 1/2y {2 log (2/K) + 1 — (2/K)} = 7*. 


Closure from impact end at T = 2k(k > 0). Let us suppose that closure commences (for 
the first time) at the impact end at 7 = 2k. With the usual notation put T = 2k + 1, 
and let X = x be the position of the closure wave front. Let the stress and velocity just 
ahead of the closure wavefront be given by 6 = > = Z(x, rT), V = V(x, 7). Since no 
previous closures have occurred Eq. (3) holds in this region ahead of the closure wave- 


front. Hence we obtain 





oV z= o> 4 
v4 d= _ 0, oy - 
Or Ox OT Ox 


W = W(r) gives the velocity of the impact end and S = = |x.» . From Eq. (11) we find 
wer Xie _ 
{W — Vix, 7)} = 2° {K — 2x, 7}. (23) 
dr 
and 


Dlxex = 2,2) + FW - VO, OI. (24) 
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But, since the portion 0 < X < x is closed and moving as a rigid body, 


> = Vix-, + (x — X) , O< X<x 
dr 


and in particular 


;, dw 
rex TX Ge 


Hence 


on 


V > 
= (,% + w%) + {2G, 9» -% ve, a}. 
dr ( dr ; 


dr 
But, from Eq. (20), dW/dr + uS = 0. Hence 


ZW + wd) + 200 2) — Voe 9) Xh = 0. 
\ 


dr ' dr 
Let F(x, 7) satisfy the conditions 


oF = pz oF = pV, F(O, 0) = W(0), 
Ox 


—_— : 
OT 


in virtue of the equation dV /dr + d2/dx = 0. 
Then we obtain on integration 


W(1 + ux) = Fix, 7). (25) 


Substituting for W from Eq. (23) we find that 





(1 + ux){K — 2x, 7} x + {(1 + ux): Vix, 7) — F(x, 7} = 0. (26) 
i 
But 
0 ee 0 = . 
- (1 + ux) {K — 2x, D}] — — {0 + wx): Vix, ») — Fx, 7} 
T Ox 
FP» oF OV 
“+ ao =a + 35y i + UX) a 


_ (oF Pec (22 eV) - 
- (# Vv) (1 + ux) ar + ax = 0. 


Hence the differential equation integrates immediately. 

Case k = 0.0 < K < 1. Closure commences from the impact end immediately. The 
configuration in the (X,7')-plane is as shown in Fig. 8. Here =(x, r) = V(x, r) = 0, 
and W(0) = 1. Therefore F(x, r) = 1. Hence, from Eq. (25), W(1 + ux) = 1 and from 


Eq. (26), K(1 + ux) dx/dr = 1. Therefore Kx(2 + ux) = 27, since x = 0 whens = 0. 
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This is the equation of a parabola whose axis is parallel to the T-axis. We find that 








Q= dx = See and S = - 
~ dr K+ px) ' ‘ K(1 + ux)” 
Also we have 
(1 + wX) 


ten-——4, OS X Sx. 
K(1 + ux)’ ' 


Hence in the closed portion 2 decreases as X decreases so that the coils remain in contact 
provided S > K. As x increases S decreases so complete closure takes place if, and only if, 


.4.>K, i« Kee St. tt O<K $1. 
Therefore, 
fu >0 0<K<1 
7 or 
lo<K<(a+™” lo<u< (kK — 1). 


If this condition is not satisfied the coils start to open up from the end X = 0 when 


] 7-2/3 
x= -(K — 1) = x* 
bu 
and then 
= K iio a ¢* 
T= Qu (K 1) = fT. 


Casek=1. 2<K< (2+e€™). 
The configuration in the (X,7)-plane is as shown in Fig. 9. Here 


L(x, 7) = 2e*"*” cosh {u(1 — x}, 


Vix, r) = —2e°**” sinh {u(1 — x)} and W(0) =e”. 
Hence, from Eq. (25), 
W(1L + ux) = 2e*"*” cosh {u(l — x)} — 1. 
Also, from Eq. (26), we obtain on integration 
Kux(2 + ux) + Qur = 4{1 — &*"*"? [(1 + ux) sinh {u(1 — x)} + cosh {u(1 — x)}]}- 


5. Discussion of results obtained. This paper has been devoted to investigating 
the occurence of coil closure and the stresses during closure in a helical compression spring 
one end of which is fixed, the other being given an impulsive velocity when the spring is 
initially unstrained and at rest. The number K = 1/vp (Co€max) arises, where vp is the 
impact velocity and Coénax is a quantity depending on the dimensions and material of 
the spring. Inelastic coil on coil impact conditions have been assumed so that the coils 
remain in contact after closure, until separated by the elasticity of the spring. The 
results may, however, give a good indication as to the actual state of affairs in which 


bouncing of the coils occurs. 
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In the first instance the case in which the impulsive velocity at the impact end was 
maintained was considered. Here complete closure inevitably occurs and it was found 
that exceptionally high stresses could be set up during closure if the coils were very close 
together prior to the passage of the closure wave. Closure commences from either the 
impact or the fixed end and complete closure occurs on that wave. 

Secondly, the case in which the impact end was retarded at a constant rate a was 
considered. The number \ = (a/vo)(l/co) occurs, where l/c, is a quantity depending on the 
dimensions and the material of the spring. The criterion for complete closure is 2\K < 1, 
so that if complete closure of the spring is to be avoided we need (2a/v,)le,,.. > 1. Note 
that le... depends on the dimensions of the spring only and not on the material. Thus 
it is seen that, unlike the first case, complete closure does not necessarily occur and for a 
given spring and prescribed impact velocity it will not occur if the retardation of the 
impact end is sufficiently large. Further, it is found that it is possible that no closure 
whatsoever will occur but that if closure does occur it commences from an end of the 
spring. If 0 < K < 1 closure commences immediately from the impact end, whilst if 
1 < K < 2 closure commences from the fixed end after the passage of the first dis- 
continuity wave, and these statements are independent of the value of \. For K > 2, 
no closure occurs if \ > 1/2. If 1 < K < 2 and only partial closure occurs from the 
fixed end after the first wave, then, if closure does not commence from the impact end 
after the second wave, no more closures occur at all. However, if closure does start from 
the impact end after the second wave complete closure of the spring does not necessarily 
occur. 

Lastly, the case in which there was a mass attached to the impact end was con- 
sidered. The number y, which is equal to the ratio of the total weight of the spring to 
the weight of the attached mass, occurs. This will, in general, be small. If0 < K < 1 
closure commences immediately from the impact end and unless K is fairly close to 1 
complete closure will occur on this wave. If 1 < K < 2 closure commences from the 
fixed end after the first wave and complete closure will occur on this wave unless K is 
fairly close to 2 (u being small). However, if 1 > 1/2, complete closure will not occur 
on this wave, however close K is to 1. It is noted that u = O in this case and \ = O in 
the second case both correspond to the first case. 

We will now discuss the stresses in the spring. Before closure occurs the stress, or 
nominal compressive force, is proportional to the nominal compressive strain and less 
than the quantity mc,v.K. When, however, coils of the spring come into contact, and 
remain so, the compressive force is equal to mcov)= where = > K. In the first case the 
compressive force in the closed portion of the spring remains constant while closure is 
taking place and, as previously mentioned, we may have a large compressive force. The 
time taken for the closure wave to traverse the length of the spring is correspondingly 
small, though. In the other two cases the value of = does not remain constant in the 
closed portion of the spring while closure is taking place. When closure commences from 
the fixed end = is the same throughout the closed portion at any given instant, but it 
decreases as the closure wavefront passes only the spring. When closure commences 
from the impact end the value of = at a given position in the closed portion of the spring 
again decreases as the closure wavefront passes along the spring. In addition, however, 
the value of = decreases uniformly as we pass along the spring from just behind the 
closure wavefront to the impact end, this uniform rate of decrease being independent 
of the time in the second case. 
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APPENDIX 
\In region 27, (r = 0) 
V = 2rx(X — 1) 
> = 2rf{(1 — AT) + nA} 
In region (27 + 1), (r > 0) 
V = {(1 — AT) + (27 + 1)AX} 
>= {(2r + 1)(1 — AT) + AX $F 2r(r + DA} 
(At T=2%+7, (k>0), O<7r<1, OF XKEr 
V = {(1 — Ar) — 2kA + (2k + I)AX} 
> = {(2k + 1)(1 — Ar) — 2A + AX} 
At T= 2k+ 17, (& > ©), 0< + < i, 7<X¥fS1 


V = 2k(X — 1) 
> = 2k(1 —AD — AP} 

At T=(hk+)+7, (k>0, O<r<1, OFS XK<(1-7 
V = {(1 — At) + (2k + DMX — 1)} 
S = {(2k + 1)(1 — Av) — (2k? + QW + AF AX} 

‘i T=(k+l)+r7r, (k>0, O<rs1, (l-7)<XS!1 
V = 2k + MX — 1) 


[ > = Ak + 1)f{(1 — Ay) — AA} 
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ON FINITE TWISTING AND BENDING OF CIRCULAR RING SECTOR PLATES 
AND SHALLOW HELICOIDAL SHELLS* 
BY 
ERIC REISSNER 
Massachusetts Institute of Technology 
1. Introduction. In the following we consider a thin circular ring sector plate under 


the action of two equal and opposite forces perpendicular to the plane of the plate, 
along the axis through the center of the ring (Fig. 1). The ring sector plate may be con- 


he 











sidered as part of a winding of a close-coiled helical spring. The problem of the stress 
distribution in the twisted ring sector with rectangular cross section has first been 
considered by J. H. Michell’ as a problem of three dimensional elasticity. A solution of 
this problem by means of the theory of thin plates will bear the same relation to Michell’s 
solution as Kelvin and Tait’s solution for the torsion of a rectangular plate bears to St. 
Venant’s solution for the torsion of a beam with rectangular cross section. 

The reason for the present note is the further observation that by treating the problem 
as a plate problem it becomes possible to analyse non-linear effects in a relatively simple 
manner by making use of the equations for finite deflections of thin plates. 

In addition to the problem of non-linear effects for an originally flat plate we also 
consider the corresponding problem for a shallow helicoidal shell, thereby obtaining 
information concerning the influence of initial pitch on stresses and deformations in the 
spring. 

*Received Aug. 11, 1952; revised manuscript received March 23, 1953. The present paper is a report 
on work done under the sponsorship of the Office of Naval Research under Contract N5-ori-07834 


with Massachusetts Institute of Technology. 
1J, H. Michell, Proc. London Math. Soc., 31, 140-141 (1889). 
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In the consideration of both these problems one is led to a study of the simultaneous 
action of the pair of axial forces described above and of a pair of couples the axes of 
which coincide with the direction of these forces (Fig. 1). We determine in particular of 
what magnitude the couples have to be in order to prevent the association of circum- 
ferential displacements of the points of the ring plate with the axial displacements 
caused by the axial forces. We find that if these circumferential displacements are 
prevented the effect of non-linearity and the influence of initial pitch are much more 
pronounced than in the absence of the couples. We also encounter a problem of in- 
stability of the ring sector plate which is associated with the presence of the couples. 

A further result which we find is to the effect that non-linearity is relatively more 
pronounced in regard to the magnitude of stresses than in regard to the axial force- 
displacement relation. A similar result is known for the problem of finite torsion of thin 
rectangular plates’. 

It seems worth noting that the present developments may also be considered as a 
contribution to a non-linear theory of dislocations. 

Apart from the question of method the present treatment of the problems of the 
ring sector differs from that by means of the theory of thin rods’ in the following respects. 

Account is taken of a non-linear effect which is significant if the ratio of width of 
cross section to thickness of cross section is sufficiently large. For the helicoidal shell 
this effect is found to be present even in the range of applicability of the linearized 
theory. In the theory of torsion of rectangular plates this effect manifests itself through 
axial normal stresses proportional to the square of the angle of twist. These normal 
stresses have no axial resultant or bending couple but give rise to a twisting couple which 
is proportional to the cube of the angle of twist”. A corresponding result is here obtained 
for the problem of the ring sector. 

While the results of the theory of thin rods are based on the assumption that the 
ratio of width of cross section to radius of center line of cross section is sufficiently 
small the present results hold for all possible values of this ratio. 

2. Differential equations and boundary conditions of the problem. 
equations for finite bending of thin plates of uniform thickness are in polar coordinate 


The differential 


form* 
+22 (wi 98) +23 (w..), a) 
vor =m(()T #428}, a 


where w is the transverse deflection, F Airy’s stress function, and D = Eh*/12(1 — v’). 


28. Timoshenko, Strength of materials, Part 2, Van Nostrand, p. 301; A. E. Green, Proc. Roy. Soc. 
London (A) 154, 430-455 (1936); 161, 197-220 (1937). 

%See A. E. H. Love, Treatise on the mathematical theory of elasticity, 4th Ed., Cambridge 1934, pp. 
414-417 for reference to work by Kirchhoff, Kelvin and Tait, St. Venant, and Perry. 

4K, Federhofer, Z. Angew. Math. Mech. 25/27, 20 (1947). 
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Stress resultants and couples are given as follows 
: °F : 10F , 1 oF 2 ( oF) 
N = : > 4 — sailing” Joags a) 9 4 = e 
, or’ NV, r Or + r ow’ Nv ~ or\r 80 (3) 
aw , vdw y ew) (3 dw , 1 dw *w) 
f l or” *; or r 06°/’ Me D r or ar r 08 +? or’/’ 
(4) 
- 0 fi oe) 
Me = (1 v)D or a (2 00 
: aoV*w : aV*w " 
V, = -—D ie Ve = —D ie? (5) 
R,=V, +22 py, Mant %, 
(6) 
. 0 af 6 | dw ,- Ow 
R, = V,+- + N,z te le 


In addition to this there occur concentrated corner forces of magnitude 2M,, at corners 
formed by lines r = const. and @ = const. 

Let r = a and r = b be the inner and outer circular edges of the plate and let 6 = 0 
and @ = a@ be the radial edges of the plate (Fig. 1). We prescribe the following boundary 
conditions, assuming for the present that the couples 7 indicated in Fig. 1 are absent, 


r=a,b:; N,=WN, =R, = M,=0 (7) 
ab ab b ab 
[ Near =| Nudr= | Near = | Medr=0 
6=0,a;)] Redr+2(M,). — 2Mu)s = P (8) 
| po 
|| Ry dr + 2rMy»). — 2rM ys), = 0 
The boundary conditions at the edges 6 = const. are taken in such form that a semi- 


inverse procedure of solution becomes possible. 

3. Solution of the boundary value problem. The fact that the solution to be obtained 
should give stress resultants and couples which are independent of the polar angle @ 
and a transverse displacement w which is proportional to @ and knowledge of the corre- 
sponding solution in three-dimensional linear elasticity suggests the following form of 
w and F, 

w= ké, F = F(r), (9) 
Substitution of (9) in the differential equations (1) and (2), with p(r, @) = 0, shows that 
(1) is identically satisfied and (2) becomes 


2 


V’V'F = En®, . (10) 
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From (10), 
F = 4Ehk*[(Inr)? + A Inr + Br + Cr Inz], (11) 
where A, B, C are constants of integration and an unessential additive constant has 


been omitted. 
Stress resultants and stress couples are now 


N, = 1 pnw 2 my + 4 + 2B + C(1 + 2 In |, 


a t zn 2 Ree : + 2B 4+ C(3+2In |, 


N., = 0, Vy, = ©, VY, = 0, (12) 


M. = 0. M, => 0, Mo => (1 = v)D J ° 
r 
k - 


r 


R,=0, R, = —2(1—»D : + N, 


Determination of the constants of integration A, B, C and k is effected by means of the 
boundary conditions (7) and (8), some of which are seen to be satisfied identically. The 
relevant remaining conditions are the first of equations (7) and the third and fifth of 
Eqs. (8). We have from (7) 


2 + 2B+ C(1+2Ina) = -2m 

(13) 
A J, . > In_b 

z2 + 2B + C(I + 2Inb) = -2 F 


The third of Eqs. (8) may be transformed as follows if account is taken of the fact that 
F'(b) = F(a) = 0, 


b ab 


| rNedr = | rF" dr =F), — | F’ dr = F(b) — F(a) = 0. 
We have then 
A In (b/a) + B(b? — a’) + C(b’ In b — a’ Ina) = —(In Dy)’ + (In a)’ 


= —In (b/a) In ab. (14) 


Equation (13) and (14) serve to determine A, B and C. The remaining fifth of Eqs. (8) 
is to be used for the determination of the relation between force P and deflection k. 
We may write 


ob ab 


[ Nor dr = [ Fr" dr = (F’r"), + | F'r~?* dr = | N,r™ dr. 
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Therewith, we have from this fifth of Eqs. (8), 


(1 —- » DK! + | py| 1+2ne — ‘+ am + 4 (1, - 1) 





ab? 2a” 2b" a 


+ 2B In” + c(n? + (In b)? — (In a) = P, (15) 


After some transformations in which we make use of the following formulas 


a(t - 1) = 20m? + (nt - ne), 
a b 





a” b? a 
2B(b* — a’) = —C[b® — a? + 2(b’ In b — a’ Ina)] — 2 In ® , (16) 
J 2a*b* b\? b? — a’ b +a’, b 
(3 (in -]) —- ——_— = hb — +=; h-/, 
lb° — a a 2 a b° —a a 
a=R-e, b=R-+¢, E = 21+ )G. (17) 


Eq. (15) appears in the form 


[1% G0] -BBL-G)T 
tf + j fi R h _ 2Gch' ] R ‘ (18) 


where 


(19) 
For sufficiently small values of x, practically up to x ~ 1/4, we have 
. 16, “5 
iia) ~ — 2, Fs (20a 
re 45 —s ) 
For larger values of x we have the limiting relation 
lim f(z) = 1 (20b) 
rl 


Equations (18) and (20) indicate the way in which the nonlinear effect depends on 
deformation and dimensions of the plate. For wide plates for which 1 — z < 1 we have 
the usual result that for linearity the ratio of deflection per unit of circumferential 
angle to plate thickness must be small compared to unity. As the plate gets narrower 
larger and larger deflections lie within the linear range, the condition for linearity being 
that (k/h)(c/R)° be somewhat smaller than unity. 

A further question of interest concerns the magnitude of the stresses in the plate. 
As we have neglected the effect of transverse shear deformation we cannot determine the 
magnitude of the transverse shearing stresses. But we can determine, within the accuracy 
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of the results of plate theory the magnitude of the horizontal shearing stress r,, = 
6M,,/h’. We have from (12) 


6k 
te = (1 —vD ae (21a) 
and from this 
EB kh 
Tro — r?,max -_ _* (2 
a) Tr. x1 +») a Ib) 
In the non-linear range we have in addition a stress ¢, = N,/h. We limit ourselves 
here to a consideration of this stress along the edges r = a and r = b of the plate. We have 
from (12), taking account of the fact that N,(a) = N,(b) = 0, 
cate b/c 4b? I 
hi «8S, 4+ 8 ~ yet hh, 
8 la be —a b —a a 
(22) 
1/2 In b/a 4a’ b 
(bd) = = EWN — 4G + C2 —- a 
on re} Lb° ' be —a t | b* — a’ mn Al 
where C is given by (16). 
Equations (22) may be written in the alternate form 
; 1 Ek’ f.(c/R) 1E 1 f,(c/R) 
(a) = EEE eC a(t rire (2% 
Oe 2 GO — ey (0) = aR + c/R)?’ 3) 
where 
(l1—2)?, l+z2z - (l+2 a 
(x) = 1-3 l — CR*(1 — x)} 1 -—- + 
; b “i-s"* [1 _ "i= 2 
sie (1+2)?, 1+ . (l—a2 ite ' 
fz) = 1- > in 2 + CR(1 + 2) E —: == i (24) 
ae 1 l+2fl+2’ l+2 { Ls, we ts)) 
CR => ~ In Int | ji - 
l—2z\ 2 x) f° 


Of the two stresses o,(a) and o,(b) the larger one is the stress o,(a) at the inner edge of 
the plate. For the sake of comparison we may write 


ofa) 1+» (<) k - 
t44) @32 fo R/h (25) 


For sufficiently small values of c/R, practically for c/R < 1/4, we have 


c 4(c \’ cel : 
9 i D D 26¢ 
15) ie R~ 4 ean 
For larger values of c/R we have the limiting relation 
lim £(% ) = ] (26b) 
e/R-1 
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Comparison of the stress ratio (25) with the load deflection relation (18) indicates 
that non-linearity affects the stresses much more strongly than the deflection characteris- 
tics of the plate. As long as c/R < 1/4 we have for example that when o,(a)/r,,(a) = 1/4 
then non-linearity is responsible for a three-percent correction only to the stress dis- 
placement relation and when o,(a)/7,»(a) = 1/2 then this correction amounts to only 
about twelve percent. 

4. Interaction between pure twisting and pure bending. A further problem within 
the present context is the problem of the effect of a moment 7 about the line of action 
of the force P on the relation between force P and deflection k. Differential equations 
and boundary conditions are as before with one exception, which consists in replacing 
the boundary condition {3 rN, dr = 0 by the condition 


b 
[ rN, dr =T. (27) 


*a 


The solution for the stress function F now consists of two parts, one due to non-linearity 
in k and with T = 0, and the other without the non-linearity in k but with 7’ ¥ 0. Since 
the first part of the solution is given in the previous section we may now restrict atten- 
tion to the second part of the solution. We write for this second part 


F = T{A* Inr + Btr’ + C*r Inz}. (28) 
The boundary conditions N,(a) = N,(b) = 0 are satisfied by setting 


A*ta™? + 2B* + C*(1 + 2Ina) = 0, 





(29) 
A*h? + 2B* + C*#(1 42nd) =0. 
The condition {2 rN, dr = T becomes 
A* In (b/a) + B*(? — a’) + C*(b’ In b — a’ Ina) = —T. (30) 
From (29) and (30) follows 
eB a# Laced, pen (Ly PBR ene), 
ab a 2 bo —a (31) 








senna ID Qab_ yt 
ca al - (eam ; 


A relation between P and k follows again from f? R, dr — 2{[M,,]; = P, in the form 














ne POs el APB — Oe , ape bs om b| 

(1 — v)Dk a8 + ri 2 Ba + 2B* In = + C*(1 + In ab) In= = P, (32) 
Introduction of A*, B* and C* from (28) gives after some transformations 

. T .(c\| _ 3PR | (<) | 

i — A+» rp (5) = acer 1 ~ \r) J? (33) 
where 

Lae Lee[lte l+2_ | 
of 2x m 1-—2z 2x In l1-—z i 


(34) 








es eee) 
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Equation (33) shows that a positive moment 7’, that is a moment which tends to close 
an open ring sector, reduces the stiffness of the ring plate in so far as the effect of the 
axial force P is concerned. In contrast to this a negative moment 7’ increases the trans- 
verse stiffness of the ring plate. For sufficiently large positive 7’ we have instability 
in the sense that we may have k ~ 0 when P = 0. It may, however, be that instability 
occurs in other modes of deformation for values of 7’ which are lower than that given 
by (33). 
When c/R is sufficiently small, we have the approximation 


YY ] ] 
fi(x) ~ 3 e zQ rt (34a) 


In the range of applicability of (34a) we may further write 


—— 2 . 
T = 3 orch, (35) 


where o7 ~ —a,(a) ~ o,(b). With this we have 


‘—_ (< ) o7cR ee! 
173 — ww 2 7 2239 a . (36 
En /\R EW’ R™~4 » 
The present section may be concluded by listing the form of the relation between 

k and P which holds when both the effect of 7’ and non-linearity is considered. We have 


then 
3(1 +») , <\(k) 7 | (5) Z) _ 3PR' | 7 (<) | - 
E +" At n) ~ 2+ SAR) en} = aGeatl'—~\R) |: 8? 


We make the following further observation. Finite transverse deflections of the ring 
plate will be associated, in the absence of a moment 7’, with circumferential displace- 
ments parallel to the undeflected middle surface of the plate. Such displacements are 
also caused by the moment 7’. This means that if the conditions of load application are 
such that circumferential displacements are prevented 7 will have a definite value, 
proportional to k? and the non-linear correction term in the transverse load-deflection 
relation for this form of load application will differ from the result (18) which holds 
when 7’ = 0. 

The problem of the determination of this modified correction term will be considered 
in the last section of this work, in conjunction with the problem of the ring plate with 
initial deflection. 

5. Pure twisting and bending of an initially deflected plate. We now consider an 
initially deflected plate, or helicoidal shell, with middle-surface equation 


W = Koé. (38) 
We assume that K is sufficiently small for the shell to be considered shallow. It seems 
that for practical purposes we may admit values of K up to about (R — c)/z. 


The differential equations for shallow shells which take the place of the flat-plate 
equations (1) and (2) have been given by Marguerre.” The necessary changes consist 


5K, Marguerre, Proc. Vth Int. Congress Appl. Mech. Cambridge 1938, p. 98. 
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in replacing the operation D(w) on the right of (2) by D(W + w) — D(W) and in re- 
placing w on the right of (1) by W + w. 
We may take as before 


w= ké, F = F(r) (9) 
which reduces the two differential equations for w and F to the one equation 


70°F = Bh aKk x. (39) 


Equations (3), (4) and (5) which define stress resultants and couples remain un- 
changed while in equations (6) w must be replaced by W + w. 
The form of the boundary conditions (7), (8) and (27) remains unchanged and this 


means that we have now 
F = 1EhR(2Kk + k’)[(Inr)? + A Inr + Br + Cr Inz) 
+ T[|A* nr + Btr + C*r Inv, (40) 


where A, B and C and A*, B* and C* are given by (16) and (31), respectively. 
In determining the relation between force P and deflection k according to the fifth 
of Eqs. (8) it must be observed that R, of (12) is now changed to 
R, = —2(1 — »)D “ +N, Ath ; (41) 


A comparison of (40) and (41) with the corresponding earlier results shows that Eq. 
(37) for P as a function of k and T is changed into 


fy 4 80+» , (2) b+ 2m + K) _ (<.) Z| 
a f( 2 (+ SNR) ne 


4 R 
: c) 7 | _ 3PR/, _ A 
= K\ + n{&) | = 9Gch* [ (5 j; & 


As long as k «< K equation (42) reduces to the linear relation 


to 





| 3(1 +») .(c\ 2K?] _ c\KT  3PR’[, _ (5) | 
A 7 oe (5) | (1 +d) ER ~ PE rR) |: 9) 


In view of the nature of the distribution of stresses ¢, which give rise to the terms account- 
ing for the initial deflection of the plate the term with K’ can be obtained only by taking 
account of plate action and is consequently not incorporated in the classical theory of 





curved beams. 

We note further that while for the initially undeflected plate the non-linear correc- 
tion term varies as the cube of the deflection we have, for the initially deflected plate, 
correction terms varying both as the square and as the cube of the additional deflection 
caused by the forces P. 

6. Consideration of circumferential displacement. We now consider in addition to 
the transverse displacement w radial and circumferential displacement components 
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u and v. Components of finite strain for the middle surface of the shallow shell are of the 


form 
_ u , 1) (AW + #)/ _ (ary | 
“= 431 or ar) 





_%,l% 1 aw + u)\ _ (20) 

ced r -f- r 00 + Dy  ( r 00 r 6 ; (44) 
=i Mire (?) | + wv) (W + w) _ aw or 

ye 7 00 ar \r ar r 30 ar 7 00)" 


For the present problems we have W = K06, w = ké, du/00 = O and y,, = 0. From 
this it follows that 


v= > or, (45) 


where w is the relative angular displacement of the ends of one winding of the plate. 
The expressions for e, and ¢, reduce to 


wo, gate ge EEE a 
leading to a compatibility relation of the form 
dr € w Kk + 3k’ ‘a 
~~ -**s° - (47) 
We introduce the relations 
Ehe = N, — wWN,, Eh, = N, — wNo (48) 


and express NV, and N, in terms of the stress function F by means of (3), taking account 
of the fact that here dF /0@ = 0. This leads to the result that 

3 2m a 2. 1],2 
| dF dF 1 ar | _ 2 _ Kk+}e (49) 


2r r 


Eh 


tg SS 

dr® dr rdr 
Into (49) we introduce F from Eq. (40), and this gives us after some cancellations the 
relation 


2r 


1 q* 2\ 07 4T we 4 
5 (2Kk + W)C +5, C* = 5, (50) 


where C and C* are defined by (16) and (31). 

If we wish the axial force-displacement relation under the assumption of vanishing 
circumferential displacement we find from (50), with w = 0, as appropriate value of the 
couple 7’, 

T (2Kk + k*)C 2Kk + k’ E +a’, b | b , 
<= — OTe | ee | eth - 2 (51 
Ch 8C* 8 bo —-a “s a |* a pa) 
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Introduction of this value of 7 into Eq. (42) gives after some transformations the follow- 
ing relation between k and P, 


a , 31+) (5) (K + (2K + B] _ 8PR [1 _ (:) | . (52) 











4 R h* ~ 2Gch* R 
where 
fz) = 1 - (15% mite) (53) 
2x l—2z 
and 


f(x) ~ : z*, ts : ; lim f,(x) = 1. (54) 
3 2 = 

It may be seen that in the range c/R < 1/4 the effect of non-linearity and of initial 
deflection is much more pronounced when w = 0 than it is when 7 = 0, a factor (4/15)- 
(c/R)* which occurs when T = 0 being replaced by a factor (c/R)* when w = 0. On the 
other hand, when c/R is sufficiently near to unity the effect of non-linearity and of pre- 
twist is the same for JT = 0 and for w = 0. 


TaBLE I. Values of functions occurring in load displacement relations and stress ratio o4/Trg . 














Zz hi fe fs ts 

0 0 0 2 0 

0.1 0.000036 0.0159 29.08 0.0135 
0.2 | 0.000574 0.0465 14.32 0.0530 
0.3 0.00328 0.0836 8.97 0.1186 
0.4 | 0.01024 0.1691 6.14 0.208 
0.5 0.0265 0.255 4.31 0.321 
0.6 0.0600 0.359 2.98 0.453 
0.7 0.1248 0.483 1.94 0.600 
0.8 0.2484 0.631 1.106 0.756 
0.9 0.4923 0.809 0.427 0.903 
1.0 1.0 1.0 0 | 1.0 
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—NOTES— 
A SUFFICIENT CONDITION FOR AN INFINITE DISCRETE SPECTRUM* 


By C. R. PUTNAM (Purdue University) 


1. In the differential equation 

x’ + fix = 0, (1) 
let f = f(t) denote a real-valued, continuous function on the half-line 0 S t < &. 
Both necessary and sufficient conditions in order that the equation (1) be oscillatory, 
so that every solution of (1) possesses an infinity of zeros on 0 S ¢ < © clustering at 
+o, are known; see, for instance, [10], [5]. The present note will deal primarily with 
the problem of obtaining a sufficient criterion in order that (1) be oscillatory, in the 
particular case that f(t) satisfies the limit relation 


f(j) — 0, as [—> @, (2) 
The following will be shown: 


(I) If f(t) satisfies (2), then the differential equation (1) is oscillatory whenever the 


inequality 
»T 8 

lim sup a( tim sup s| lim sup | f() dt + | f(t) |(t —y)/(S — y) |'™* at}) > (3) 
h-+0 \ S-@ T-< “8S “0 


holds for every fixed number y 2 0. 


Obviously, the inequality is satisfied in case the innermost “lim sup” is +. (It 
should be pointed out here that if (2) holds and if limr..f7 f(é)dt fails to exist either 
as a finite limit or as — ~, then (1) is surely oscillatory; [3], p. 389. Cf. also [11] and [2].) 
It is noteworthy that the criterion furnished by (I) remains valid if the assumption 
(2) is replaced by certain other conditions; cf. the remark at the end of section 2. 

The result (I) will have an implication concerning the discrete spectrum of the 
boundary value problem on 0 S t < © ([7]) determined by the differential equation 


a’ ++ f(d)r = 0 (4) 
and the boundary condition 


z(0) cosa + x’(0) sna = 0, OSa<rn, (5) 


in the particular case that (2) holds. For every fixed a, the relation (2) implies that 
the half-line X = O belongs to the spectrum; [2]. Moreover, only a discrete spectrum 
(isolated eigenvalues) can exist for \ < 0 and there exists a finite number of eigenvalues 
X < 0, or an infinity of such eigenvalues clustering only at \ = 0, according as (1) is 
not or is oscillatory; [7], p. 252. Accordingly, (I) provides a sufficient condition that 
(4) and (5), in the case (2), determine a boundary value problem with an infinity of 
negative eigenvalues clustering at \ = 0. (It should be pointed out that (4) may have 





*Received Jan. 28, 1953. 
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a non-trivial solution z(#0) of class L*[0, ~) for a positive value \ even when (2) holds; 
cf. the constructions of [8], pp. 394-395 and [9], pp. 268-269.) 

In section 3 below, a specific application of the criterion (I) to the spectral problem 
mentioned above will be made. The boundary value problem to be considered will be 
of the type arising in the quantum mechanical treatment of the two particle problem. 
In this case, there is a singularity in the coefficient function of the differential equation 
at the origin but the nature of the problem remains essentially identical with that 
considered in connection with (4) above; cf. [4], pp. 154, 163. 

That the assertion of (I) can become false if the strict inequality > of (3) is relaxed 
to = can easily be shown by an example. In fact, if f(t) = Ct-*, where, say,0 <t < @, 
it is readily seen that (3) reduces to the inequality C > }. However, if C = }, so that f(é) 
becomes 3¢~”, then (1) possesses the non-oscillatory solution x = ¢'”, for0 <t < @. 
(Needless to say, the criterion (I) actually requires the continuity of the function f(é) 
only for large values ¢ and the fundamental interval 0 S t < @ may be replaced by 
any half-line 7 S t < o.) 

2. The proof of (I) will depend upon an application of an oscillation criterion ob- 
tained in [5]. (A somewhat similar application was made in [6] in the case that f(t) was 
periodic.) It was shown in [5] that (1) is oscillatory if and only if, for every T = 0, 
there exists on 7 < ¢ < @ a continuous function z = z(t) satisfying z(7) = 0 and 
possessing a piecewise continuous derivative 2z’(¢) such that each of the integrals 


=) PT ) 


[ ca, f xa, and | saat 
JT JT JT 
is finite, while 


a wo 


| (w? — fa’) dt <0. (6) 


ee i 


It will now be shown that relations (3) and (2) imply (6), and so, (I) will follow. 
For every positive number 7 and for every number n > 3, consider numbers T;, < 
T, < T, such that T = T, < T,, and define the function x = z(t) on T, St < © as 


follows: 


w= T;)", for T,s¢taT; ; 

\(T2 = T;)", for T: Ss t Ss T; ; 
z(t) = 

(7. - Tye T;) ‘(T. — %), for T. S2¢8 7. ; 

0, for 1224 <X ®, 


If T, — T, = Aand T, — T; = B, direct calculation readily shows that the requirement 
(6) reduces to 
aTs aT 


A fdt+A™ | . f(t — T,)*" dt 


oT, ’T; 


Ts 
4+ AB”? [ f(T, — 0? dt > n°(Qn — 1)? + AB". 
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Suppose that 7, (hence A) is determined and then, for an arbitrary positive number 
e, choose B so as to satisfy AB™* < e«. Next choose 7’; (hence 7',) so large that 


| Z a 
AB? | #(T, — O° dt| <e. (7) 
j “Ts | ° 


That this can be done is clear from (2) and from the fact that | 7, — ¢ |? B-? S 1 when 
T; Sts T, (= T; + B). It now follows that (6) is surely satisfied if 
aTs Ts 
Af sat+ arm | ft—7,)" dt > n'Qn — 1) + 2%. (8) 
“Ts “T; 
Actually it will be shown that (8) holds for (certain) large values A, 7, (that is, for 
certain A, T, — ~) as a consequence of (3). 
Let L(S) be defined by 


»T 


L(S) = lim sup | f(b dt, (9) 


T32 7S 


and for a fixed value T, , choose 7 so large that 
Ts 
[  fat> LT.) — 73’. 
JT; 


It is clear that for a fixed value T, , AT," > 1 as T, — ~, and so, for certain large A, 
T,, the relation (8) will hold if 7’, is fixed and 


aTs 
T.L(T:) + T, | f(t — T,)/A]™ dt > n°(Qn — 1)" + 2 
7T; 


holds for certain large 7, . However, this last relation will certainly hold for certain 
large T, (7, fixed) if 


S 


lim sup | sx(s) +58 f\|(@—T7T)/(S — T,) |” at| > n°(2n — 1)" + Qe. (10) 


pi 
S-@ 0 


Since ¢ can be chosen arbitrarily small, the 2e appearing on the right side of (10) may 
be deleted. If now h = 2n — 1, one sees that (6) holds for every T (= 7, = 7) if the relation 


a 


h lim sup s| 115) + | fi) |(t-y/S -—y |'* at| > (1 + h)’/4 (11) 
holds for some h > 0 (where possibly h = h(y)) and is an arbitrary non-negative number. 
However, relations (3) and (9) clearly imply that (11), for every y 2 0, is valid for some 
positive h and thus the assertion (I) is proved. 

Remark. The assumption (2) was used in order to obtain (7). It is clear though from 
the proof given above that (3) will imply (6) if 


.T2+B jf, Js 0 
AB* | = f-()(T, +B-tdt|<~6 fo = 


| JTs 


} 


0, f>0o0 
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can be obtained for A, B fixed and for (certain) large 7’, . Thus, the criterion furnished 
by (I) will be valid if, for instance, the restriction (2) is replaced by any one of the 
three assumptions (i) f(f) 2 0, or even (ii) f(t) ~0ast— o, or 


(iii) [ i at | ce. 





3. In order to obtain an application of (I), consider the radial portion of the separated 
form of the quantum mechanical wave equation of the two particle problem, namely, 


R”’ + cA — Vir) — UL + 1)/er*)R = 0; 


cf. [4], p. 150. Here, c, \ and / are constants and the prime denotes differentiation with 
respect to r. It will be assumed that V(r) — 0 asr — o; ef. [4], p. 152. It is clear from 
an earlier remark that, as far as concerns the application of (I), only the continuity of 
the function V(r) for large r, say for 1 S r <~, is needed. (As is customary, it will 
be assumed that the singularity of V(r) at r = 0 is of a suitably restricted type; cf., 
e.g., [4], pp. 152, 163.) It follows from (I) and the calculation of section 1 that the equation 


R” + e—V(r) — Ul + 1)/er’)R = 0 (12) 


is oscillatory if, for all y 2 0, 


T Ss 
lim sup n( tim sup sf tim sup [ —cV(r) dr+ [ —cV(r) | (@r — y)/(S — y) |*" ar }) 
S-@ T-@ “8S “1 


A~+0 


>i+Ul+1). (13) 


Thus, if the last inequality is satisfied, there exists an infinity of negative energy levels 
clustering at X = 0. In the case of the hydrogen atom, V(r) = kr-*, where k < 0, so 
that the bracketed portion of the left side of (13) is + © and hence (12) is oscillatory 
for all values /; cf., e.g., [4], pp. 157 ff. In many cases however, the potential V(r) appears 
to be unknown (cf. [4], p. 156 and [1], pp. 30 ff.) and the relation (13) offers a property 
of V(r) guaranteeing the existence of an infinite (negative) discrete spectrum. 
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THE TORSION AND STRETCHING OF SPIRAL RODS (II)* 


By H. OKUBO (Tohoku University, Sendai, Japan) 


The torsion and stretching problems of spiral rods were discussed in a preceding 
paper.’ There, the equations of equilibrium were expressed in terms of displacements 
that were independent of the position of the section perpendicular to the axis of a spiral 
rod. The differential equations for the displacements were integrated for the particular 
sase where the helix angle was small, and the corresponding displacements and stresses 
were obtained. In the calculations, however, the displacements were preliminarily 
assumed in special forms, and consequently solution was valid for some special problems. 
In the previous paper, the displacements for the stretching problem were assumed in 
forms that reduce to those for a uniform tension in the limit case when the helix angle 
approaches zero. But when a spiral rod with axis that does not pass the centroid of the 
cross section, is pulled axially, the displacements must be in forms that reduce to those 
for a uniform tension combined with a uniform bending moment in the limit case when 
the helix angle approaches zero. Hence, the validity of the previous solution was re- 
stricted to the problem for a spiral rod with axis through the centroid of the cross section. 

As in the preceding paper, we take the axis of helix as the axis of z, and denote the 
displacements in x’, y’, z directions by wu’, v’ and w respectively, in which 2x’ and y’ are 
the axes perpendicular to each other and fixed to a section of the rod perpendicular to z. 
We take for the displacements the expressions 


a, B’ 
ul’ =u, — yx’ a(x” — y”) — ay’z + re (1 — cos kz — kzsin kz), 
ra 1 | 
, | 
‘ . are , 8 : i 
vy =v, — yy’ — y't'y’ + ar’ — z2 (sin kz — kz cos kz), (1) 


w= w, + —(x' sin ke — y’ cos kz + y’) + Bz, 
where u; , v; , w, are the functions of x’, y’ and are independent of z, k is a constant which 


specifies the helix angle, a, 8, 8’ are arbitrary constants, and 


9 


- 


l 
= { oe = ( — , = ——., 
Y 1 — pp, Y 5 (1 — vB, ama CO 


From (1) we have the cubical dilatation 


Ou, Ov, 


le —— kD, W,) + pB’x’ + ps, 
Ox oy 
where 
ne ee Q 
D,=y ry x ay” 


*Received January 13, 1953. 
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The equations of equilibrium for this case can be expressed in the forms 


>,” + p{ Viu, + k* D, (uy) = 2k? D.(0;) a 5 ky (a ~. y”’) — B’} = 0, 





re] ‘ 7 
ay + p{ Vw. + k’D,(v,) + 2k’?D,(u,) + k*y’z’y’} = 0, (2) 
kDA(Ad) — p{Viw, + k’D,(w,) + k’w, — kp’y’} = 0,) 
where 
; 3” 3’ 0 7) rf) 0 
Vi = >> 73, £ ) aig oo — 2 SAY 2 SO}. 
ei Ox” + Oy ~ and +1 (y Ox . Oy Nu Ox " 3) 


The differential equations (2) are independent of z. Solving the simultaneous equa- 
tions for u, , v, and w, , we can find the displacements from (1). The displacements 
for a straight rod are readily obtained from (1), by taking the limit case when k ap- 


proaches zero; thus 


, , 
“= —yx — aye — . (x? — y’) — c “| 
v= —yy + arz — y'xry, [ (3) 


| 
w= w, + Be + B’2z, 
where uv, , ¥; are assumed to vanish when k approaches zero. 
Assume that w, and a in (3) also vanish when k approaches zero; the corresponding 
stresses become 
X,= Y, =X, =X, = Y, = 0, Z, = u(3 — p(B + B’z). (4) 
This is the solution for a straight rod submitted to simple tension combined with a 
uniform bending moment. If we assume that 8 and @’, instead of w, and a, vanish when 
k approaches zero, the corresponding stresses become 
, ; ’ ? Ow, : Ow, 
X,= Y,=2Z2,=Z, =0 xX,= (2 — ay) Y,= (42 + az). 5 
. 2 a ’ be ox y ’ KB Oy ( ) 
This is the solution for the torsion problem of a straight rod. 
Consider now a spiral rod of small k, pulled by a pair of axial forces. Assuming that 
a, wW, are small quantities of the order k and u, , v, are of the order k*, and neglecting the 
small quanties of the higher order, the equations of equilibrium (2) can be written as 


follows: 


) 


0 , ‘ - . 
ee + pViu + 5 Pky (a — y’”) — pp’ = 0, 


(6) 


—) 


0A Galast 
gy? + PV + pkiy'2’y’ = 





= 


Viw, — 2kp’y’ = 
Let us take for w, the expression 


o, = i(fs = fa) + 1kp’y’®, (7) 
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where f; is an arbitrary function, f; is a function conjugate with fs , and 


, 


f- = fd), fs = f.(0), r= 2’ + iy’, ¢t=2' — iy’. 
The displacement w, satisfies the third Eq. (6). Substituting this expression of w, into 
the first and second equations (6), the equations for u, and v, become 


. Ou, av, au, | 
1 ») a > + ae -j- | 
(+7 Ox ~ dx’ dy’ - P Oy ~ | 
is | 

= k(fs + fs + fs’ + Sf3') — 5 pk y'2’ (4 —p+py)k py” (8) 
| 
Ov, ‘ Oy (1+ ») 222 | 
ox * Ox OY ~* OY 


° 44 ra writ pr l, 2) 7.2 , 
= ik(fs — fs + Sf3’ — Sf3) —54+PD-P )k’B’x y’. | 
Integrating the differential equations (8), we find 
u= b +. Fa = xr'(fe oa fe) +k | f3¢ d¢ +k [ Fer dé 


aoe 


+ Fe {3 — pa! — 64 + p — pay” + (6 — p— Py" },¢ = 9) 
v, = if, —f) + ix(f. - fo) +i1+ 2n)\ | fed — | 4 a), 


where f; , f, are arbitrary functions of ¢. The corresponding stresses become 


XL = Wt HAL GALI + lle + fd) + MES + FD | 


ro 1g _ ahatarety’™ 
+ 53 (8 — PL + pik'a’s g (9 perc), 


Y2, = —Qufit fit 2(f2 + fd) + (2+ pi(fe + fo) 
4 (3 — pl — p)k’p’x"” + (1 — py'k's'2'y"}, 
Z, = —24\( — p)(fo + fo) + k(ES, + 53) ; (1 — p38 — p)k*e'x - | (10) 


1 ! 
3 (7 + 2p — p)k’B’z'y” — 5 (3 — p(B+ Bix’) ¢ 


( ' o k - 
X’, = 2 fl — fl tri —f) + (1+ Die — fe) + 5 (SS3 — Sf) | 


+ 8 4 + p — pk B'x"y’ — 3 (6 — p— PK BY, 


X! = pli(fs — f3) — ay’ + hy’2’y’}, 
( 


| | 
, } py 7 14,02 ky’ /2 72 | 
Yy’ = — WSs + fi — az’ — kp'y” += (@" —-y Mf | 
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Take for the bounding curve of the section, the expression 
F(z’, y’) = 0. (11) 


The conditions for the lateral surface of the rod to be free from external forces are 


x oF + X!, oF — kDAF)X! = 0, 
Ox Oy 
x, OF 4 vs, 2% _ eam! =0,, (12) 
Ox oy 
” a 
x2 + 7: > — uk(3 — p)\(6 + B’z’)D.AF) = 0. 





Consider a spiral rod stretched by a pair of axial forces P, and imagine a small portion 
of the rod cut by two parallel planes perpendicular to the axis of rod, as shown in Fig. 1. 


|Z 





Fie. 1. 
The equilibrium condition of the surface tractions for this portion is 


I Z, az' dy’ = P, 


I/ Z,x' dx' dy’ = 0, (13) 


we 





I (y’X! — x’ Y!) dx’ dy’ = 0, 


where the integrals are taken over the cross section. 
The arbitrary functions f; , f. and f; are determined so as to satisfy the boundary 


condition (12), and the constants a, 8 and #’ are obtained from the condition (13). 
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Substituting the expressions of X/ , Y! in (10) into the third Eq. (12), it becomes 


] _ 
fa + (3 — p)ks} 5 (58) 


bd 


L a — apna’ @ — 1 5 na’ © xy"? =0 
— — 3p)kp’ —— — = (5 — p)k z'y’’ = 0, 
12 ( Pp h ds 4 ( ) P h ds lt 7] 

where ds is the element of are of the bounding curve of the cross section. It follows that 
the equation 


fs +fs = 4 fa + (3 — p)kB}(x” + yy”) 


, Bie 
+ ‘ (11 — 3p)k@’x" + Fr (5 — p)kp’x’y”’ + const., (14) 


holds on the bounding curve, from which f; is determined. For the convenience of further 
calculations, we shall rewrite the stresses X/. , Y/. and X/. in (11) in the expressions as 
ps) , Y 


X!. = : 5 + Quk(efs + of) + - (3 — pl + p)k’p’x” — ra — p)k’p’x'y”, 
ou “a : 


where 
x= -—g-A-—teta, & = Wht), = Wh. 


Substituting these expressions for stresses into the first and second Eqs. (12), by 


virtue of (14), we obtain 


, dP 


] Ox Ox - d . - > 
a e. i. ¢ 4- uh - [F3(¢) — F3(¢)] + uk’B a 


ds \@x 0 y’ 





1 . 2n/& d e to’ a ‘ \ + Pe d ( te 2 — 
= 9 (3 . p) uk B'¢ ds “88? — { (3 — p) mk Buy ds 88 + 36) = 0, 
(16) 
1 (dx .0x | l Fae rr 207 db 
= 7 9 ? a “ } = pl en [F 2 ( ¢) 7 fk 3\ ¢) + pk 8B ds 
I wom ic i ak a - 
— ~(3 — p)uk’B’t — (62) + 7 (3 — puk’B’x"y’ =| (85 + $) = 0, 
2 ds 4 ds 


where 


F; = 2¢fs, 
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@= — i (11 + p)(3 — px” + 6(1 — p)’x”?y” — (9 — 2p + py") 


+ 5 (3 — pC + px” — (7 + py" Iz’y’. 


From (16), f; and f, are obtained. 
As an example of the procedure, consider an elliptic spiral rod whose cross section is 


(x’ — c/a? + y”/b = 1, (17) 


as shown in Fig. 2. 











Transform now the elliptic section in the ¢-plane into a unit circle in the ¢-plane by 
the equation 


t=ct+a’'tt+ bt’, (18) 


where a’ + b’ = a,a’ — b’ = b. 


Take for the function f, , the expression 
fs=e(t+st') +e(? +80’) +¢(f +s’), (19) 


where S = b’/a’. The unknown coefficients c, , c. and c; are readily obtained from the 
boundary condition (14). 

Remembering the condition that g, and ¢g, are analytic at any point of the section in 
the ¢-plane, we take for the functions the expressions 


a 4 
g = > AM! + 8°t’), o = DB + s8't"), (20) 
n=1 n=1 


where A, and B, are real constants. These unknown constants are obtained from the 
boundary condition (16), and the other unknown constants a, 8 and @’ are finally obtained 
from the conditions (13). 

In the case of stretching, the predominating stress is the normal stress Z, and the 
shearing stresses X/ and Y! follow it, but the latter are smaller quantities of the order k. 
The other stresses are of the order k* and are very small quantities when k is small. 
When the section of the spiral rod is a circle of radius unity, then a = 6 = 1, and it 
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follows a’ = 1, b’ = 0, s = 0. The main stress Z, for the circular section, referred to the 
polar coordinates with the pole at the center of the circle, is 


Z, = w(3 — p(B + Bc) — 21 — p)B, — 4uke,ec 


+ [ue’(3 — p) — 4(1 — p)B, — 4uk(2c.c + ¢,) |r cos 0 


— 2(/3(1 — p)B; + Quk(3c3c + 2c,)]r° cos 26 — 4[2(1 — p)B, + 3ykes)r° cos 30 = (21) 


l 


aan 12 uk’B’(1 — p)(3 — p)\(c +r cos 6)" 


1 


+ = pk’?3'(7 + 2p — p (ce +1 cos 6)r’ sin’ 6. 
1 I 


4 


From the third Eq. (13), 


we have 


a = 3(1 — p)kp'e. (22) 


Hence, a twist almost proportional to k arises when a spiral rod with a circular section 


Fic. 


3. 
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The shearing stress along the bounding circle, when k = 0 25. 








1954] H. OKUBO 495 


is pulled axially.” From the remaining Eqs. (13), the other unknown constants 6 and 
B’ are obtained. 

If r be the shearing stress along the bounding circle of the section, it becomes 
} =—y{[2c, — 4k6’c*(1 — p) — 1pkB’] cos 6 + 4c, cos 26 + (6c; + 4k6’) cos 36}. (23) 


5 
The shearing stress along the periphery has been calculated from (23) for various values 
of c, assuming the Poisson’s ratio and k to be 0.3 and 0.25, respectively, and is shown 
in Fig. 3. As is seen from the figure, the shearing stress becomes large at both ends of two 
diameters parallel to the coordinate axes (z’, y’), and attains its maximum value at the 
outer end of the diameter on the x’-axis. The distribution of the normal tension Z, on 
the axis of x’, obtained from (21), is given in Fig. 4. For the sake of comparison, the 


~ 
r 














\ C45 
‘\ 





Fic. 4, The distribution of Z, on the z’-axis, when k = 0.25. 


corresponding distribution of Z, for a straight rod is also shown by dotted lines in the 
same figure. 
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23’ is a function of k and c, but it remains almost constant for the variation of k, when k is small. 
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ON THE GAPS IN THE SPECTRUM OF THE HILL EQUATION* 
By C. R. PUTNAM (Purdue University) 


1. Let f = f(t) be a real-valued, continuous, periodic function of period 1, so that 
fi) ~ x: c, exp (27 int), (c.., = 6), (1) 


and consider the Hill equation 
ze’ +A + fz = 0, (Areal; ’ = d/dt). (2) 


It is known that (if f # 0) there exists a sequence of closed intervals J,:\, < 4 < 
(region of stability), where \, < AY < xs, and k = 1, 2, --- , with the property that 
(2) has some solution zx # 0 which is bounded on — ~ < ¢ <= if and only if \ belongs 
to the closed set S = >> J, ; cf. [7], p. 14. The complementary set of S consists of a half- 
line — © < <, and the sequence of open intervals J, :\” < \ < Aye, k = 1, 2, --- 
In several recent papers, various lower bounds for fhe value X, , the least point of the set 
S, in terms of the Fourier coefficients c, of f(t), have been obtained; [11], [5], [3]. The 
present note will be devoted to the problem of obtaining estimates (upper bounds) of 
the lengths \,., — \* of the “gaps” J, of the set S in terms of these Fourier coefficients. 
It follows from [4], p. 613, that the length of every gap J, is surely not greater than 


lim sup f(t) — lim inf f() < 4 y 2 a © (3) 
t@ ta n=1 
— oo, for these gaps are known; [2]. In fact, since 


In addition, asymptotic estimates, as 2’ 
< o, the lengths \,., — d* of the intervals J, 


f(t) is uniformly continuous on 0 < ¢ 
tend to zero as \z4; ~~; loc. cit., p. 850. Furthermore, additional regularity conditions 
on f(é) result in more refined estimates. It should be pointed out here that the investiga- 


tions of [2] related to singular boundary value problems ([8]) on the half-line 0 < ¢ < © 
determined by (2) and a linear, homogeneous boundary condition at t = 0, and were not 
confined to the special case that f(t) be periodic. 

Let m(A), for —© < \ <~o, be defined to be the distance from \ to the set S con- 


sidered above, so that 


mr) = g.l.b. |A — pl, win S. (4) 
It will be shown in section 2 below that m(A) satisfies the inequality 
m(r) < 2 ree iat provided A> -G. (5) 
n=1 


e ~ . k . 
As a consequence of (4) and (5), one readily sees that the lengths \,., — ° of the gaps 


J, satisfy 


7 2 é 1/2 ! 1 
Aisi ~A S 2( 2 > Cr ) : provided 5 Aner + M) > -—a. (6) 


_ 


It will remain undecided whether (6) actually must hold for all gaps J, , so that the 
first inequality of (6) would hold without the proviso of the second inequality. In any 


*Received Aprilji3, 1953. 
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case, it is readily seen that the estimate of (6), when it applies, is an improvement over 
that of (3), namely 4 = » 1 & fh. 
In this connection, it should be pointed out that Kato [3], by an adaptation of a 


relation used by Wintner [11], has obtained the inequality 


1 @ 
1 2-4 — (2) >» | Cp J 
n=1 


for the least point A, of the set S. (Wintner had previously shown that A, > —¢) — 2: 
>=, | ¢, |?.) Consequently, it is easily seen that the first inequality of (6) is surely valid 
for all gaps J, if, for instance, the inequality 
| wo @o : 1/2 
(1) Lila l < (2 par P) 
holds. (If one normalizes f so that its mean value is zero, hence c, = 0, this last inequality 
is equivalent to f, f° dt < 256). 

Before proceeding to the proof of (5), it can be noted that the first inequality of 
(5) surely becomes false if the restriction \ > —c, is dropped. In fact, if f(t) = co , so 
that (2) becomes the differential equation of the harmonic oscillator, then >. *-, | ¢, |? = 0, 
and (5) yields the known result that m(A) = 0 for’ > —c) . However, m(A) > 0 for 
\ < — ¢, since S is the half-line —cg < AX <@&. 

2. The proof of (5) will depend upon certain results obtained in [6]. Let g, (8), go(t), --- , 
denote a sequence of functions possessing continuous second derivatives on 0 < t <@, 
satisfying 

g.(0) = g(0) = 0, (7) 
and such that g,(¢) — 0 uniformly on every finite ¢-interval [0, 7]. Then, if g, and L(g,) 
(where L(x) = x’’ + fx) are of class L’ [0, ~), the inequality 


m*(d) lim inf | g, dt < lim inf | (L(g,) + Agn)” dt (8) 


holds. This follows readily by a method analogous to that given in [6], p. 580. (It is to be 
noted that the set S considered above is identical with the invariant spectrum (Weyl 
[8], p. 251) associated with the differential equation (2); [9], [1]. Moreover, the investi- 
gations of [6] related to the Wey] theory of singular boundary value problems, alluded to 
in section 1.) 

Next, let » > 0, and let g, = y,h, where h = sin (ut) or h = cos (ut), and the 
Yn = y(t) are functions possessing continuous second derivatives on 0 < t < @. In addi- 
tion, suppose that y,(0) = y/(0) = 0, so that (7) certainly holds, and that y, and L(y,) 
belong to L*(0, ©). Finally, suppose that the y, are such that the “lim inf” appearing 
on the left side of the inequality (8) can be replaced by “lim” for both h = sin (n't) 
and h = cos (ut). 

It follows from (8) that 


m*(r) lim | y2h? dt < lim inf | (yo +A — ut fly jh + 2y%h’)’ dt. 
If now the y, satisfy 


| yo dt— 0, | yi”? dt— 0, (n — co), 


“0 “0 
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it is seen that 


m’(r) lim | y,h* dt < lim inf | (A — uw t+ f)’ysh? dt. (9) 
Since (9) holds for both functions h, addition of the two corresponding inequality rela- 
tions yields, in view of the fact that lim inf A + lim inf B < lim inf (A + B), the in- 
equality 


m*(d) lim | y, dt < lim inf | (X\— ut f)’y dt. (10) 

Let T > O and define the function Y,(t) on 0 < t <@ so that the graph of Y(t) 

on 0 < ¢ < T consists of.three line segments joining, in order, the four points (0, 0), 

(1, T~*), (T — 1, T™*), and (T, 0). On T < t <~@, let Y,(t) = O. It is clear that the 

corners of this function can be smoothed out so as to obtain a function y7(t) satisfying 

the conditions imposed upon the y, above. Furthermore, it is clear that if y, = yr, , 

where T = T,— © as n —©, one can arrange that the functions y, be such as to make 
(10) imply 


m*(X) < lim inf S™’ | (XY — w+ f) dt, (u > 0). (11) 
It is clear that the inequality u > O in (11), and not merely wu > 0, can be allowed.) 
1 } : 


Now suppose that \ > —c, and choose » > 0 so that \ — wp = —c,. Then (11), (1), 
and the Parseval relation yield 


so that the relation (5) is now proved. 
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AN ADDITION TO PORITSKY’S SOLUTIONS OF A 
DIFFERENTIAL EQUATION OF TORSION* 


By J. C. WILHOIT, JR. (Stanford University) 
In the theory of circular shafts of variable cross section, the only stress components 
which do not vanish are the shearing stresses, 


_ Fay 


= = 3 
. r Or’ 


Py = -5 or (1) 


P 
r Ox 


the stress function y is a function of r and x alone, and the coordinate system is cylindrical 
(r, x, 6). Equations (1) satisfy the equilibrium equations and will satisfy compatibility if 


a(l ov) a (! ov) se 
or (3 or + ax \r* dx] 0 (2a) 
or 

vy say , Hy _ 

or” or Or + ox” 0. (2b) 


In problems relating to cones or spheres, it will be convenient to express Eq. (2a) in 
spherical coordinates R, w, ¢, as 


fs. 3. mb. Oh. 5. 24 ” 
OR E (1 — yp)’ OR + Ou E (i-p) an) 0, (3) 


where F is the radius in spherical coordinates and 4» = cos w = 2/R. 
The product solutions of (3) are: 


y= (1 —p)’Pi.WR", (4a) 
¥ = (1 — w’)’PL(w)R", (4b) 
¥=(1—2)'Q.WR”, (4c) 
¥ = (1 — 2')’QU WR", (4d) 


where Pn and Qn are the Legendre functions of the first and second kind. 
The solutions (4) do not give functions y varying as R, and only Eq. (4d) gives a 


solution varying as R*(n = —1). These missing solutions are: 
¥=(1+4)R vy = R's, (5a) 
y=(1+4)R v= Ru. (5b) 


Solutions (4) are given by equations (79) and (80) of reference [1]. The solutions 
(5a) and (5b), however, do not agree with those given in Eqs. (107) and (108) of reference 
[1]. 

Although one is usually interested in torsion of shafts with stress free surfaces, it is 


*Received April 15, 1953. 
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interesting to consider the type of loading which results from Eqs. (1) if this requirement 
is relaxed. For a cone, consideration of the equilibrium of a small volume element at the 
surface of the cone will show that the surface stress consists of rings of shear with a 


magnitude 


_ldy| 1 dy | ) 
Poe = GR (6) 


R(1 Ese u’) dR Sener ag 





pweonst. 

Equation (6) indicates that for a given shear distribution on the surface of the cone, 
the boundary condition is on (d¥/dR),-..n.... rather than on y. It is evident that although 
Eqs. (4) and (5) give a complete set of functions varying as any positive or negative 
power of FR, this set does not include any function which has a derivative behaving as 
1/R. The two functions which have a derivative behaving as 1/F are, omitting arbitrary 


multipliers. 


y; = E of. (3u — ad ) log 1 2s +- log (1 —#) +. log r|, 


a 4 l1— up 2 
atte 6 is ite. Gee)... » « 8), | 7 
.=|[5 9 3 i, 5 log (1 pw) + 3 log R |. (7) 


To obtain these functions the following procedure was used: 
1. Equation (3) is multiplied by R* and differentiated with respect to R giving 


ca Se i Se a ove - . 
OR E OR (4: (1 — yp) ve} + Ou FE ~~ ae 0. () 


2. The desired solutions are those with ~, behaving as 1/R. Thus zg must be of the 


form. 
f(u) ’ 
= _ 9 
Wr R . { ) 
Equation (8) becomes an ordinary differential equation in f(u), and shows that 
/ uw 1 
=| Alu — > = ) 
vr | a() “) +22 (10) 
and therefore 
y= | (. _ e) + B| log R + g(u), (11, 
where A and B are constants of integration 
3. With 
¥, = BlogR+ gly), yy, = A(n = H) log R + go(u), (12) 


substitution into Eq. (3) will give ordinary differential equations for the unknown 
functions g,(u) and g.(u). These may be integrated, giving 


+ 


glu) = 3 By? + 2B | (1 — w?) log + dy, 
4 & l—up 
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ou) = Au — 5 f (1s?) log (1 — 22) du. (13) 


If the integration indicated in Eq. (13) is performed and the values of g,(u) and g2(u) 
used in Eq. (12), the expressions given in Eq. (7) are obtained. 
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BOOK REVIEWS 


Description of a magnetic drum calculator. By The Staff of the Computation Laboratory. 
Harvard University Press, Cambridge, 1952. 318 pp. $8.00 


This book is one of a series put out by the staff of the Computation Laboratory of Harvard Univer- 
sity and describes the Mark III calculator. This machine was completed in March of 1950 and was then 
moved to the Naval Proving Ground at Dahlgren, Virginia. 

The book itself is a detailed description of this machine and is of principal interest to those persons 
who are immediately associated with the machine. It combines both an engineering and a mathematical 
description of the device. 

The text is extremely well illustrated both with photographs and with schematics of the principal 
organs of the machine. 

To illustrate the coding of problems for the machine there is a chapter which contains among other 
things the programming for four illustrative examples. 

The text is undoubtedly an invaluable aid to those immediately concerned with the operation and 


programming of problems for the Mark III calculator. 
H. H. Goupstine 


Introduction to the theory of plasticity for engineers. By Oscar Hoffman and George Sachs. 
McGraw-Hill Book Company, Inc., New York, Toronto, London, 1953. xib + 276 


pp. $6.50. 


As one would expect on the basis of the authors’ well deserved reputation their book contains an 
excellent exposition of the technological application of plasticity. Major emphasis is placed on the approxi- 
mate solutions to problems of rolling, extruding, drawing, etc. in which the material is assumed to be 
ideally plastic and the true three dimensional character of the flow is not taken into account. Appreciable 
space is also devoted to the elementary classical problems of the thick-walled shell and tube and to the 
rotating cylinder and disk. Although the text opens with a discussion of stress and strain tensors and 
considers a stress space, the discussion of stress-strain relations and experimental data is brief and is 
essentially confined to the maximum shear stress and octahedral shear stress criteria. The extensive 
modern literature on stress-strain relations in the plastic range is, in the main, ignored. Also, except for a 
short section on two-dimensional plastic flow problems, little of the classical mathematical theory of 
plasticity is treated. No mention is made of plastic waves nor of the theorems or applications of limit 
analysis and design. The latter would be especially useful in evaluating some of the results obtained in 


the approximate solutions which are treated. 
D. C. DrucKER 
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Tensor calculus. By Barry Spain. Oliver and Boyd, Edinburgh and London, Interscience 
Publishers, Ine., New York, 1953. viii + 125 pp. $1.55. 


In this readable little book, the author manages to cover a surprisingly large amount of material— 
tensor algebra and differentiation, and introductions to differential geometry, elasticity, and relativity. 
This is accomplished by using a very concise style of writing in which some aspects of a subject are pre- 
sented only formally. The book should be of value to both mathematicians and engineers. 

The following material is covered in the text. Tensors are defined by their transformation laws, and 
addition, multiplication, contraction, and the quotient law are discussed. The metric tensor and the 
principal directions of a second order symmetric tensor in n-dimensional Riemannian space are studied. 
After formally introducing the Christoffel symbols, the covariant derivative is defined and its properties 
are determined. Geodesics, parallelism, and the curvature tensor of n-dimensional Riemannian space are 
discussed. In addition, a brief survey of three-dimensional differential geometry is given (Frenet formulas 
for curves in space, the normal vector of a surface, the second fundamental tensor of surface, etc.). In the 
section on elasticity, orthogonal transformations, rotations, infinitesimal strain, the compatibility rela- 
tions, the stress tensor, Hooke’s law, and the equilibrium relations are studied. A brief introduction is 
given to curvilinear coordinates and isotropic tensors. Finally, an introduction is given to the special 
theory of relativity and then the general theory, including the Schwartzschild line-element and the 


Einstein and De Sitter universes, is discussed. 
N. CoBuRN 


The theory of homogeneous turbulence. By G. K. Batchelor. The University Press, Cam- 
bridge, 1953. xi + 197 pp. $5.00. 


This book gives an excellent account of the modern developments in the statistical theory of homo- 
geneous turbulence. The author has purposely chosen to omit all work requiring a Lagrangian description, 
since the methods used are very much different from those in the Eulerian description. Problems like 
diffusion are therefore omitted in the present book. It is also clear from the title of the book that a treat- 
ment of shear flow is not to be found in the present volume. The omission of these two important aspects 
of turbulence might discourage some people chiefly interested in applications. However, it is generally 
agreed among workers in turbulence that a treatment of the statistical theory of homogeneous turbulence 
would lead to the best understanding of our established knowledge and basic concepts. 

Chapter I gives a general description of the problems involved in the statistical theory of turbulence, 
and a brief history of the subject, beginning with Taylor’s work of 1935. Chapter II describes the methods 
used in the mathematical representation of the field of turbulence, including a discussion of the method of 
taking averages and the relation between correlation and spectral tensors. Chapter III applies these 
concepts to the treatment of specific velocity correlations. 

Discussion of the dynamical aspects of turbulence begins with Chapter IV, where a set of linear 
problems are collected and treated in some detail. In Chapter V, the general aspects of the decay of 
homogeneous turbulence are treated, including a discussion of the final period of deday, where lineariza- 
tion is again justified. 

The universal equilibrium theory of A. N. Kolmogoroff is described in Chapter VI. General dis- 
cussions of the energy transfer is included here. The specific spectrum proposed by Townsend is, however, 
postponed to Chapter VII where Heisenberg’s form of the energy spectrum (as calculated by Chan- 
drasekhar) is also discussed. The concept of self-preservation during the process of decay is treated in 
Chapter VII. The author presents some limited experimental information before introducing the general 
concept, although the general ideas were conceived and used theoretically for the prediction of law of de- 
cay before these detailed experimental information were available. It is quite likely that the general 
concepts will hold in other cases not yet examined experimentally. 

The final chapter gives a discussion of the probability distribution of the velocity fluctuation and the 
consequences drawn from the approximately Gaussian character of the joint velocity distribution at 
several points. The book ends with an excellent bibliography. 
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Advances in applied mechanics. Volume III. By Richard von Mises and Theodore von 
Karman. Academic Press Inc., New York, 1953. x + 324 pp. $9.00. 


This volume contains the following articles. 


G. F. Carrier: Boundary layer problems in applied mechanics. The term “boundary layer problem,” 
borrowed from hydrodynamics, is here applied to problems leading to non-dimensional differential equa- 
tions in which the coefficient of the highest derivative is small compared to the coefficients of lower 
derivatives. The solution of such a problem may often be obtained by matching a “boundary layer 
solution” to an “interior solution’’. Several examples are given that illustrate this technique. 


O. Zaldastani: The one-dimensional isentropic fluid flow. The discussion is restricted to analytical solutions 
of the linear equations obtained by applying a Legendre transformation to the equation of continuity and 
the momentum equation for one-dimensional isentropic flow. Particular attention is given to the possible 
occurrence of singularities. The procedure of unfolding the characteristic plane that the author ascribes 
to Ludford in this connection is actually due to 8. Christianovich (Mat. Sbornik 1, 511, 1936). Interaction 
of simple waves, expansion of a monatomic gas into a vacuum, and flow in a closed tube are discussed as 
examples. 

F. N. Frenkiel: Turbulent mean diffusion: mean concentration distribution in a flow field of homogeneous 
turbulence. On account of the importance of turbulent diffusion to chemical engineers, meteorologists, 
and other readers not expected to be familiar with the general theory of turbulence, the article begins 
with a concise discussion of the statistical description of a turbulent field. The main portion of the article 
is concerned with diffusion in a field of homogeneous isotropic turbulence in which the decay of turbulence 
is neglected. Diffusion from a point source and diffusion from a line source are discussed. An extension 
of the results to a case of non-isotropic turbulence is discussed in the appendix. 


H.F. Ludloff: On aerodynamics of blast. The article treats the diffraction of weak shocks around wedges of 
arbitrary angle, the pressure and density fields behind blasts advancing over arbitrary flat surfaces, and 
the head-on encounter of a shock with a wall that is almost parallel to the shock front. Most space is 
devoted to the last two problems and the discussion is based on previously published work by the author, 
C. 8. Gardner, and L. Ting. 


G. Guderley: On the presence of shocks in mixed subsonic-supersonic flow patterns. A supersonic region that 
is embedded in a subsonic flow field will, in general, contain shocks. The first part of the article discusses 
this phenomenon from the standpoint of the mathematician and the second, from that of the physicist. 
The mathematical problem is attacked by studying the change that a small deformation of the boundaries 


produces in a transonic potential flow. It is shown that the resulting boundary value problem does not, in 
general, admit a continuous solution. The particular solutions that indicate the impossibility of a potential 
flow are then discussed physically. They represent wavy flow patterns that may lead to physically unac- 


ceptable overlapping portions of the flow field. 

L. Rosenhead: Vortex systems in wakes. This, by far the shortest article in the volume, reviews the litera- 
ture on vortex systems in wakes for the range of Reynolds numbers up to about 2500 and deals mainly 
with papers of a theoretical nature. The author concludes that “‘most of what has been done is interesting, 
suggestive, and qualitatively correct, but if greater precision is needed, it is necessary to reinvestigate 
in fine detail much of what has already been done.” 

H. Geiringer: Some recent results in the theory of an ideal plastic body. Most of this article is devoted to a 
thorough mathematical discussion of plane plastic flow. Arbitrary isotropic yield functions and plastic 
potentials are considered, and it is not assumed that the yield function and the plastic potential are 
identical. Much of the material may strike the superficial reader as well-known, but a closer inspection 
reveals a wealth of results that render more precise, amplify, or generalize, classical theorems. Unfortu- 
nately, the space available for this review does not allow detailed comments on this great number of 
small but valuable contributions to the mathematical theory of plasticity. 


A. I. Bellin: Non-autonomous systems. The term “non-autonomous”’ is used for a system of one degree of 
freedom whose acceleration a depends explicitly on the time ¢ and not only on the displacement z and the 
velocity v. Particular attention is given to the case where a(z,v,t) is a periodic function of t, and the 
behavior of the system is discussed with the aid of a three-dimensional phase space with the coordinates 
and ¢. Criteria for the stability of periodic solutions are presented. Important special cases are dis- 
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cussed in detail. 


























